P.C. Matthews 


Ç £ i | s 


Springer Undergraduate Mathematics Series 


Advisory Board 


P.J. Cameron Queen Mary and Westfield College 
M.A.J. Chaplain University of Dundee 

K. Erdmann Oxford University 

L.C.G. Rogers University of Cambridge 

E. Süli Oxford University 

J.F. Toland University of Bath 


Other books in this series 


A First Course in Discrete Mathematics I. Anderson 

Analytic Methods for Partial Differential Equations G. Evans, J. Blackledge, P. Yardley 

Applied Geometry for Computer Graphics and CAD, Second Edition D. Marsh 

Basic Linear Algebra, Second Edition T.S. Blyth and E.F. Robertson 

Basic Stochastic Processes Z. Brzeźniak and T. Zastawniak 

Complex Analysis J.M. Howie 

Elementary Differential Geometry A. Pressley 

Elementary Number Theory G.A. Jones and J.M. Jones 

Elements of Abstract Analysis M. Ó Searcóid 

Elements of Logic via Numbers and Sets D.L. Johnson 

Essential Mathematical Biology N.F. Britton 

Fields, Flows and Waves: An Introduction to Continuum Models D.F, Parker 

Further Linear Algebra T.S. Blyth and E.F. Robertson 

Geometry R. Fenn 

Groups, Rings and Fields D.A.R. Wallace 

Hyperbolic Geometry J.W. Anderson 

Information and Coding Theory G.A. Jones and J.M. Jones 

Introduction to Laplace Transforms and Fourier Series P.P.G. Dyke 

Introduction to Ring Theory P.M. Cohn 

Introductory Mathematics: Algebra and Analysis G. Smith 

Linear Functional Analysis B.P. Rynne and M.A. Youngson 

Mathematics for Finance: An Introduction to Financial Engineering M. Capiriksi and 
T. Zastawniak 

Matrix Groups: An Introduction to Lie Group Theory A. Baker 

Measure, Integral and Probability, Second Edition M. Capiriksi and E. Kopp 

Multivariate Calculus and Geometry, Second Edition S, Dineen 

Numerical Methods for Partial Differential Equations G. Evans, J. Blackledge, P.Yardley 

Probability Models J.Haigh 

Real Analysis J.M. Howie 

Sets, Logic and Categories P. Cameron 

Special Relativity N.M.J. Woodhouse 

Symmetries D.L. Johnson 

Topics in Group Theory G. Smith and O. Tabachnikova 

Vector Calculus P.C. Matthews 


P.C. Matthews 


Vector Calculus 


With 63 Figures 


G Springer 


Paul C. Mathews, PhD 
School of Mathematical Sciences, University of Nottingham, Urijversity Park, 
Nottingham, NG7 2RD, UK 


Cover illustration elements reproduced by kind permission of: 

Aptech Systems, Inc., Publishers of the GAUSS Mathematical and Statistical System, 23804 S.E. Kent-Kangley Road, 

Maple Valley, W A 98038, USA. Tel.: (206) 432-7855 Fax (206) 432-7832 email:info@aptech.com URL:www.aptech.com 

American Statistical Association: Chance Vol 8 No 1, 1995 article by KS and KW Heiner ‘Tree Rings of the Northern Shawangunks’ 
page 32 fig 2 Springer-Verlag: Mathematica in Education and Research Vol 4 Issue 3 1995 article by Roman E Maeder, Beatrice 
Amrhein and Oliver Gloor ‘Illustrated Mathematics: Visualization of Mathematical Objects” page 9 fig 11, originally published as a 
CD ROM ‘Illustrated Mathematics’by TELOS: ISBN 0-387-14222-3, german edition by Birkhauser: ISBN 3-7643-5100-4. 
Mathematica in Education and Research Vol 4 Issue 3 1995 article by Richard J Gaylord and Kazume Nishidate “Traffic Engineering 
with Cellular Automata” page 35 fig 2. Mathematica in Education and Research Vol 5 Issue 2 1996 article by Michael Trott’The 
Implicitization of a Trefoil Knot’ page 14. 

Mathematica in Education and Research Vol 5 Issue 2 1996 article by Lee de Cola ‘Coins, Trees, Bars and Bells: Simulation of the 
Binomial Process’ page 19 fig 3. Mathematica in Education and Research Vol 5 Issue 2 1996 article by Richard Gaylord and Kazume 
Nishidate ‘Contagious Spreading’ page 33 fig 1. Mathematica in Education and Research Vol 5 Issue 2 1996 article by Joe Buhler and 
Stan Wagon ‘Secrets of the Madelung Constant’ page 50 fig 1. 


British Library Cataloguing in Publication Data 
Mathews, P.C. 

Vector calculus.- ( Springer undergraduate mathematics series) 

1. Vector analysis 2. Calculus of tensors 

LTitle 

§15.6°3 

ISBN 3540761802 
Library of Congress Cataloging-in-Publication Data 
Matthews, P.C. (Paul Charles), 1962- 

Vector calculus/ P.C. Matthews. 
cm— (Springer undergraduate mathematics series) 

Includes index. 

ISBN 3-540-76180-2 (pbk. : acid-free paper) 

1. Vector analysis. I. Title. II. Series 
QA433.M38 1998 97-4192 
§15’.63—dc21 
Apart from any fair dealing for the purposes of research or private study, or criticism or review, as permitted under the Copyright, 
Designs and Patents Act 1988, this publication may only be reproduced, stored or transmitted, in any form or by any means, with 
the prior permission in writing of the publishers, or in the case of reprographic reproduction in accordance with the terms of licences 
issued by the Copyright Licensing Agency.Enquiries concerning reproduction outside those terms should be sent to the publishers. 


Springer Undergraduate Mathematics Series ISSN 1615-2085 
ISBN3-540-76180-2 Springer-Verlag London Berlin Heidelberg 
Springer Science+Business Media 

springeronline.com 


The use of registered names, trademarks etc. in this publication does not imply, even in the absence of a specific 
statement, that such names are exempt from the relevant laws and regulations and therefore free for general use. 


The publisher makes no representation, express or implied, with regard to the accuracy of the information 
contained in this book and cannot accept any legal responsibility or liability for any errors or ommissions that may 
be made. 


Mathews: Vector Calculus 
© 1998, Springer-Verlag London Limited 
All rights reserved. No part of this publication may be reproduced, stored in any 


electronic or mechanical form, including photocopy, recording or otherwise, 
without the prior written permission of the publisher. 


First Indian Reprint 2005 
Second Indian Reprint 2008 


ISBN 978-81-8128-295-8 


This editicn is manufactured in India for sale in India, Pakistan, Bangladesh, 
Nepal and Sri Lanka and any other country as authorized by the publisher. 


This edition is published by Springer (India) Private Limited, 
A part of Springer Science+Business Media, Registered Office: 906-907, 
Akash Deep Building, Barakhamba Road, New Delhi - 110 001, India. 


Printed in India by Rashtriya Printers, Delhi. 


Preface 


Vector calculus is the fundamental language of mathematical physics. It pro- 
vides a way to describe physical quantities in three-dimensional space and the 
way in which these quantities vary. Many topics in the physical sciences can 
be analysed mathematically using the techniques of vector calculus. These top- 
ics include fluid dynamics, solid mechanics and electromagnetism, all of which 
involve a description of vector and scalar quantities in three dimensions. 

This book assumes no previous knowledge of vectors. However, it is assumed 
that the reader has a knowledge of basic calculus, including differentiation, 
integration and partial differentiation. Some knowledge of linear algebra is also 
required, particularly the concepts of matrices and determinants. 

The book is designed to be self-contained, so that it is suitable for a pro- 
gramme of individual study. Each of the eight chapters introduces a new topic, 
and to facilitate understanding of the material, frequent reference is made to 
physical applications. The physical nature of the subject is clarified with over 
sixty diagrams, which provide an important aid to the comprehension of the 
new concepts. Following the introduction of each new topic, worked examples 
are provided. It is essential that these are studied carefully, so that a full un- 
derstanding is developed before moving ahead. Like much of mathematics, each 
section of the book is built on the foundations laid in the earlier sections and 
chapters. In addition to the worked examples, a section of exercises is included 
at the middle and at the end of each chapter. Solutions to all the exercises are 
given at the back of the book, but the student is encouraged to attempt all 
of the exercises before looking up the answers! At the end of each chapter, a 
one-page summary is given, listing the most essential points of the chapter. 
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The first chapter covers the basic concepts of vectors and scalars, the ways 
in which vectors can be multiplied together and some of the applications of 
vectors to physics and geometry. 

Chapter 2 defines the ways in which vector and scalar quantities can be 
integrated, covering line integrals, surface integrals and volume integrals. Again, 
these are illustrated with physical applications. 

Techniques for differentiating vectors and scalars are given in Chapter 3, 
which forms the essential core of the subject of vector calculus. The key concepts 
of gradient, divergence and curl are defined, which provide the basis for the 
following chapters. 

Chapter 4 introduces a new and powerful notation, suffix notation, for ma- 
nipulating complicated vector expressions. Quantities that run to several lines 
using conventional vector notation can be written extremely compactly using 
suffix notation. One of the main reasons for writing this book is that there 
are very few other books that make full use of suffix notation, although it is 
commonly used in undergraduate mathematics courses. 

Two important theorems, the divergence theorem and Stokes’s theorem, are 
covered in Chapter 5. These help to tie the subject together, by providing links 
between the different forms of integrals from Chapter 2 and the derivatives of 
vectors from Chapter 3. 

Chapter 6 covers the general theory of orthogonal curvilinear coordinate 
systems and describes the two most important examples, cylindrical polar co- 
ordinates and spherical polar coordinates. 

Chapter 7 introduces a more rigorous, mathematical definition of vectors 
and scalars, which is based on the way in which they transform when the 
coordinate system is rotated. This definition is extended to a more general 
class of objects known as tensors. Some physical examples of tensors are given 
to aid the understanding of what can be a difficult concept to grasp. 

The final chapter gives a brief overview of some of the applications of the 
subject, including the flow of heat within a body, the mechanics of solids and 
fluids and electromagnetism. 
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Vector Algebra 


1.1 Vectors and scalars 


This book is concerned with the mathematical description of physical quanti- 
ties. These physical quantities include vectors and scalars, which are defined 
below. 


1.1.1 Definition of a vector and a scalar 


A vector is a physical quantity which has both magnitude and direction. There 
are many examples of such quantities, including velocity, force and electric field. 
A scalar is a physical quantity which has magnitude only. Examples of scalars 
include mass, temperature and pressure. „~ 

In this book, vectors will be written in bold italic type (for example, u is a 
vector) while scalar quantities will be written in plain italic type (for example, 
a is a scalar). There are two other commonly used ways of denoting vectors 
which are more convenient when writing by hand: an arrow over the symbol 
(W) or a line under the symbol (u). 

Vectors can be represented diagrammatically by a line with an arrow at the 
end, as shown in Figure 1.1. The length of the line shows the magnitude of the 
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vector and the arrow indicates its direction. If the vector has magnitude one, 
it is said to be a unit vector. Two vectors are said to be equal if they have the 
same magnitude and the same direction. 


Fig. 1.1. Representation of a vector. 


Example 1.1 


Classify the following quantities according to whether they are vectors or 
scalars: energy, electric charge, electric current. 

Energy and electric charge are scalars since there is no direction associated 
with them. Electric current is a vector because it flows in a particular direction. 


1.1.2 Addition of vectors 


Two vector quantities can be added together by the ‘triangle rule’ as shown 
in Figure 1.2. The vector a + 6 is obtained by drawing the vector a and then 
drawing the vector b starting from the arrow at the end of a. 


a+b 


Fig. 1.2. Addition of vectors. 


The vector —a is defined as the vector with magnitude equal to that of a 
but pointing in the opposite direction. 
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By adding a and —a we obtain the zero vector, 0. This has magnitude zero 
and so does not have a direction; nevertheless it is sensible to regard 0 as a 
vector. 


1.1.3 Components of a vector 


Vectors are often written using a Cartesian coordinate system with axes z, y, z. 
Such a system is usually assumed to be right-handed, which means that a screw 
rotated from the z-axis to the y-axis would move in the direction of the z-axis. 
Alternatively, if the thumb of the right hand points in the z direction and the 
first finger in the y direction, then the second finger points in the z direction. 

Suppose that a vector a is drawn in a Cartesian coordinate system and 
extends from the point (21, y1, 21) to the point (T2, Y2, z2), as shown in 
Figure 1.3. Then the components of the vector are defined to be the three 
numbers a; = T2 — £1, @2 = y2 — yı and az = 22 — zı. The vector can then be 
written in the form a = (a1, @2, a3). 


a ee an Z,) 


(X Yez)? 


a 


Fig. 1.3. The components of the vector a are (£2 — £1, Y2 — Y1, Z2 — 21). 


By introducing three unit vectors e;, e2 and e3, which point along the 
coordinate axes z, y and z respectively, the vector can also be written in the 
form a = a, €; + a2€2 + age3. Using this form, the sum of the two vectors a 
and b is a +b = a ,e) + ane, + a3e3 + b1€1 + b2€2 + b3e3 = (a, + b; ey + (a2 + 
by )e2 + (a3 + b3)ẹ3. It follows that vectors can be added simply by adčing their 
components, so that the vector equation c = a + b is equivalent to the three 
equations cı = a; + b1, C2 = a2 + be, cz = ag + b3. 
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The magnitude of the vector is written |a|. It can be deduced from Pythago- 
ras’s theorem that the magnitude of the vector can be written in terms of its 
components as |a| = ya? + aż + a. 

The position of a point in space (x,y, z) defines a vector which points from 
the origin of the coordinate system to the point (x,y,z). This vector is called 
the position vector of the point, and is usually denoted by the symbol r, with 
components given by r = (2, y, z). 


Example 1.2 


The vectors a and b are defined by a = (1,1,1), b = (1, 2,2). Find the magni- 
tudes of a and b, and find the vectors a + b and a — b. 

The magnitude of the vector a is |a| = VI? + 12 + 12 = V3. The magnitude 
of b is |b] = V1? + 2? + 2? = 3. The vector a +b is (1,1,1) + (1,2, 2) = (2,3, 3) 
and a — b = (0, —1,—1). 


< 
1.2 Dot product 


The dot product or scalar product of two vectors is a scalar quantity. It is written 
a:b and is defined as the product of the magnitudes of the two vectors and 
the cosine of the angle between them: 


a-b= |a||b| cos 8. (1.1) 
A number of properties of the dot product follow from this definition: 


e The dot product is commutative, i.e. a-b=b.a. 

e If the two vectors a and b are perpendicular (orthogonal) then a - b = 0. 

e Conversely, if a -b = 0 then either the two vectors a and b are perpendicular 
or one of the vectors is the zero vector. 

ea-a=|al’. 

e Since the quantity |b| cos @ represents the component of the vector b in the 
direction of the vector a, the scalar a -b can be thought of as the magnitude 
of a multiplied by the component of b in the direction of a (see Figure 1.4). 

e The dot product is distributive over addition, i.e. a-(b+c)=a-b+a-c. 
This follows geometrically from the fact that the component of 6 + c in the 
direction of a is the same as the component of b in the direction of a plus 
the component of c in the direction of a (see Figure 1.5). 
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0 


<=——— |blcos 8 ———> a 


Fig. 1.4. The component of b in the direction of a is |b| cos 0. 


Fig. 1.5. Geometrical demonstration that the dot product is distributive over addi- 
tion. 
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A formula for the dot product a- 6 in terms of the components of the two 
vectors a and b can be derived from the above properties. Considering first 
the unit vectors €1, €2 and e3, it follows from the fact that these vectors have 
magnitude 1 and are orthogonal to each other that 


é€,°e; = 1, e2 -6&2 = 1, €3:e3 = 1, € -e€5 = 0, €2:e3 = 0, €3 -€i = 0. 


The dot product of a and b is therefore 


a-b = (aie: + a2€2 + azez) . (bie: + b2e3 + b3e3) 
= abe) - €; + a2b2€2 - €2 + a3b3€3 - €3 
= a,b) + azb: + a3b3. (1.2) 
Example 1.3 


Find the dot product of the vectors (1,1,2) and (2, 3, 2). 
(1,1,2) - (2,3,2)=1x2+1x34+2x2=9. 


Example 1.4 


For what value of c are the vectors (c, 1,1) and (—1, 2,0) perpendicular? 
They are perpendicular when their dot product is zero. The dot product is 
—c + 2 + 0 so the vectors are perpendicular if c = 2. 


Example 1.5 


Show that a triangle inscribed in a circle is right-angled if one of the sides of 
the triangle is a diameter of the circle. 


a> 


a a 


Fig. 1.6. Geometrical construction to show that a is a right angle. 


Introduce two vectors a and b as shown in Figure 1.6. Since these two 
vectors are both along radii of the circle they are of equal magnitude. The two 
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sides c and d of the triangle are then given by c = a+b and d = a —b. The dot 
product of these two vectors is c-d = (a@+b)-(a—b) = |a|*—a-b+b-a—|b|? = 0. 
Since the dot product is zero the vectors are perpendicular, so the angle a is a 
right angle. This is just one of many geometrical results that can be obtained 
using vector methods. 


1.2.1 Applications of the dot product 


Work done against a force 


Suppose that a constant force F acts on a body and that the body is moved a 
distance d. Then the work done against the force is given by the magnitude of 
the force times the distance moved in the direction opposite to the force; this 
is simply — F - d (Figure 1.7). 


d 


Fig. 1.7. The work done against a force F when an object is moved a distance d is 
-F -d. 


Equation of a plane 


Consider a two-dimensional plane in three-dimensional space (Figure 1.8). Let 
r be the position vector of any point in the plane, and let a be a vector 
perpendicular to the plane. The condition for a point with position vector r to 
lie in the plane is that the component of r in the direction of a is equal to the 
perpendicular distance p from the origin to the plane. The general form of the 
equation of a plane is therefore 


r:a = constant. 


An alternative way to write this is in terms of components. Writing r = (2, y, z) 
and a = (a1, 4@2,43), the equation of a plane becomes 


aiT + azy + azz = constant. (1.3) 
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O 


Fig. 1.8. The equation of a plane is r - a = constant. 


EXERCISES 


1.1 Classify the following quantities according to whether they are vec- 
tors or scalars: density, magnetic field strength, power, momentum, 
angular momentum, acceleration. 

1.2 If a = (2,0,3) and b = (1,0, —1), find |a|, |b], a+b, a—b and a-b. 
What is the angle between the vectors a and b? 

1.3 If u = (1,2,2) and v = (—6,2,3), find the component of u in the 
direction of v and the component of v in the direction of u. 

1.4 Find the equation of the plane that is perpendicular to the vector 
(1,1,—1) and passes through the point z = 1, y = 2, z = 1. 

1.5 Use vector methods to show that the diagonals of a rhombus are 
perpendicular. 

1.6 What is the angle between any two diagonals of a cube? 

1.7 Use vectors to show that for any triangle, the three lines drawn from 
each vertex to the midpoint of the opposite side all pass through the 
same point. 
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1.3 Cross product 


The cross product or vector product of two vectors is a vector quantity, written 
a x b. Since it is a vector, its definition must specify both its magnitude and 
direction. The magnitude of a x b is |a||b| sin 0, where @ is the angle between the 
two vectors a and b. The direction of a x bis perpendicular to both a and bina 
right-handed sense, i.e. a right-handed screw rotated from a towards b moves in 
the direction of a x b (Figure 1.9). We may therefore write a x b = |a||b| sin 8 u, 
where u is a unit vector perpendicular to a and b in a right-handed sense. 


axb 


a 


Fig. 1.9. The cross product of a and b is perpendicular to a and b, in a right-handed 
sense. 


The cross product has the following properties: 


e The cross product is not commutative. Because of the right-hand rule, a x b 
and b x a point in opposite directions, so a x b = —b x a. 

e If the two vectors a and b are parallel then a x b = 0. 

eaxa—QO. 

e The magnitude of the cross product of a and b is the area of the parallelogram 
made by the two vectors a and b (Figure 1.10). Similarly the area of the 
triangle made by a and b is |a x b|/2. 

e The cross product of a and b only depends on the component of b perpen- 
dicular to a. This is apparent from Figure 1.10 since the component of b 
perpendicular to a is |b| sin 8. 

e The cross product is distributive over addition, i.e. ax (b+c) =axb+axce. 
This is demonstrated geometrically in Figure 1.11, where the vector a points 
into the page. The vectors b, c and b + c do not necessarily lie in the page, 
but from the previous point the cross products of these vectors with a only 
depend on their projections onto the page. The effect of taking the cross 
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ae E 


a 


Fig. 1.10. The area of the parallelogram is the length of its base, |a|, multiplied by 
its height, |b| sin 8. 


product with a on any vector is to project it onto the page, rotate through 
n /2 clockwise and then multiply by |a|. Thus the triangle made by the vectors 
b, c and b + c becomes rotated and scaled as in Figure 1.11 but remains a 
triangle. 


axe 


a X (b+c) 


Fig. 1.11. Geometrical demonstration that the cross product is distributive over 
addition. The vector a points into the page. 


A formula for the cross product a x b in terms of the components of the two 
vectors a and b can be derived in a similar manner to that carried out for the 
dot product. Consider first e; x e2. Since these two vectors have magnitude 
l and are perpendicular, sin? = 1 and the magnitude of e x ez is 1. The 
direction of e; x ez is perpendicular to both e; and ez in a right-handed sense, 
SO €1 X €z = @3. 

It follows that the unit vectors e,, e2 and e3 obey 


€, Xe, = 0, e2 xez = 0, 23 xe3 = 0, €i X€ = €3, €2 X €3 = Cj, €3 X€] = Ep. 


The cross product of a and b is therefore 
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axb = (ae; + a2€2 + a3e3) x (b1e1 + b2€2 + b3e3) 
= @)b2e1 X €2 + a1b3 e) X €3 + dob; eo X €] 
+aob3 €2 X €3 + a3b1 €3 X €e; + a3b2 €3 x ez 
= (azb — a3b2)eı + (a3bı — aıbz3)ez2 + (aıbz — a2b1)ez3. (1.4) 


This can also be written as the determinant of a 3 x 3 matrix as follows: 


€; €2 €3 
axb= âi Qq a3 
bi bz b 


Example 1.6 


Find the cross product of the vectors (1,3,0) and (2, —1, 1). 
(1,3,0) x (2,-1,1) = (8 — 0,0 — 1, -1 — 6) = (3, -1, -7). 


Example 1.7 


Find a unit vector which is perpendicular to both (1,0,1) and (0,1, 1). 

A perpendicular vector is (1,0,1) x (0,1,1) = (—1,—1,1). To make this a 
unit vector we must divide by its magnitude, which is V3, so the unit vector 
perpendicular to (1,0, 1) and (0,1, 1) is (—1, —1, 1)/v3. 


Example 1.8 


What is the area of the triangle which has its vertices at the points P = (1, 1,1), 
Q = (2,3,3) and R = (4,1, 2)? 

First construct two vectors that make up two sides of the triangle. The 
vector from P to Q is a = (1,2,2) and the vector from P to R is b = (3,0, 1). 
The cross product of these vectors is a x b = (2,5, —6). The area of the triangle 
is then |a x b|/2 = 65/2 = 4.03. 


1.3.1 Applications of the cross product 


Solid body rotation 


Suppose that a solid body is rotating steadily about an axis. What is the 
velocity vector of a point within the body? 

Consider a body rotating with angular velocity (2 (this means that in a 
time t the body rotates through an angle 92t radians). Since there is a rotation 
axis, a vector §2 can be defined, with magnitude |§2| = N and directed along 
the rotation axis. Since this vector could point in either direction, the following 
form of the right-hand rule is used to define the direction of 92: a screw rotating 
in the same direction as the body moves in the direction of §2. Alternatively, if 
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the fingers of the right hand point in the direction of the rotation, the thumb of 
the right hand points in the direction of $2. This means that for a body which 
is rotating to the right, 92 points upwards (Figure 1.12). 


Fig. 1.12. Motion of a rotating body. 


Now consider the motion of a point at a position vector r, which makes 
an angle 0 with the rotation axis. The speed at which this point moves is Nd, 
where d is the perpendicular distance from the point to the rotation axis. Since 
d = |r|sin@ (Figure 1.12), the speed of motion is v = |r| sin @. Note that this 
is equal to |92 x r|. Now consider the direction of the motion. In Figure 1.12, 
where both f2 and r lie in the plane of the page, the direction of motion is 
into the page, perpendicular to both £2 and r and so in the direction of 2 xr. 
Therefore the velocity vector of the point at r is 


v=N xr, (1.5) 


since this vector has both the correct magnitude and the correct direction. 


Equation of a straight line 


The equation of a straight line can be written in terms of the cross product 
as follows. Suppose that a is the position vector of a particular fixed point on 
the line, and that u is a vector pointing along the line (Figure 1.13). Then any 
point r on the line can be reached from the origin by travelling first along the 
vector a onto the line and then some multiple of the vector u along the line: 


r =a + àu, (1.6) 
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where À is a parameter. This is referred to as the parametric form of the equa- 
tion of a line. 


O 


Fig. 1.13. The equation of a line is r = a + Au. 


To obtain a form of (1.6) that does not involve the parameter A, the term 
involving the vector u must be eliminated. This can be done by taking the 
cross product of (1.6) with u. This gives r x u = a x u. Since the vector a x u 
is a constant, it can be relabelled b, giving the second form for the equation of 
a straight line: 

rxu=b. (1.7) 


Physical applications of the cross product 


There are many physical quantities that are defined in terms of the cross prod- 
uct. These include the following: 


e A particle of mass m has position vector r and is moving with velocity v. Its 
angular momentum about the origin is h = mr x v. 

e A particle of mass m moves with velocity u in a frame which is rotating with 
angular velocity §2. Due to the rotation, the particle experiences a sideways 
force called the Coriolis force, F = 2mu x N. Since the Earth is rotating, 
this force influences motion on the surface of the Earth. The effect deflects 
particles to the right in the northern hemisphere and is strongest for motions 
on large scales such as ocean currents and weather systems. 

e A particle with electric charge q moves with velocity v in the presence of a 
magnetic field B. This results in a force, called the Lorentz force, equal to 
qv x B. This is the force which is responsible for the operation of an electric 
motor. 
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1.4 Scalar triple product 


The scalar triple product of three vectors a, b and c is defined to be a- (b x c). 
In fact the brackets here are unnecessary: (a - b) x c is meaningless since (a - b) 
is a scalar and so cannot be crossed with the vector c. Therefore the expression 
a b xc is well defined. 

The formula for the scalar triple product in terms of the components of 
the three vectors a, b and c can be obtained using the formula for the cross 
product (1.4): 


a-bxc= a ;b2C3 = a163Cc2 + azb3c, = a2b1c3 + a3bic2 = a3b2c1. (1.8) 


The scalar triple product has a number of properties, listed below. The first 
four follow directly from (1.8). 


e The dot and the cross can be interchanged: 
a-bxc=axbO-c. 
o The vectors a, b and ¢ can be permuted cyclically: 
a-bxc=6b-cxa=c-axb. 
e The scalar triple product can be written in the form of a determinant: 


ai a2 Q3 
a:-bxe= bi ba b3 
C C2 Cg 


e If any two of the vectors are equal, the scalar triple product is zero. 

e Geometrically, the magnitude of the scalar triple product is the volume of 
the three-dimensional object known as a parallelepiped formed by the three 
vectors a, b and c (Figure 1.14). This can be shown as follows. The area of the 
parallelogram forming the base is |bx c|. The height is the vertical component 
of a, which is the magnitude of the component of a in the direction of b x c. 
This is |a- b x c|/|b x c|, so the volume is the area of the base multiplied by 
the height, which is |a-b x c|. Similarly, the volume of the tetrahedron made 
by the vectors a, b and c is |a- b x c|/6. 


The scalar triple product of a, b and c is often written [a, b, c]. This notation 
highlights the fact that the dot and the cross can be interchanged. 
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b 


Fig. 1.14. The volume of the object formed by the three vectors a, b and c is |a-bx c|. 


Example 1.9 
Find the scalar triple product of the vectors (1,2,1), (0,1,1) and (2,1,0). 

First find the vector (0,1,1) x (2,1,0) = (—1,2,—2). Now dot this with 
(1, 2,1), giving the answer 1. 


Example 1.10 

Show that if three vectors lie in a plane, then their scalar triple product is zero. 
If a, b and c lie in a plane, then the vector b x c is perpendicular to the 

plane and hence perpendicular to a. Since the dot product of perpendicular 

vectors is always zero, it follows that a -b x c= 0. 


Example 1.11 
A particle with mass mand electric charge q moves in a uniform magnetic 
field B. Given that the force F on the particle is F = qv x B, where v is the 
velocity of the particle, show that the particle moves at constant speed. 

The equation of motion of the particle is written using Newton’s second 
law, force equals mass times acceleration. The acceleration of the. particle is 
the rate of change of the velocity, written ù, so the equation of motion is 


qux B=mù. 


Now taking the dot product of both sides of this equation with v, the scalar 
triple product on the left-hand side gives zero since two of the vectors are equal. 
Hence 
0=mv-v= le -v)/2= ne (lvl?) /2 
dt dt i 


so the speed of the particle, |v|, does not change with time. 
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1.5 Vector triple product 


The vector triple product of three vectors is a x (b x c). The brackets are 
important here, since a x (b x c) # (a x b) x c. Since only cross products are 
involved, the result is a vector. An alternative expression for a x (b x c) can 
be obtained by writing out the components. Since 


bxc= (b2c3 _ b3c2)e, + (b3c; = bic3)e2 + (bic ~ boc )e3, 
the first component of a x (b x c) is 


[a xX (b x c)|1 = a2(bi c2 = bec; ) = a3(b3c, = bic3) 


= b (a2c2 + a3C3) — C] (a2b2 + a3b3). 
By adding and subtracting the quantity a,b,c,, this can be written 


[a x (b x e)l] = bi (aicy + A2C2 + a3C3) — Cj (a,b, + azb + a3b3) 


ba-c—ca-b. 


Similar equations hold for the second and third components, SO the vector triple 
product can be expanded as 


a x (bx c) = (a - c)b — (a - b)c. (1.9) 
From this result it also follows that 


(a x b) x c = —c x (a x b) = -(c - bja + (c - a)b. (1.10) 


Example 1.12 


Under what conditions are a x (b x c) and (a x b) x c equal? 
By comparing (1.9) with (1.10), the two are equal if —(a- b)c = -(c - b)a. 
This can alternatively be written b x (a x c) = 


Example 1.13 


Find an alternative expression for (a x b) - (c x d). 
Since the dot and cross can be interchanged in a scalar triple product, 


(a x b) - (c x d) a: (b x (c x d)) 
a-((6-d)c — (b-c)d) 


(a -c)(b-d) — (a - d)(b - ©). 
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1.6 Scalar fields and vector fields 


A scalar or vector quantity is said be a fteld if it is a function of position. 
An example of a scalar field is the temperature inside a room; in general the 
temperature has a different value at different points in space, so the temperature 
T is a function of position. This is indicated by writing T(r), where r is the 
position vector of a point in space, r = (x,y,z). Other examples of scalar fields 
include pressure and density. An example of a vector field is the velocity of the 
air within a room. 

In general, a scalar field T is three-dimensional, i.e. it depends on all three 
coordinates, T = T(z,y,z). Such fields are difficult to visualise. However, if 
the scalar field only depends on two coordinates, T = T(2,y), then it can be 
visualised by sketching a contour plot. To do this, the line T(z, y) = constant 
is plotted for different values of the constant. For example, consider the scalar 
field T'(z,y) = x? + y?. The contour lines are the lines z? + y? = constant, 
which are concentric circles centred at the origin, as shown in Figure 1.15(a). 


(a) 


Fig. 1.15. (a) Contours of the scalar field T(z,y) = z? + y”. (b) The vector field 
u(z,y) = (y,2). 


Vector fields in two dimensions can also be visualised by a sketch. In this 
case the simplest procedure is to evaluate the vector field at a sequence of points 
and draw vectors indicating the magnitude and direction of the vector field at 
each point. An example of this procedure is the drawing of wind speeds and 
directions on weather maps. For example, consider the vector field u(x, y) = 
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(y,z). At the point (1,0), u = (0,1), so at this point a vector of magnitude 1 
pointing in the y direction is drawn. Similarly, at (0,1), u = (1,0) and at (1,1), 
u = (1,1). By considering a few additional points, a sketch of the vector field 
can be built up (Figure. 1.15(b)). 


Summary of Chapter 1 


e A vector is a physical quantity with magnitude and direction. 

e A scalar is a physical quantity with magnitude only. 

e In Cartesian coordinates a vector can be written in terms of its components 
as either a = (a), @2,a3) Or @ = Q1€ı + a2€2 + a3€3, where e;, €2 and e3 are 
unit vectors along the z-, y- and z-axes respectively. 


The magnitude of the vector a is |a| = ya? + a2 + aż. 
e The dot product or scalar product of a and 6 is a scalar, 


a- b= |a||b| cos = a,b; + agbz + a3b3. 


This can also be thought of as |a| multiplied by the component of b in the 
direction of a. Applications of the dot product include the work done when 
moving an object acted on by a force and the equation of a plane. 

e The cross product or vector product of a and b is a vector, a x b, with 
magnitude |a||b| sin@, perpendicular to a and b in a right-handed sense. In 
component form, 


a x b= (a2b3 — a3bz)e1 + (a3b, = ab3)e2 + (a; b2 = azbı)ez. 


The magnitude of a x b is |a| multiplied by the component of b perpendicular 
to a, which is the area of the parallelogram made by a and b. Applications 
of the cross product include the equation of a straight line and the rotation 
of a rigid body. 

e The scalar triple productisa-bxc=axb- c=b. cxa. 

e The vector triple product is a x (b x c) = (a c)b — (a - b)c. 

e A scalar or vector quantity is a field if it is a function of position. 
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EXERCISES 


1.8 Find the equation of the straight line which passes through the points 
(1,1,1) and (2,3,5), (a) in parametric form; (b) in cross product 
form. 

1.9 Using vector methods, prove the sine rule, 


sin A = sin B 7 sin C (1.11) 


and the cosine rule, 
c? = a? + b? — 2abcosC (1.12) 


for the triangle with angles A, B, C and sides a, b, c in the figure 
below. 


1.10 (a) Show that the set of vectors and the operation of vector addition 
form a group. (The set of objects a, b,c,... and the operation x form 
a group if the following four conditions are satisfied: (i) for any two 
elements @ and b, ax b is in the set; (ii) (a xb) xc = ax (b x c); (iii) 
there is an element J obeying ax I = I x a = a; (iv) each element a 
has an inverse a~! such that axa~! = a`! xa = I.) 

(b) Do the set of vectors and the dot product form a group? 

(c) Do the set of vectors and the cross product form a group? 
1.11 Simplify the following expressions: 

(a) |a x b|? + (a - b)’; 

(b) a x (b x (a x b)); 

(c) (a — b) - (b — c) x (c—a); 

(d) (a x b)- (b x ce) x (ec xa). 

1.12 The vector x obeys the two equations z-a = 1 and x xa = b, where 
a and b are constant vectors. Solve these equations to find an ex- 
pression for x in terms of a and b. Give a geometrical interpretation 
of this question. 
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1.13 Find the equation of the line on which the two planes r -a = 1 and 
r-b=1 meet. 

1.14 (a) Express the vector a x b in the form aa + 8b + ye, assuming 
that the vectors a, b and c are not coplanar. 

(b) Hence find an expression for (a x b - c)? that does not involve 
any cross products. 

(c) Hence find the volume of a tetrahedron made from four equilat- 
eral triangles with sides of length 1. 

1.15 A particle of mass m at position r and moving with velocity v is 
subject to a force F directed towards the origin, F = —f(r)r. Show 
that the angular momentum vector h = mr x v is constant. 

1.16 Sketch the scalar field T(z, y) = z? — y. 

1.17 Sketch the vector field u(x, y) = (x + y, —2). 
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Line, Surface and Volume Integrals 


2.1 Applications and methods of integration 


This chapter is concerned with extending the concept of integration to vector 
quantities and to three dimensions. Before embarking on these more compli- 
cated types of integration, however, it is useful to review the concept of inte- 
gration and some standard techniques for evaluating integrals. It is important 
that the reader is familiar with these methods, since this will be assumed in 
the following sections. 


2.1.1 Examples of the use of integration 


Example 2.1 


A rod of length a has a mass per unit length p(x) that varies along the length 
of the rod according to the formula p(x) = 1+z. What is the total mass of the 
rod? 

Consider dividing the rod into N small sections, each of length dz;. The 
mass of each section is p(z;)dz;. The total mass M of the rod is the sum of the 
masses of all these sections, 
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N 
M= > p(z;) dzi. 


i=1 


The integral of p(x) is defined to be the limit of this sum as N — oo: 


å N 
|, Pade = Jim Y olz)de:. 


i=1 


The total mass M is therefore 


a 
m= | l +zdr = [z+2?/2], =a+a?/2. 
o 


Example 2.2 


A vehicle starts from rest and accelerates uniformly up to a speed of 10 m/s 
over a time of 20 s. What is the total distance travelled during this time? 

The vehicle starts from rest and reaches a speed of 10 m/s after 20 s, so 
its speed at a time t is v(t) = t/2 m/s. In a small time interval dt the distance 
travelled is v(t) dt = t/2 dt. The total distance S travelled in the total time of 
20 s is therefore 


20 
S =| t/2dt = [#?/4], = 100m. 
0 


2.1.2 Integration by substitution 


In this method for the evaluation of integrals, a complicated integral is trans- 
formed to a simpler one by a substitution or change of variable. In some cases 
the choice of the change of variable is easy to find, but in others it can be 
difficult to spot the most sensible substitution. Often, there is more than one 
possible substitution. Three examples of the application of this method are 
given below. 


Example 2.3 


Evaluate f r/ V1 — T dz. 

Here, the difficulty is caused by the v1 — z in the denominator. This sug- 
gests that the appropriate substitution is u = 1 — z, so z = 1l — u and 
dx = —du. The integral becomes f —(1 — u)/Vudu = f —1/vu + vudu = 
—2u!/? + 2u3/2/3 +c, where c is an arbitrary constant of integration. The result 
can be expressed in terms of the original variable z as -2V1 — zT (2 + 1)/3 + c. 
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Example 2.4 
Evaluate f V1 — 2? dz. 


For integrals involving the quantity /1 — x?, the appropriate substitution is 
z = sin (or x = cos@, which would do equally well). With this choice, V1 — x? 
becomes cos@ and dr = cos dŷ. The integral then simplifies to f cos? 6 dd. 
Integrals of this type, which occur very frequently, are evaluated using the 
trigonometric formula cos? @ = (1 + cos20)/2, so the value of the integral is 
(20 + sin 20)/4 + c. In terms of z this result can be written 


[Vi-2 dz = (sints + zV/1— 2?) /2+ c. 


Example 2.5 
Evaluate fj z? VI- 2? dz. 


Again the substitution z = sin ĝ is used. Since the limits of integration are 
given, these can also be expressed in terms of the new variable 8. When z = 0, 
0 = 0 and when z = 1, 0 = 7/2, so the integral becomes ie * sin? 8 cos? 6 dé. 


This can be simplified to ie 1/4 sin? 20 dé = Hi 1/8 (1 — cos 40) dð = 7/16. 


2.1.3 Integration by parts 


Integration by parts is an important and useful technique, used when an integral 
involves a product of two terms. The integration by parts formula is derived 
from the product rule for differentiation. Given two functions of z, u(x) and 
v(x), the rule for the derivative of their product is 

d(uv) dv 


du 
de u F Yir (2.1) 


Integrating this expression and rearranging the terms gives the integration by 


parts formula: 7 F 
v u 
fugar=uv— |v a. (2.2) 


As with the case of integration by substitution, some experience is helpful in 
determining whether this formula will be useful in evaluating an integral, and 
exactly how to split the integral into the two parts. In general, it is best to 
choose u to be a function which becomes simpler when differentiated. The 
following two examples illustrate the use of the method of integration by parts. 
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Example 2.6 


Evaluate f xsin z dz. 
In this example we choose u = z, du/dz = sinz, so v = — cos z. Applying 
the formula (2.2) gives 


[ssinzde = -zcosz+ | coszdz = —z COs T + sinz +c. 


Note that it is essential to make the right choice for u and v. If we had chosen 
u = sinz, dv/dzr = z then the resulting integral would have involved z? cos z 
which is more complicated than the integral we started with. 


Example 2.7 


Evaluate f exp az cos z dz. 
For this case two applications of (2.2) are necessary. Choosing u = exp az, 
du/dzx = cosz, 


f epas cosz dx = expazsinz — J cexpar sins dz +c. 
Now apply (2.2) to the integral on the r.h.s. with u = aexpaz, dv/dz = sin z: 
[ expazcoszds = exp az sin T + aexpar cost — fe expar cos z dz +c. 


At this stage it may appear that no progress has been made, since the original 
integral has reappeared on the right-hand side. However, by rearranging the 
terms, 


(1 + a”) f epas cos z dz = expazsinz + aexpaz coss +c 
and so the value of the original integral is 
[exper cos z dz = (expazsinz + aexpaz cosa + c)/(1+ a’). 
In this case the choice of u and v does not matter: the result can also be obtained 
by choosing u = cos z, du/dz = expaz, provided that a similar choice, with u 


chosen to be the trigonometric term and du/dz chosen to be the exponential 
term, is made for the second application of (2.2). 
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2.2 Line integrals 


2.2.1 Introductory example: work done against a force 


As an introductory example of a line integral, consider a particle moving along 
a curved path C through space. The particle is acted on by a force F(r), which 
is a vector field. What is the total amount of work done as the particle moves 
along the curve C? 

To answer this question, first divide the curve C into a large number of small 
pieces. Consider the work done when the particle moves from the position r 
to r + dr (Figure 2.1). On this small section of the curve C, the work done is 
—F - dr. The total amount of work done W as the particle moves along C is 
therefore the sum of the contributions from all the small segments of the curve, 


N 
W=) > -Fi- dri. (2.3) 


t=1 


The line integral of F along the curve C is defined by 


N 
f. F -dr = fm, F, dri. (2.4) 


The vector dr is often referred to as a line element. 

Note that the direction of the integral along the curve C must be specified. 
If the direction of the curve is reversed, all the line elements dr are reversed 
and so the value of the integral is multiplied by —1. 


dr 


r+ar 


C O 


Fig. 2.1. A small section of the curve C is represented by the line element dr. 
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2.2.2 Evaluation of line integrals 


Line integrals are evaluated by using a parameter, t for example, together with 
a formula giving the value of the position vector r in terms of t. This can be 
regarded as an example of integration by substitution, since 


[Fara |r- Ea (2.5) 
A di 


For example, suppose that the curve C is given in terms of t by 
tSt, y =t, z = 2t, (2.6) 


and t lies in the range 0 < t < 1. Then as ¢ varies between 0 and 1, the position 
vector r = (x,y,z) moves along a curve C in space connecting the points 
(0,0,0) and (1,1,2). Suppose that the vector field F is given by F = (y, 2, z). 
To evaluate the line integral, both F and dr/dt must be written in terms of t. 
Substituting (2.6) into the expression for F gives F = (t,t, 2t?), and 


dr (dx dy dz\ _ 
dt = (Gd. Z) = (1,1, 4¢). 


The line integral can now be evaluated: 


1 1 
| F-ar= | (62) (11,48) dt = f ot + 8 dt = 3. 
C. 0 0 


In this first example, the parametric form of the curve C was given. If the 
curve C is given in a different form, then a parametric form must be constructed 
so that the line integral can be evaluated. For example, suppose now that 
F = (y,z,z) as before, but C is the straight line connecting the origin to the 
point (1, 2,3). The way in which the curve C is parametrised is not unique, so 
we can make the arbitrary choice z = t. Since x varies between 0 and 1 along 
the line, this is also the range for t. The end point of C is (1,2,3), so y and 
z must be given by y = 2t, z = 3t, and so dr = (1,2,3) dt. The value of the 
integral is therefore 


1 1 
[Fare f (2t,t,3t) -(1,2,3)at = | 13t dt = 6.5. 
C 0 0 


Line integrals sometimes occur over curves that are closed, i.e. when the 
starting point and end point of the curve are equal. In this case the integral is 
written using the symbol ¢, which indicates that the integral is along a closed 
curve. For example, consider the integral of F = (y,z,z) around the closed 
curve given by z = cos@, y = sin@, z = 0, where 0 < @ < 27. Here, as 0 varies, 
the curve C describes a circle in the z, y plane. The line element dr is expressed 
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in terms of the parameter 0 as dr = (dz, dy,dz) = (— sin, cos 8,0) dO, so the 
value of the line integral is 


2 
$ Far = f (sin 8, cos 8, 0) - (— sin 8, cos 8, 0) d8 
C 0 


2n 
f — sin? 6 + cos? 0 d8 
0 


2n ; 
i cos 20 d0 = [1/2 sin 20)?" = 0. 
0 


The line integral of a vector field F around a closed curve C is often called the 
circulation of F around C. 
Example 2.8 


Evaluate the line integral of the vector field u = (ry, z*, x) along the curve 
given by z=1+t,y=0,z=?7,0<t<3. 

First write u and dr in terms of t: u = (0,#4,1 + t), dr = (1,0, 2t) dt. The 
value of the integral is therefore 


3 3 
[ F-ar= | (0,t8,1 +4): (1,0,24) de = f 2t + 2t? dt = t? + 2t3/3]° = 27. 
C 0 0 
Example 2.9 


Find the line integral of F = (y, —z,0) along the curve consisting of the two 
straight line segments (a) y= 1,0<2< 1, (b) z=1,1<y<2. 

Here, the contributions from the two line segments must be taken separately. 
On section (a), using z as the parameter, we have 


[av-2.0) 7 (dz, 0, 0) =L 
0 


Similarly on section (b) we have 


2 
f (y, -1, 0) i (0, dy, 0) ==], 
1 
Therefore the total value of the integral is 0. 


Example 2.10 


Find the circulation of the vector F = (y, —z,0) around the unit circle, z? + 
y? = 1, z = 0, taken in an anticlockwise direction. 

The circle is-written in terms of a parameter ĝ as z = cos@, y = sin@, z = 0, 
0 <8 < 27. The value of the integral is 


2n 27 
$ F dr = (sin 6, — cos 8,0) - (— sin 8, cos 8, 0) dô = f —1d@ = —2r. 
C 0 0 
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2.2.3 Conservative vector fields 


A vector field F is said to be conservative if it has the property that the line 
integral of F around any closed curve C is zero: 


$ F-dr=0. (2.7) 
C 


An equivalent definition is that F is conservative if the line integral of F 
along a curve only depends on the endpoints of the curve, not on the path 
taken by the curve, 

F -dr = F dr (2.8) 
C'i C2 
where C and C» are any two curves that have the same endpoints but different 
paths (Figure 2.2). 


A 


Fig. 2.2. A vector field is conservative if the line integrals along two different curves 
Cı and C: from A to B are equal. 


The equivalence of these two definitions can be demonstrated as follows. 
‘Consider two curves C and Cp that start from the point A and end at the 
point B (Figure 2.2). Let C be the closed curve that starts from the point A, 
follows the curve C; to the point B and then follows the curve C% in the reverse 
direction to return to A. Then 


$ F-dr = F-dr— | F-dr (2.9) 
C C1 C2 

since the effect of the reversed direction of the integral along C% is to change 
the sign of the integral. From this equation it follows that if the integral around 
C is zero, then the integrals along Cı and C2 are equal, and similarly, if the 
integrals along Cı and C} are equal, then the integral around C is zero. 
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Conservative vector fields are of great importance, since many physical ex- 
amples of vector fields are conservative. Consider for example the Earth’s grav- 
itational field, g. A particle of mass m experiences a force mg, so the work 
done in moving the particle along a path C from point A to point B is just 
minus the line integral of mg along C'. However, we know physically that the 
work done only depends on the position of the points A and B - in fact the 
work done is simply the difference in the potential energy of the particle at A 
and B. Equivalently, if the particle is moved around but returns to its starting 
point, the total work done is zero. Therefore, the Earth’s gravitational field is 
an example of a conservative vector field. 


Example 2.11 
By considering the line integral of F = (y, xz? — z,0) around the square in the 
z, y plane connecting the four points (0,0), (1,0), (1,1) and (0,1), show that 
F cannot be a conservative vector field. 

This line integral consists of four parts (Figure 2.3). On the first section 


C; 


Ci x 


Fig. 2.3. The line integral around the square is split into four straight sections 


Cı, from (0,0) to (1,0), dr = (dz,0,0) and y = 0 so F = (0,2? — z,0) and 
F -dr = 0. On the second section C2, dr = (0, dy,0) and z = 1 so F = (y,0,0) 
and again F -dr = 0. On the third section C3, dr = (dz,0,0) and y = 1 so 
F =(1,2z? — 2,0) and the contribution to the line integral is 


0 
J ldz = -1. 
I 


Finally, on the fourth section C4 of the square F - dr = 0 again, so the total 
value of the integral is —1. Since the integral around the closed circuit is non- 
zero, F cannot be conservative. 
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2.2.4 Other forms of line integrals 


The line integral fo F -dr is the most important type of line integral, but there 
are two other forms of the line integral which can occur. These are 


f ar and [ Pxar 
C C 


where @ is a scalar field and F is a vector field. Note that in each of these cases 
the result of the integral is a vector quantity. These integrals can be evaluated 
using a parameter, as in the following examples. 

Example 2.12 

Evaluate the line integral 


| z+y’dr, 

C 

where C is the parabola y = z? in the plane z = 0 connecting the points (0,0,0) 
and (1,1,0). 


The curve can be written in terms of a parameter t as z = t, y = t?, z = 0, 
0 <t <1, so dr = (1, 2t, 0) dt. The value of the integral is therefore 


[a+var 
C 


f t + t*)(1, 2t, 0) dt 
0 


1 1 
e] (/ t+ tat) + (f 2? +24 dt) = 0.7e, + eo. 
0 0 
Example 2.13 


Evaluate the line integral 


| Exar, 
C 


where F is the vector field (y, z,0) and C is the curve y = sin z, z = 0, between 
z=Oandz=n7. 

The curve can be written as z = t, y = sint, z = 0,0 < t < am. Then 
F = (sin t, t,0) and dr = (1, cost, 0) dt, so F x dr = (0,0, sin t cost — t) dt and 
the integral is 


N 
| Fxar=es | sint cost — t dt = 1/2 [sin? t — t°]; e3 = -7° /2 e3. 
C 0 
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EXERCISES 
2.1 Evaluate the line integral 
J Pdr where F=(5z?,22,2 + 2y) (2.10) 
C 


and the curve C is given by z = t, y = t?, 2 =t,0<t<1. 

2.2 Evaluate the line integral of the same vector field F' given in (2.10) 
along the straight line joining the points (0,0,0) and (1,1,1). Is F 
a conservative vector field? 

2.3 Find the line integral of the vector field u = (y?, z, z) along the curve 
given by z = y =e” from z = 0 toz = 1. 

2.4 Find the line integral far x dr where the curve C is the ellipse 
T? /a? + y?/b? = 1 taken in an anticlockwise direction. What do you 
notice about the magnitude of the answer? 


E 


2.3 Surface integrals 


2.3.1 Introductory example: flow through a pipe 


Suppose that fluid flows with velocity u through a pipe. What is the total 
volume of fluid passing through the pipe per unit time (Figure 2.4)? This 


ee ——— 
u 
—— —— 


Fig. 2.4. Fluid flows with velocity u along a pipe. The rate at which it crosses the 
surface S at the end of the pipe is an example of a surface integral. 


volume flow rate is often called the fluz of fluid through the pipe, or the flux 
of fluid across the surface S that forms the end of the pipe. We will consider 
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this question in several different cases, beginning with the simplest case and 
progressing to more complicated examples. 

Suppose first that the pipe has cross-sectional area A and that the velocity 
u, which in general can be a function of space and time, is a constant and is 
directed parallel to the walls of the pipe, with speed |u| = Uo . In this case, 
the fluid moves along the pipe as if it were a solid block. In a time t, the fluid 
moves a distance Uot, so a ‘block’ of fluid of volume UptA emerges from the 
end of the pipe. The flow rate Q, or flux, of fluid through the pipe is therefore 
this volume divided by the time t, giving Q = UoA. 


Fig. 2.5. Enlargement of a small surface element dS forming part of the surface S. 


Now suppose that the flow is again directed parallel to the walls of the 
pipe but that the speed of the flow depends on the position within the pipe, so 
lu| = Up(z,y), and that the pipe has a square cross-section with the walls at 
x = 0,1 and y = 0,1. Now consider a small surface element with area dS on 
the surface S, which is a small rectangle with sides of length dz and dy located 
at the point (x,y) on the surface S, so that dS = dz dy (Figure 2.5). Following 
the argument of the previous paragraph, the flux dQ of fluid across this surface 
element dS is dQ = Uo(x,y) dS = Uo(z, y) dz dy. To calculate the total flux Q 
across the surface S, we need to add up the contributions from all the small 
surface elements dS. This sum of contributions becomes an integral, but since 
the surface is two-dimensional the resulting integral is a surface integral or 
double integral, representing integration in both the z and y directions: 


Q= If Uo(z, y) dS = If. Uo(, y) dz dy. (2.11) 


In the above example, the fluid flow direction is perpendicular to the surface 
S. Consider now the case where the vector field u and the surface S are both 
arbitrary. In general, S may represent a curved surface. Again we consider a 
small surface element dS (Figure 2.6) and compute the flux of u across dS. Now 
if u is perpendicular to dS, this flux is just |u|dS, but if u is not perpendicular 
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Fig. 2.6. The unit normal vector n is perpendicular to the surface element dS. 


to dS, only the component of u perpendicular to dS contributes to the flux 
across dS. To extract this component it is necessary to introduce a normal 
vector n to the surface dS, with the properties that n is perpendicular to dS 
and that n isa unit vector, |n| = 1. The component of u perpendicular to 
dS is then the component of u in the direction of n, which is just u-n (see 
Section 1.2), so the flux across the surface element dS is u- ndS. The total 
flux across the surface S is given by the surface integral 


Q= || u-nds= ff u-ndzdy. 


Note that the direction of the normal vector n has not been uniquely spec- 
ified: the vector n could have been chosen to point in the opposite direction 
in Figure 2.6, and this would change the sign of the answer. Therefore the 
direction of n must be specified when a surface integral is written down. 

Surface integrals often occur over surfaces which are closed. In this case, the 
normal to the surface which points outward is used (Figure 2.7). To indicate 
that the surface is closed, a circle is sometimes drawn through the integral, in 
a similar manner to the notation used for line integrals around closed curves: 


ff unas. 
S 


2.3.2 Evaluation of surface integrals 


Surface integrals can be evaluated by carrying out two successive integrations. 
Consider first the evaluation of (2.11) where S is the square surface given by 
0<2<1,0<y <1. The surface integral is then 


o= f| venas = f i Uo(z,y)dz dy. 
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Fig. 2.7. For a closed surface, the convention is that the normal points outward. 


There is some potential ambiguity here, since it is not immediately clear which 
of the integral signs refers to the z integration and which to the y integration. 
The convention adopted is that the integrals are ‘nested’, so that the first 
integral sign represents the y integral and the second one represents the =z 
integral. The double integral is then interpreted as 


Q= [ F Volz, y) dz) dy. 


It is also assumed that the inner integral, the z integral in the above equation, 
is to be evaluated first. For example, suppose that Uo(z, y) = (z — 2”)(y — y’”). 
The inner, z, integral is evaluated first, and within this inner integral y is 
regarded as a constant. Carrying out the first integral gives 


J | (z — 2?)(y — y?) dz dy 
J [z?/2 — 23/3] (y - y?) dy 


Q 


J 1/6 (y — y*) dy. 


The double integral has now been reduced to a single integral which can be 
evaluated in the usual way: 


Q= [1/6 (y—v?)ay = 1/6 [p?/2 - 4°13] = 1/36. 


Consider now the case where the surface is circular. In the original example 
of pipe flow, this corresponds to a cylindrical pipe. If the surface is circular, it 
is best to use polar coordinates to evaluate the integral. Polar coordinates (r, 0) 
are related to Cartesian coordinates (z, y) by x = r cos 0, y = r sin 8. The angle 
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ð is measured in radians and ranges from 0 to 27. The area element dS can be 
computed by considering a small angle dô and a small change in radius dr as 
shown in Figure 2.8. If dr and d@ are both small then the corresponding area 


| rdð 


dr 


Fig. 2.8. The area element dS in polar coordinates. 


element dS is almost rectangular, with length dr in the r direction and r dé in 
the @ direction, so dS = r dô dr. Suppose now that the radius of the surface S 
is 1 and that Up = 1 — r°. The value of the surface integral is then 


Q= f was = f [a-rao 


The inner, 0, integral is carried out first, with r temporarily regarded as a 
constant. Since there is no dependence on @ in the integral, this inner integral 
just gives a factor of 27, so 


1 
Q= i 2r(1 — r°)rdr = 2r [r?/2— r*/4], =m /2: 


Finally, consider the case where the surface S is curved. The surface can 
be written in terms of two parameters, v and w, so that a position vector r 
lying in the surface is written r = r(v, w). Now consider a small change in the 
value of v, to v + dv. The vector r(v + dv, w) also lies on the surface, so the 
difference between these two vectors, r(v + dv, w) —r(v, w) = (Or /Ov) dv must 
be a vector lying in, or tangent to, the surface, and similarly for the vector 
(Or /Ow) dw. To evaluate the surface integral we need an expressions for n dS, 
but we know from Section 1.3 that this is simply the cross product of the two 
vectors, since the cross product of two vectors gives a vector perpendicular to 
both and with a magnitude equal to the area of the parallelogram created by 
the two vectors. Therefore the surface integral can be written as 


Or Or 
[funds = || uF Ewa. 


For example, consider the integral of the vector field u = (z,z, —y) over the 
curved surface of the cylinder x? + y? = 1 lying between z = 0 and z = 1. The 
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Fig. 2.9. A point on the cylindrical surface z? + y? = 1 can be denoted by the 
parameters z and @. 


two parameters describing the surface are the height z and the angle @ around 
the cylinder (Figure 2.9). In terms of these parameters the position vector is 
r = (x,y,z) = (cos@, sin 8, z), so 


Or ; Or _ 
Or Or F 
and 30 x Əz = (cos®ð, sin 0,0). 


The value of the integral is therefore ' 


[pers 


1 p2r f 
fl | (cos 8, z, — sin 8) - (cos 6, sin 8, 0) dé dz 
0 Jo 


1 2n 
f / cos? 0 + z sin 6 d0 dz 
o Jo 


1 
= f ndz=M. 
0 
Example 2.14 


Evaluate the surface integral of u = (y,z”,z”), over the surface S, where S is 
the triangular surface on z = 0 with y > 0, z > 0, y +z < 1, with the normal 
n directed in the positive z direction. 

In this example n = (1,0,0) and so u-n.= y. When the surface is not 
rectangular, care must be taken when setting the limits of the integration. If 
we choose to do the z integral first, then for any given value of y, the range 
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0 
0 1 9 


Fig. 2.10. For each value of y, z runs from the line z = 0 up to the line y + z = 1. 


of values for z is 0 < z < 1 — y (Figure 2.10). The outer y integral then has 
limits of 0 and 1. This corresponds to covering the triangular area with vertical 
strips. The value of the integral is 


1 pl-y 1 isj 1 
f f ydzdy= | [yz], ay = | y — y? dy = 1/6. 
o Jo 0 0 


The integral could also have been evaluated by doing the y integral first, in 
which case the limits for y are 0 < y < 1 — z and the limits for z are 0 < z < 1. 
This ordering corresponds to covering the region of integration in Figure 2.10 
with horizontal strips. The value of the integral is the same: 


i [va dz = [ [y2/2)3-* Tz [o — 2)?/2dz = 1/6. 


Example 2.15 


Find the surface integral of u = r over the part of the paraboloid z = 1—2? —y? 
with z > 0, with the normal pointing upwards. 

Since the surface is curved, a description of the surface in terms of two 
parameters is needed. Using simply z and y, a point on the surface is (z,y,1-— 
z? — y?) and the two tangent vectors in the surface, obtained by differentiating 
with respect to z and y, are (1,0, —2x) and (0, 1, —2y). Taking the cross product 
of these two vectors, the quantity n dS is (2x, 2y, 1) dz dy. Note that this has 
a positive z component, so is directed upwards as required. Taking the dot 
product with u gives u-ndS = 2r? + 2y? + zdz dy = 1 + z? + y? dz dy. The 
limits on the integral are determined in a similar manner to Example 2.14. The 
edge of the surface is given as z = 0, which is the circle z? + y? = 1. Choosing 
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Fig. 2.11. For any value of z, y ranges from the lower half of the circle z? + y? = 1 
to the upper half of the circle. 


to do the y integral first, the y integration is carried out with z fixed, so the 
range of values for y is -v1 — z? < y < V1 — r? (Figure 2.11). The range of 
values for z is —1 < x < 1. The value of the integral is 


1 V1—z? l , 
J J 1 + x° + y? dy dr 
~1 J—V1—22 


iz ly +2y +t] ae 


[ (8/3+4/327)V1 — r? dz 


4n/3+ 7/6 = 37/2, 


where the final integral has been evaluated using the results of Section 2.1.2. 


2.3.3 Other forms of surface integrals 


The surface integral of u : n is the most important type of surface integral. 
However, as in the case of line integrals, other types of surface integrals can be 
defined, for example 


Jta. Jf tnas and J| xas, 


where f is a scalar field and v is a vector fieid. These integrals are evaluated 
using the methods of the previous section. l 
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Example 2.16 
If S is the entire z, y plane, evaluate the integral 


r= || e77"? dS, 
S 


by transforming the integral into polar coordinates. 
In polar coordinates (r,0), z? + y? = r? and dS = r dô dr. The ranges of 
the variables to cover the whole plane are 0 < r < œ and 0 < 0 < 2r, so 


CO pf2r 5 CO r 2700 
I aj f e`” rdôdr =f 2re" rdr=r [-e7" | = 7. 
0 o 0 0 


This answer can be used to show that [9 e-** dr = V/7, a result which cannot 
be obtained by standard methods of integration. 


2.4 Volume integrals 


2.4.1 Introductory example: mass of an object with 
variable density 


Suppose that an object of volume V has a density p. If p is a constant, the 
mass M of the object is simply M = pV. Now suppose that the object has a 
density which is a function of position, p = p(r). What is the total mass of the 
object? 

The argument proceeds in a similar manner to the examples of line and 
surface integrals. The volume V is divided into N small pieces with volumes 
ôV;, i =1,...,N, which are called volume elements. Within each of the volume 
elements the density is approximately constant (assuming that p is a continuous 
function of position) and so the mass M; of the volume element at position r; 
is M; = p(r;) 6V;. The total mass of the object is the sum of the masses of all 
the volume elements, 


N 
M = 5 p(ri) ÔV;. 
i=l 


The volume integral of p over the volume V is defined to be the limit of this 


sum as N > oo: N 
J| ea = aim, Do alr) aV: (2.12) 
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Volume integrals can also be used to compute the volumes of objects, in which 
case p = 1 in the above example. Note that volume integrals usually occur as 
integrals of scalar quantities. However, the volume integral of a vector field u, 


fff uw 


can be defined in a similar way. 


2.4.2 Evaluation of volume integrals 


Volume integrals are evaluated by carrying out three successive integrals. The 
same rule for the evaluation of the triple integral applies as for double integrals: 
the inner integral is evaluated first. The main difficulty in this process is in 
determining the correct limits for the integrals when the shape of the object is 
complicated. It is often helpful to sketch the region of integration in order to 
find the limits on the integrals. Also useful is the rule that in general, the limits 
on an integral can depend only on the variables of integrals that lie outside that 
integral. For example, if the integrals are evaluated in the order x, y, z then 
the limits on the y integral may depend on z but not on z. 


Example 2.17 


A cube 0 < z,y, 2, < 1 has a variable density given by p= 1+2z+y+2z. What 
is the total mass of the cube? 


The total mass is 
w= fo 
V 
l pl pl 
[ | | 1+z+y+zdzdyaz 
o Jo Jo 


1 pl 
f f [£ + 1?/2 + zy +22], dy dz 
o Jo 


[ [ er+u+aa 


1 
J [3y/2 + y?/2 + yz], dz 
0 


fare dz 
[2z + 27/2], = 5/2. 


Note that in this example the integrals were carried out in the order z, y, z, but 
any other choice of ordering is equally valid. 
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Example 2.18 


Find the volume of the tetrahedron with vertices at (0,0,0), (a,0,0), (0, b, 0) 
and (0,0, c). 


Fig. 2.12. For any given values of z and y, z ranges from the plane z = 0 up to the 
plane z/a + y/b + z/c = 1. This is indicated by the vertical column. 


A sketch of the tetrahedron is shown in Figure 2.12. The faces of the tetra- 
hedron are the planes x = 0, y = 0, z = 0 and the plane which passes through 
the three points (a,0,0), (0,5,0) and (0,0,c). The equation of this plane is 
z/a+y/b+z/c = 1, which can be deduced from the general formula for the 
equation of a plane (1.3). Suppose that we choose to do the z integral first. 
This integral is carried out for fixed values of z and y, so the range of z is 
from the plane z = 0 to the plane z = c(1 — z/a — y/b). Choosing to do the 
y integral next, y ranges from 0 to the line that passes through (a,0,0) and 
(0,b,0), which is y = 0(1 — z/a). Finally the range of z is from 0 to a. The 


volume V is therefore 
a pb(l—z/a) pe(l—z/a—y/b) 
f J f dz dy dz 
o Jo 0 


a pb(1—z/a) 
f J c(1 — z/a — y/b) dy dz 
o Jo 


ef [y(1 — z/a) - y? /20] 07 dx 
0 


cb f° 2, _ abc 
Sfa z/a) dr = —. 


V 


Note that this result is consistent with the formula for the volume of a tetra- 
hedron given in terms of the scalar triple product in Section 1.4. 
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Summary of Chapter 2 


e There are several different types of integral, but each one should be inter- 
preted as the limit of a sum. 
| F -dr, 
C 


e A line integral, written 
represents the sum of the elements F - dr along the curve C. Applications 
include the total amount of work done when a particle moves in the presence 
of a force that is a function of position. 

e Line integrals are evaluated by writing the vector F and the curve C in terms 
of a parameter, t. 

e If the line integral of F around any closed curve is zero, F is said to be 


conservative. 
J | u:nds, 
S 


e The surface integral, 
represents the flux of u across the surface S; this can be thought of as the 
volume of fluid flowing with velocity u across the surface S per unit time. 
The normal vector n is a unit vector that is perpendicular to the surface S. 
e Surface integrals are evaluated by carrying out two successive integrations. 


represents the sum of pdV over all the volume elements dV contained within 
V. If p is the density, the volume integral gives the total mass of the object 
with volume V. 

e Volume integrals are evaluated by carrying out three successive integrations. 
Care must be taken over setting the limits of the integrals and over the order 
in which the three integrals are evaluated. 

e In both surface and volume integrals, the inner integral is evaluated first. 

e There are other forms of line, surface and volume integrals, but the forms 
displayed above are the most commonly occurring. 
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EXERCISES 


2.5 Evaluate the surface integral of u = (ry,2,r + y) over the surface 
S defined by z = 0 withO < xz < 1,0 < y < 2, with the normal n 
directed in the positive z direction. 

2.6 Find the surface integral of u = r over the surface of the unit cube 
0 < z,y, 2 <1, with n pointing outward. 

2.7 The surface S is defined to be that part of the plane z = O lying 
between the curves y = x? and z = y”. Find the surface integral of 
u:n over S where u = (z,zy, x?) and n = (0,0,1). 

2.8 Find the surface integral of u -n over S where S is the part of the 
surface z = z + y? with z < 0 and z > —1, u is the vector field 
u = (2y + z, —1,0) and n has a negative z component. 

2.9 Find the volume integral of the scalar field ¢ = z? + y? + z? over the 
region V specified byO<2<1,1l<y<2,0<2z<3. 

2.10 Find the volume of the section of the cylinder z? + y? = 1 that lies 
between the planes z = z + 1 and z = -z — 1. 

2.11 A circular pond with radius 1 m and a maximum depth of 1 m has 
the shape of a paraboloid, so that its depth z is z = 1 — z? — y?. 
What is the total volume of the pond? How does this compare with 
the case where the pond has the same radius and depth but has the 
shape of a hemisphere? 


I 


Gradient, Divergence and Curl 


3.1 Partial differentiation and Taylor series 


This chapter introduces important concepts concerning the differentiation of 
scalar and vector quantities in three dimensions. These concepts form the core 
of the subject of vector calculus. In this preliminary section, the methods of 
partial differentiation and Taylor series are reviewed. 


3.1.1 Partial differentiation 


Consider a scalar quantity f which is a function of three variables, so f = 
f(x,y, z). Then the partial derivative of f with respect to z is defined to be the 
derivative of f with respect to z, regarding y and z as constants. To indicate 
that f is a function of more than one variable, the partial derivative is written 
using a curly d, 0. More formally, the definition of the partial derivative is 


Of _ . f(x+dz,y,z) — f(x,y,z) 
Or Pre ox (21) 


Second derivatives, such as 0? f /ðz? and mixed derivatives, such as 8? f /OyOz 
can also be defined: the mixed derivative means that f is differentiated with 
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respect to x regarding y as a constant, and then differentiated with respect to 
y regarding z as a constant. An important property of this mixed derivative, 
or cross-derivative, is that the order of the two derivatives does not matter, i.e. 


Of _ f 
Ardy Oydr’ 
provided that these second partial derivatives exist and are continuous. 
The main applications of partial differentiation are in the following sections. 
One additional application is in finding the maximum or minimum of a function 
of more than one variable: when a function is at a maximum or a minimum, 
all of its partial derivatives are zero. 


(3.2) 


Example 3.1 
The function f(z,y,z) is defined by f = z? + zysinz — yz. Find the partial 
derivatives of f with respect to z, y and z and verify the result (3.2) that the 
order of partial differentiation does not matter. 

Differentiating f with respect to z, y and z in turn gives the three partial 
derivatives 

ð : ð 

LEE E A ee E oF = sy cosz -y 


Ox Oy 
Differentiating ðf / dy with respect to x gives 
6? f ð (ðf ; 
ðrðy = Or (3) = Sin z, 
and similarly, differentiating Of /ðx with respect to y gives 
of ð (ðf : 
OyOz = ay (52) = sn z. 


In the same way it can be confirmed that the ordering of the cross-derivatives 
in z and z or in y and z does not matter. 


Example 3.2 


A rectangular box has height a, length b and breadth c, and is open at the top. 
If the volume of the box is fixed, deduce how a, b and c should be related to 
minimise the surface area of the box. 

Let the fixed volume of the box be V, so abe = V. Two sides of the box have 
area ab, two have area ac and the base has area bc, so the total surface area is 
A = 2ab + 2ac + bc. Using the constraint abc = V, a can be eliminated so that 
A = 2V/c+2V/b+ bc. At a maximum or a minimum, the partial derivatives of 
A with respect to b and c must both be zero. This gives the two simultaneous 
equations 
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—2V/b? +c=0, —2V/c’ +b=0. 


From the first equation, c = 2V/b? and substituting this into the second equa- 
tion gives the solution for b, b? = 2V or b = (2V)?/3, and similarly c = (2V)}/3. 
The condition abc = V gives a = V/((2V)?/3) = (2V)!/3/2. Therefore in the 
arrangement that minimises the surface area, b = c and a = 6/2, so the height 
of the box is half its length. 

Note that we have not shown that this is a minimum and not a maximum. 
However, common sense suggests that since the area would be very large if the 
box were tall and thin, the solution found probably does represent a minimum. 


3.1.2 Taylor series in more than one variable 


The Taylor series for an infinitely differentiable function f(z) of a single variable 
1S 


fle) = f(a) +(2-a)Leqy4 EE a.. 
See =o" cs a). (3.3) 
This can also be AA 
ôf = ôr £ pa 2i oe (3.4) 


where ôr = (xz — a) is a small perturbation and ôf = f(z) — f(a) is the 
corresponding perturbation in the value of the function. 

Taylor series can also be constructed for functions of more than one variable. 
For a function f(z,y) of two independent variables, the analogous formula is 
of (62)? f | (dy)? f 6° f 
ôf = ba + by Sh + z! “al aat a 2 Sa + röva +... (3.5) 

In the following sections we will make use of the Taylor series for a function 
f(x,y, z) of three variables, but in all cases only the linear terms, that is, only 
those that only involve a single power of dz, dy or 6z, will be needed: 


fabs SE tovt + az 2ra r (3.6) 


Taylor series can be useful for approximating functions, as in the following 
example. 
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Example 3.3 


Find an approximate value for the function f(z, y,z) = 2x + (1 +y)sinz at the 
point z = 0.1, y = 0.2, z = 0.3. 

The function f takes the value 0 at the point (0,0,0). Near to this point, 
the function can be approximated by its Taylor series expansion. To do this, 
the three partial derivatives of f evaluated at (0,0,0) are required. These are 

of_, of _. Of _ 


aa : oe 3, 7 (l1 +y)cosz=1. 


Hence the Taylor expansion (3.6) is 
Of = 2r +z +... 


which at the point (0.1,0.2,0.3) gives the approximate value f œ% 0.5 (the 
correct value to four decimal places is 0.5546). 


3.2 Gradient of a scalar field 


In Section 1.6 the concept of a scalar field was introduced as a scalar quantity 
which is a function of position in space. A scalar field f can be visualised in 
terms of the level surfaces or isosurfaces on which f is constant. The gradient 
of the scalar field f is a vector field, with a direction that is perpendicular to 
the level surfaces, pointing in the direction of increasing f, with a magnitude 
equal to the rate of change of f in this direction (Figure 3.1). 

The gradient of a scalar field f can be written as grad f, but the gradient is 
so important that a special symbol for grad, V, is used, so grad f = V f. This 
symbol is sometimes referred to as ‘del’ or ‘nabla’. 

The gradient of f can also be defined in a Cartesian coordinate system in 
terms of the partial derivatives of f: 

Vf= see, + sees + Les, (3.7) 
We will now show that these two definitions are equivalent, by showing that the 
vector Vf defined in (3.7) satisfies the two conditions of being perpendicular 
to the level surfaces and with magnitude equal to the rate of change of f in 
this direction. 
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f = constant Vf 


Fig. 3.1. The gradient of f is a vector perpendicular to the surface f = constant. 


Consider an infinitesimal change in position in space from r to r + dr. This 
results in a small change in the value of the scalar field f from f to f + df, 
where, from (3.6), 


df = SI ar + Lay + cL dz 
_ (of Of df 
= & oe ay 3z) (dx, dy, dz) 


Now suppose that dr lies in the surface f = constant (Figure 3.2). In this 
case the change in the value of f must be zero, so we have df = Vf -dr = 0. 
Now in general Vf # 0 and dr # 0, so the two vectors Vf and dr must be 
perpendicular. Since dr is in the level surface f = constant, the vector V f 
must be perpendicular to the level surface. 

So the vector V f defined by (3.7) has the correct direction; it remains to be 
shown that it has the correct magnitude. This is achieved by using (3.8) with 
dr = nds, where n is the unit normal to the level surface and s is a distance 
measured along the normal. In this case, df = Vf -nds = |Vf|ds, since Vf 
and n are parallel and |n| = 1. Hence the magnitude of V f is 


Vf] = T, (3.9) 


which is the rate of change of f with position along the normal. 

From (3.8) it follows that V f can be used to find the rate of change of f in 
any direction. To find the rate of change of f in the direction of the unit vector 
u, set dr = uds where ds is the distance along u. Then df = V f - u ds and so 
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f= constant Vf 


Fig. 3.2. Grad f is perpendicular to any vector dr lying in the surface f = constant. 


df 
— = Vf. u. 1 

1. f-u (3.10) 
This is the rate of change of f in the direction of the unit vector u, and is 


called the directional derivative of f. This can also be written as 
of = |V f| cos 8, (3.11) 
ds 


where @ is the angle between Vf and the unit vector u. Since —1 < cos@ < 1, 
it follows that the magnitude of Vf is equal to the maximum rate of change 
of f with position. 

The symbol V can be interpreted as a vector differential operator, 


0 ð ð 
V= (5 By’ 2) ; (3.12) 


where the term operator means that V only has a meaning when it acts on 
some other quantity. 

The gradient has many important applications. These include finding nor- 
mals to surfaces and obtaining the rates of change of functions in any direction, 
as in the following examples. 


Example 3.4 


Find the unit normal n to the surface z? + y? — z = 0 at the point (1,1,2). 

Define f(x,y,z) = z* + y? — z = 0, so the surface is f = 0. Then Vf = 
(22, 2y, —1). At the point (1,1,2), Vf = (2,2, —1). This is a vector normal to 
the surface. To find the unit normal we need to divide by the magnitude, which 
is (2? + 2? + 1?)}!2 = 3 so n = V f/|V f| = (2/3, 2/3, -1/3). Note that the 
unit normal is not uniquely defined: the vector —n = (—2/3, —2/3,1/3) is also 
a unit normal to the surface. 
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Example 3.5 


Find the directional derivative of the scalar field f = 22+ y +2? in the direction 
of the vector (1,1,1), and evaluate this at the origin. 

The gradient of f is Vf = (2,1,2z). To find the directional derivative, we 
must take the dot product with the unit vector in the direction of (1,1, 1) which 
is u = (1,1,1)/V3. The directional derivative is then Vf -u = (3 + 2z)/¥V3. 
At the origin, z = y = z = 0, the directional derivative takes the value v3. 


3.2.1 Gradients, conservative fields and potentials 


There is a very important link between the gradient of a scalar field and the 
concept of a conservative vector field defined in Section 2.2.3. Recall that a 
conservative vector field is one in which the line integral along a curve connect- 
ing two points does not depend on the path taken. The connection between 
gradients and conservative fields is given by the following theorem. 


Theorem 3.1 


Suppose that a vector field F is related to a scalar field ¢ by F = V@ and 
V > exists everywhere in some region D. Then F is conservative within D. 
Converseiy, if F is conservative, then F can be written as the gradient of a 
scalar field, F = V¢. 


Proof 


Suppose that F = V@. Then the line integral of F along a curve C connecting 
two points A and B is 


F -dr =| V¢d-dr. 
C C 
Using (3.8) this can be written as 


f| F-ar= f dé = [$]? = ¢(B) - 4(A), 
C C 


where the line integral has been evaluated simply using ¢ as the parameter. 
Since this result only depends on the end points of C, F' is conservative. 

Conversely, suppose that F is conservative. Then a scalar field ¢(r) can be 
defined as the line integral of F from the origin to the point r: 


r 
lr) = F-dr. (3.13) 


Since F is conservative, the value of ¢ does not depend on the path taken from 
0 to r, so ¢ is well defined. From the definition of an integral it then follows 
that an infinitesimal change in ¢ is given by 
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dġ = F -dr. 


Comparing this with (3.8) shows that F : dr = Vọ dr. This must be true for 
any choice of dr and so F = V¢. 0 


If a vector field F is conservative, the corresponding scalar field @ which 
obeys F = V¢ is called the potential for F. Note that the potential is not 
unique, since an arbitrary constant can be added to ġ without affecting Vo. 
This arbitrary constant corresponds to the arbitrary choice of the origin for the 
lower limit in the integral in the definition (3.13). 


Example 3.6 
Show that the vector field F = (22 + y,z,2z) is conservative. 

F is conservative if it can be written as the gradient of a scalar field ¢. This 
gives the three equations 


a¢_ 2 


ð 
E E ðy T, ar = 2 


Ox Oz 


Integrating the first of these equations with respect to z gives ¢ = z? + ry + 
h(y, z) where h is an arbitrary function of y and z, analogous to a constant of 
integration. The second equation forces the partial derivative of h with respect 
to y to be zero, so that h only depends on z. The third equation yields dh/dz = 
2z, so h(z) = z* +c, where c is any constant. Therefore all three equations are 
satisfied by the potential function 


o=2*+2y+2" 


and F is a conservative vector field. 


3.2.2 Physical applications of the gradient 


The gradient of a scalar field appears in many physical contexts. Two examples 
are given below. 


e Let p denote the pressure within a gas. Then there is a force F acting on any 
volume elament ôV due to the pressure gradient, given by F = —Vpdv. 

e A material has a constant thermal conductivity K and a variable temperature 
T(r). Because of the temperature variation, heat flows from the hot regions 
to the cold regions. The heat flux g is a vector quantity, q = -K VT. 
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EXERCISES 


3.1 Find the gradient of the scalar field f = zyz, and evaluate it at the 
point (1,2,3). Hence find the directional derivative of f at this point 
in the direction of the vector (1, 1,0). 

3.2 Find the unit normal to the surface y = z + 23 at the point (1, 2,1). 

3.3 Show that the gradient of the scalar field ¢ = r = |r| is r/r and 
interpret this result geometrically. 

3.4 Find the angle ee the surfaces of the sphere z? + y? + 2° = 2 
and the cylinder z? + y? = 1 at a point aye they intersect. 

3.5 Find the gradient of the scalar field f = yx? + y° — y and hence find 
the minima and maxima of f. Sketch the contours f = constant and 
the vector field Vf. 

3.6 If a is a constant vector, find the gradient of f = a-r and interpret 
this result geometrically. 

3.7 Determine whether or not the vector field F = (sin y, x, 0) is conser- 
vative. 

3.8 Consider the vector field F = (y/(z? + y?), —2/(z? + y?), 0). Show 
that F can be written as the gradient of a potential ¢. Show also 
that the line integral of F around the unit circle z? + y? = 1 is 
non-zero. Explain why this result does not contradict Theorem 3.1. 


3.3 Divergence of a vector field 


This section introduces the first of two ways of differentiating a vector field, the 
divergence. The second way of differentiating a vector field, the curl, is defined 
in Section 3.4: Each of these quantities is defined in terms of an integral. 

The divergence of a vector field u is a scalar field. Its value at a point P is 
defined by 


6V—-0 
where ôV is a small volume enclosing P with surface 6S and n is the outward 
pointing normal to 6S. Physically, this corresponds to the amount of flux of 
the vector field u out of ôV divided by the volume ôV (Figure 3.3). 
As in the case of the gradient, this physical definition leads to an equiva- 
lent definition in terms of the components of u = (wu, u2,u3) in a Cartesian 


divu = lim y wp,” -ndS, (3.14) 
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Fig. 3.3. Definition of the divergence. The small volume ôV has surface ôS and 
outward normal n. 


coordinate system. To derive this alternative form, take the volume ôV to be a 
small rectangular box with sides of length dz, dy and 6z and centred on a point 
(x,y,z) (Figure 3.4). It is assumed that the components of u have continuous 
partial derivatives. 


bx 


Fig. 3.4. Rectangular box used to obtain an expression for div u in Cartesian coor- 
dinates. 


Since the rectangular box has six faces, there are six different contributions 
to the surface integral in (3.14). Consider first the contribution from the face 
labelled Sı in Figure 3.4. This face is perpendicular to the z-axis, so the unit 
outward normal is (1,0,0) and hence u:n = u1. The centre of face S$, is at the 
point (x + 6z/2,y,z) and the area of the face is dydz, so the contribution to 
the surface integral from this face is 


I u: ndS x u(x +ôxr/2,y,z)ôyðz, 
Sı 


where we have used the fact that since the surface is small, the surface integral 
can be approximated as the value of u -n multiplied by the area of the surface. 
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A similar argument can be used to approximate the contribution to the 
surface integral from S2, which is located at (x — 6z/2,y, z). The unit outward 
normal for S2 is (—1,0,0) so u : n = —u; and the contribution to the surface 
integral is 


J u: ndS x —u,(zr — 62/2, y, z) dy dz. 
S2 


Adding the contributions from these two surfaces and making use of the defi- 
nition of the partial derivative (3.1), the combined contribution to the surface 
integral is 


I u-ndS = COETDE ue - F,v,2)] dy ôz 
Si+Se 2 2 


ou ôx dy dz 
Our sy, 


Ox 


Hence the contribution to div u defined in (3.14) from surfaces Sı and S2 is 
Ou, /ðx. Note that this is now exact since the divergence is defined by taking 
the limit ôV — 0. Similarly, the contribution to div u from the two surfaces 
perpendicular to the y-axis is Qu, /Oy and that from the surfaces perpendicular 
to the z-axis is du3/0z. These are found simply by permuting the z-, y- and z- 
axes. Finally, adding all six contributions together gives the definition of divu 
in terms of the Cartesian components of u: 
: ðuı Our Ou3 
divu = ar Oy az" (3.15) 
The divergence of u can also be written in terms of the differential operator 
V defined in (3.12), since 
= ou ĝu (2 0 ð 
The form V - u will be used to indicate the divergence of u in the remainder 
of this book. 


2 


And 


) ` (u1, U2, U3) =V-u. (3.16) 


Example 3.7 


Find the divergence of the vector field u = r. 
The components of u = r are u = (x,y,z). The divergence of u is therefore 


Ox 4 Oy Oz 
ðr Oy T 


div u = =3. (3.17) 
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3.3.1 Physical interpretation of divergence 


The physical definition of the divergence (3.14) gives an intuitive meaning in 
terms of the flux of the vector field out of a small closed surface. This can 
also be interpreted as the rate of ‘expansion’ or ‘stretching’ of the vector field. 
Consider for example the simple vector field u = (2z,0,0). This vector field 
only has a component in the z direction and it is sketched in Figure 3.5(a} It 
is useful to think of vector fields as representing the motion of a gas. Figure 
3.5(a) then represents a gas which is expanding, and the divergence of u, from 
(3.15), is 1 everywhere. The vector field v = —u = (—z, 0,0) is shown in Figure 
3.5(b). This vector field is contracting, and its divergence is V -v = —1. Finally 
consider the vector field w = (0,2,0), sketched in Figure 3.5(c). This vector 
field is neither expanding nor contracting, and its divergence is zero. A vector 
field w for which V - w = 0 everywhere is said to be solenoidal. 


3.3.2 Laplacian of a scalar field 


Suppose that a scalar field ¢ is twice differentiable. Then the gradient of ¢ is a 
differentiable vector field V@, so we can take the divergence of V@ and obtain 
another scalar field. This scalar field, V - V¢ is called the Laplacian of ¢ and 
has its own symbol, V?¢, so 


V -Vo= V"¢. 


The Laplacian of ¢ is often referred to as ‘del squared ¢’. The formula for V?¢ 
can be found by combining the formulae for div (3.15) and grad (3.7), 


ð (0¢ ð (0d ô (0¢ 
2 oon aD — ———= =a -=n aum 
ee Ox (5) ss Oy (2$) as Oz & 
ð? h g? $ ð? $ 
= —+——+—., Al 
Ox? : Oy? a Oz? (re) 
Thus the Laplacian of ġ is just the sum of the second partial derivatives of ¢. 
The Laplacian can also act on a vector quantity, in which case the result is a 
vector whose components are the Laplacians of the components of the original 
vector: 
V’u = (V°u, V7u2, V7us) . (3.19) 
The Laplacian is a very important quantity, occurring in many physical 
applications including heat transfer and wave motion. These applications will 
be considered in Chapter 8. The equation V2¢ = 0 is known as Laplace’s 
equation. 
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(a) 


(b) 


(c) 


Fig. 3.5. The three vector fields (a) u = (z,0,0), (b) v = (—z, 0,0), (c) w = (0, 2,0). 
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3.4 Curl of a vector field 


The curl of a vector field u is a vector field. Its component in the direction of 
the unit vector 7 is 


1 
nm-curlu = jim, 55 f, u: dr, (3.20) 


where ôS is a small surface element perpendicular to n, dC is the closed curve 
forming the boundary of ôS and 6C and n are oriented in a right-handed sense, 
as shown in Figure 3.6. Note that this has a similar form to the definition of 
divergence (3.14), but with a line integral instead of a surface integral. 


5S \c 


ôC 


Fig. 3.6. Definition of the curl of a vector field. The small surface 6S is enclosed by 
the curve ôC and has unit normal vector n. 


To obtain an expression for curlw in terms of the components (u1, uz, u3) 
of u, choose n = e3, the unit vector in the z direction, to determine the z 
component of curlu. The surface 6S then lies in the z, y plane and can be 
chosen to be a small rectangle with sides of length dz, dy centred on the point 
(x,y,z), with area 6S = ôx dy. The right-hand rule means that the line integral 
in (3.20) must be taken in the anticlockwise direction. The line integral then 
has four sections, as shown in Figure 3.7. It is assumed that the vector field u 
is differentiable with continuous partial derivatives. 

Consider first the section Cı of the line integral, which has its centre at 
the point (x,y — dy/2,z). Here, the line integral is directed in the positive z 
direction, so u- dr = u, dz. Since the length dz is small, the contribution to 
the line integral is approximately 


/ u-dr=u(z,y — dy/2,z) dz. 
Cı 


Similarly, on C3, centred at (z,y + dy/2,z), the integral is directed in the 
negative z direction, so 
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C3 


C, oy 0(X,y,2) A Cy 


X 


Fig. 3.7. A rectangle of four line segments is used to find an expression for curl u in 
Cartesian coordinates. 


/ u -dr = —u;(z,y + dy/2, z) dz. 
C3 
Adding these two contributions together gives 
| u-dr = (ui(z,y — dy/2,z) —uy(z,y + by/2, z)) dz 
Ci+Cs 


——— dy dz. 
Oy ies 
Proceeding in a similar way with the line integrals along C2, located at 
(x +62/2,y,z) with dr directed in the positive y direction, and C4, located at 
(x — 6x/2,y,z) with dr directed in the negative y direction, we obtain 


f u:dr z% (uz(x+ôzr/2,y,z)— u(x — ôz/2,y,z)) dy 
C2+C4 


Our 


Adding together all four contributions, dividing by 6S and taking the limit 
ôS — 0 gives the z component of curl u: 
Our _ Ou) 
Oy 
The other components can be found by permuting z, y and z cyclically (x > 
Y > z => T, UW > U2 > u3 > U1), giving 


e3-curlu = (3.21) 
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eu uz ðu» Ou) Ou3 u» =) 
curlu = ( Oz’ Oz Oz’ Or Ody) (3:22) 


Notice that there is a similarity between this formula and that for the cross 
product of two vectors (1.4). Curl u can also be written in terms of a determi- 
nant, 


ej €> ©3 
curl u = 2 Z 


Ui U2 U3 


provided that the determinant is expanded so that the partial derivatives act 
on the components of u. This can also be written as the cross product of the 
differential operator V and the vector u, so 


curlu = V x u. 
The notation V x u will be used henceforth. 


Example 3.8 
The vector field u is defined by u = (zy, z + z,y). Calculate V x u and find 


the points where V x u = 0. 
The components of V x u are found using (3.22): 


Gxu = (2 -26ta (zy) _ dy Peta Aen) 


Oy Oz ' Oz ðr’ ðr Oy 
(1 — 1,0-0,1-2) = (0,0,1 — z). 


Hence V x u = 0 on the plane z = 1. 


3.4.1 Physical interpretation of curl 


From the physical definition of V x u given in (3.20) and Figure 3.6 it is clear 
that V x u is related to the rotation or twisting of the vector field u. 
Consider the three simple vector fields shown in Figure 3.5. For the first 
of these, u = (z,0,0), the vector field is expanding but there is no sense of 
rotation, and computing the curl gives V x u = 0. A vector field u for which 
V x u = 0 everywhere is said to be irrotational. Similarly, the second example, 
v = (—z, 0,0), is also irrotational. For the third example, w = (0, x, 0), Figure 
3.5(c), V x w = (0,0, 1), so there is a component of V x w in the z direction, 
out of the page. The vector field w has a rotation associated with it in the 
following sense. Think of w as the velocity of a fluid. Then a small particle 
placed in this fluid will rotate in an anticlockwise sense as it moves with the 
fluid, since at any point the velocity component in the y direction to the right 
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of the particle is greater than that on the left. This rotation is about an axis 
in the z direction, which is in the direction of V x w. The vector V x w can 
therefore be related to the rotation of a small particle placed in the velocity 
field w: the rate of rotation depends on the magnitude of V x w, and the axis 
of rotation is in the direction of V x w. In the context of the motion of a fluid, 
the curl of the velocity field is often referred to as the vorticity of the fluid. This 
relationship between rotation and curl is made more precise in the following 
section. 


3.4.2 Relation between curl and rotation 


Consider a rigid body rotating with angular velocity 2. Then, as discussed in 
Section 1.3.1, the velocity v at any point can be written as v = £2 x r, where 
the vector §2 is directed along the axis of rotation (Figure 1.12). 

Choosing the z-axis in the direction of 9 gives £2 = (0,0, 2) and hence 
v = (0,0, 82) x (z,y,z) = (—Ny, Nz,0). Computing the curl of this velocity 
field gives 

(Nz) Ay) A(Nr) NRY) _ 

ae On ade dy J- (00520); 
Hence V x v = 282, i.e. the curl of the velocity field is equal to twice the 
rotation rate. 


vxv=( 


3.4.3 Curl and conservative vector fields 


Suppose that a vector field u is related to a scalar field ¢ by u = V ọ. Recall 
from Theorem 3.1 that this means that u is conservative. Now consider the 
curl of u, 


- e (2 (2%) _ 2 (ə ð (08) _ 8 (24 
wae = C (3) dz | a G (a) Ox (a=) 
ð (d0¢ ð (0¢ 
“=. = (sr) Oy (5) 
Since the ordering of the cross derivatives of @ does not matter (see Section 
3.1.1), all of these terms cancel out, giving the result that 
Vx V¢=0. (3.23) 


Thus any vector field that can be written as the gradient of a scalar field is 
irrotational. 
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The converse of this result is also true, so that any irrotational vector field 


is conservative. This result will be proved in Section 5.2. In combination with 
Theorem 3.1, this means that the following three statements are equivalent: 


e u can be written as the gradient of a potential: u = V¢. 
¢ u is irrotational: V x u = 0. 


ə u is conservative: the line integral of u around any closed curve is zero. 


(a) (b) (c) 
Vf ~_ a ~ 
a ~ ( f; 
x \ ~> 


Fig. 3.8. Physical picture associated with (a) gradient, (b) divergence, (c) curl. 
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Summary of Chapter 3 


The gradient of a scalar field f, written grad f or Vf, is a vector field perpen- 
dicular to the surfaces f = constant, pointing in the direction of increasing 
f. with magnitude equal to the rate of change of f in this direction. 

The components of Vf are the partial derivatives of f: 


_ (Of Of of 
NIE (Z, dy’ 5): 


o The directional derivative of f in the direction of the unit vector u is Vf - wu. 
e A vector field F' is conservative if and only if it can be written as the gradient 


of a scalar field, F = Vọ. The function ¢ is called the potential for F. 
The divergence of a vector field u, written div u or V- u, is a scalar field, 
1 Our Our uz 
V-u= lim — = — + — + —. 

u= yo v$., DOR ag. age Be 
The divergence of u corresponds to the amount of stretching or expansion 
associated with u. If V -u = 0, u is said to be solenoidal. 
The Laplacian of a scalar field ¢, written V7¢, is defined as V?¢ = V . Vo. 
The curl of a vector field u is a vector field, curl u or V x u. Its component 
in the direction of a unit vector n perpendicular to the surface element 6S is 


n-Vxu= jim asf, u- dr. 


In component form, 


_ Ou3 Our Ou) Ou3 Our 3) 
vxu= (3 Oz’ Oz Or’ Or By) 


Physically, V x u corresponds to the rotation or twist of u. If V x u = 0, u 
is said to be trrotational. 

If u = Vd, then u is irrotational. 

The physical pictures corresponding to grad, div and curl are sketched in 
Figure 3.8. 
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EXERCISES 


3.9 Find the gradient V¢ and the Laplacian V?¢ for the scalar field 
b= T? + cy + y2?. 
3.10 Find the gradient and Laplacian of 


@ = sin(kz) sin(ly) exp( V k? + 122). 


3.11 Find the unit normal to the surface zy? + 2yz = 4 at the point 
(—2,2, 3). 

3.12 For $(2,y,2z) = £? +y? +z? +ry—3rz, find Vo and find the minimum 
value of ¢. 

3.13 Find the equation of the plane which is tangent to the surface z? + 
y? — 2z° = 0 at the point (1,1, 1). 

3.14 Find both the divergence and the curl of the vector fields 
(a) u = (y, 2,2); 

(b) v = (xyz, 27,2 — y). 

3.15 Show that both the divergence and the curl are linear operators, i.e. 
V -(cu+dv) = cV -u+dV v and V x (cu+dv) =cV xu+dV xv, 
where u and v are vector fields and c and d are constants. 

3.16 For what values, if any, of the constants a and b is the vector field 
u = (ycosz + azz,bsin z + z, x? + y) irrotational? 

3.17 (a) Show that u = (y?z,—z?siny + 2zyz, 2z cosy + yz) is irrota- 
tional. 

(b) Find the corresponding potential function. 
(c) Hence find the value of the line integral of u along the curve 
z =sinnt/2,y=t? -t,z=t4,0<t<1. 


4 


Suffix Notation and its Applications 


4.1 Introduction to suffix notation 


This chapter introduces a powerful, compact notation for manipulating vector 
quantities. In the previous chapters, many of the vector expressions are awk- 
ward and cumbersome. This applies particularly to those expressions involving 
the cross product and the curl, such as the scalar triple product (1.8), the 
derivation of the alternative expression for the vector triple product (1.9) and 
the demonstration that V x V¢ = 0 (3.23). Through the use of a new nota- 
tion, suffiz notation, such complicated expressions can be written much more 
concisely and many results can be proved more easily. 

In this section, some simple vector equations are written using suffix nota- 
tion. Consider first the equation c = a + b. This vector equation is equivalent 
to the three equations for the components of c, c; = a; + b; for i = 1,2,3. In 
suffix notation, the equation is simply written 


Ci =a, + bi 


and it is understood that this equation holds for 1 = 1, 2 and 3. The suffix i is 
called a ‘free suffix’. The choice of this free suffix is arbitrary, so the equation 
could equally well be written c; = aj +b; or ck = ak +b. However, for simplicity 
and clarity the suffix 7 will be used for the free suffix in a vector equation in 
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this book. Note that the same free suffix must be used for each term in the 
equation. This is an important rule which must be followed for an expression 
in suffix notation to be meaningful: the free suffix must match in each term in 
the expression. 

Now consider the dot product of two vectors, a -b = a,b, + a2b2 + agbs. 
This can be written more compactly as 


3 
a-b= > ajbj. 


j=l 
In suffix notation, this is written simply as 
a-b= a;b;, (4.1) 


where the repeated suffix j implies that the term is to be summed from j = 1 
to 7 = 3. This is known as the summation convention: whenever a suffix is 
repeated in a single term in an equation, summation from 1 to 3 is understood. 
The repeated suffix is referred to as a ‘dummy suffix’, and must appear no more 
than twice in any term in an equation. The choice of the dummy suffix does 
not matter, so we can write a : b = a;b; = axb,, since each of these expressions 
is equivalent to a,b, + agbe + ag3bz3. 

Next, suppose that an expression involves two dot products multiplied to- 
gether, (a - b)(c - d). In order to indicate which vector is dotted with which, a 
different dummy suffix must be used for each of the dot products: 


(a-b)(c- d) = ajbjerdk. 


Here, both j and k are repeated, so the summation convention implies summa- 
tion over both j and k. Again the choice of dummy suffix is arbitrary, so for 
example we could have written (a-b)(c-d) = a,bjcmdm. However, it is essential 
that no suffix appears more than twice in any term, since this would lead to 
ambiguity. 

Note that it is the suffices that indicate which vector is dotted with which, 
not the ordering of the components of the vectors. In fact, since the components 
are just multiplied together, the ordering of terms is arbitrary, so the expression 
Ckajdkbj also means (a : b)(c i d). 


Example 4.1 


Write the suffix notation expression ajb;c; in ordinary vector notation. 
The suffix j is repeated and is therefore a dummy suffix to be summed over. 


So ajbicj means 
3 


> a bic; = (a $ c)b;, 
j=l 
which is the į component of the vector (a - c)b. 


4. Suffix Notation and its Applications 67 


Example 4.2 


Write the vector equation 
u + (a: b)v = |a|*(b-v)a 


in suffix notation. 
First introduce a free suffix 2, which is understood to run from 1 to 3, to 
write the equation in component form; also write |a|? as a dot product: 


ui + (a: b)v; = (a-a)(b- v)ai. 


Now introduce a dummy suffix, which is repeated and therefore summed from 
1 to 3, for each of the dot products: 


ui + ajbjvi = ajajbkUkQi. 
Note that two different dummy suffices are used on the right-hand side to avoid 
ambiguity. 


Example 4.3 


Show that the product of two N x N matrices A and B, C = AB can be written 
in suffix notation as C;; = AikBkj. Hence show that the trace of the matrix 
AB (defined as the sum of the elements on the diagonal) is the same as the 
trace of BA. 

Ci; is the element in the ith row and jth column of the matrix C. The rule 
for matrix multiplication is that if C = AB then the element Ci; is obtained 
by taking the ith row of A and the jth column of B and multiplying these 
together term by term, so 


N 
Cij = An Bij + AnBa;...+ AinBnj = Y Aik Brj = Air Bay, 
k=1 
where the repeated index k implies the sum from 1 to N. 
The trace of the matrix C is the sum of the elements on the diagonal, 
Tr(C) = Cii + Coa... + Cnn = Cj;. 
The trace of AB is 
Tr(AB) =Tr(Ai, Byj) = Ajk Bri- 


Similarly the trace of BA is 
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Tr(BA) = ByeAgj 

A,j Bj, (since order of terms does not matter) 
A;kBg; (relabelling j + k) 

= Tr(AB). 


Note that this proof makes use of the fact that the choice of label used for a 
dummy suffix is arbitrary, so the labels 7 and k can be interchanged. 


4.2 The Kronecker delta ĝ;j 


The Kronecker delta is written 6;; and is defined by 


1 i=j, 
ui = $ if i # j. eo) 
The suffices 7 and j can each take the values 1, 2 or 3, so 6;; has nine elements. 
From the above definition it follows that three of these. are equal to 1 (61; = 
622 = 633 = 1) while the remaining six elements are equal to 0 (12 = 613 = 
621 = ĝ23 = 031: = 632 = 0). 6;; is an example of an object called a tensor. 
Tensors are described in detail in Chapter 7, but for the time being it is simplest 
to think of 6;; as the 3 x 3 identity matrix, 


10 0 
6;={ 01 0]. 
001 


From the definition it is clear that 6;; is symmetric, i.e. 6;; = 4;:. 
Consider now the expression 6,;a;. Notice that the suffix j is repeated, so 
by the summation convention, summation from j = 1 to 3 is understood. Hence 


3 
4j2; = 5 0:52; = 6;101 + dj2a2 + 65343, 
j=1 
and it is clear that the result depends on the value of ¢. If è = 1, then ô; = 1 
while dj2 = di2 = 0, so the right-hand side simplifies to a. Similarly, if i = 2 
the result is az and if i = 3 the result is a3. In other words, the right-hand side 
simplifies to a;, giving the important equation 


450; = Qj. (4.3) 
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From the symmetry of 6;; it follows that 6,;a; = aj. Because of this property 
of 4;;, it is sometimes referred to as the ‘substitution tensor’, since its effect 
when multiplied by a; is to replace the j with i. 

There is a relationship between 6;; and the dot product. The dot product of 
a and b can be written a -b = 6;;a,b;. In this expression, both the 7 and the j 
suffices are repeated, so by the summation convention, both are to be summed 
from 1 to 3, giving a total of nine terms. However, because of the definition of 
6:;, only three of these terms (the ones with i = j) are non-zero, so 


3 3 
6;; a:b; = >»: X 5:5 a:b; = aıbı -+ azb + a3b3 =a-b. 


t=1 j=1 


This result can also be demonstrated using (4.3), since 6;;b; = bi, soa@-b = 
a,b; = aidi; bj = dijaibj. 


Example 4.4 


Evaluate ĝjj. 
Since the suffix 7 is repeated, the summation convention implies that this 
expression must be summed from j = 1 to 3, so 


3 
05; = ` 05; = 61, + 22 + b33 = 3. (4.4) 
j=1 


Example 4.5 
Simplify Oi;O5k- 
Here the suffix j is repeated, and must therefore be summed over: 


3 
bigSjk = X 5ig5jk = 5:1 61k + Ôi2ô2k + izsk. 


j=1 


The result depends on the values of 2 and k. If, for example, 1 = 1 and k = 2, we 
have 1x0+0x1+0x0= 0; but ifi = l and k = 1, we have 1x1+0x0+0x0= 1. 
It is apparent that if i and k are different the result is 0, but if i and k are equal 
the result is 1. This result is therefore simply 6;, so the solution is 6;;6;, = dix. 
Notice that this result is consistent with the substitution rule described above: 
the effect of d;; on 6;, is to replace the j with i. 
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4.3 The alternating tensor €;5, 


This section introduces the quantity which is used for writing cross products in 
suffix notation. This will prove extremely useful for manipulating expressions 
involving the cross product and the curl. 

The alternating tensor is written €;;; and is defined by 


O if any of i, j, k are equal, 
Eijk = $ +1 if (i,j,k) = (1, 2,3), (2,3,1) or (3,1,2), (4.5) 
—1 if (i,j,k) = (1,3,2), (2,1,3) or (3,2,1). 
Since €;;; has three suffices, each of which can take any of the three values 1, 
2 or 3, cijk has 27 elements. However, from the above definition, all but six 
of these are zero. The six non-zero elements are €123 = €231 = €312 = 1 and 


€132 = €213 = €321 = —1. 
There are two important symmetry properties of eij which follow directly 
from its definition: 


e cijk is unchanged if the suffices are reordered by moving them to the left and 
putting the first suffix third (a cyclic permutation of the suffices), i.e. 


Eijk = €jki = €kij- (4.6) 
e The sign of cij changes if any two of the suffices are interchanged, e.g. 
€ijk = —€jik- (4.7) 
The relationship between €;;, and the cross product is as follows: 
(a x b)i = €ijkajbk. (4.8) 


In this equation, both j and k are repeated, so they are dummy suffices and 
must be summed over. To check that this agrees with the previous definition 
of a x b, consider first the case i = 1. The right-hand side is then 


3 3 
€15h0j;0% = ) ) €1j40; 04. 


j=l k=1 


Since €i; is only non-zero when all three of its suffices are different, only the 
two terms j = 2, k = 3 and j = 3, k = 2 are non-zero in the double sum. 
Hence the right-hand side reduces to €123a2b3 + €1323b2 = a2b3 = a3b2. This 
agrees with the previous definition (1.4) for the first component of the cross 


4. Suffix Notation and its Applications 71 


product. It can be seen by cyclic permutation of indices that the second and 
third components also agree. 
There is also a relation between €;;; and the determinant of a 3 x 3 matrix. 


This can be written 
|M| = eijk MıiıM2;M3k. (4.9) 

A related formula is 
Engr |M | = Eijk Mpi Mqj Mer. (4.10) 


An expression for the scalar triple product. a -b x c can be deduced in suffix 
notation as follows: 


a:-bxc= a;(b x C); = Q5€:jh05Ck = EijkaibjCk. (4.11) 


A comparison of this neat and elegant expression with the cumbersome formula 
in terms of the components (1.8) shows the power of suffix notation. The prop- 
erties of the scalar triple product can also be deduced using suffix notation, as 
in the following examples. 


Example 4.6 


Use suffix notation to show that a-b x c=a~x)b-c. 


a-bxc = ejxa;b;cx 
= €xijaibjc, (using Eijk = €Ekij ) 
= (ax b)kCk 
= axb-c. 


Example 4.7 


Show that a-bx c=b-cxa. 
The demonstration of this result is very similar: 


a-bxc = EijkaibjCk 
= EjkiâibjCk (using Eijk = Ejki) 
= bjéjkiCka; (just rearranging terms) 
= b;(cexa); 


= b-cxa. 


Example 4.8 


Evaluate €;;%. 
Since €;;, = 0 if any of i, j, k are equal, it follows that cii} = 0. 
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Example 4.9 
Evaluate €:;% €ijk- 

In this expression all three suffices i, j and k are repeated, and must there- 
fore be summed over, giving a total of 27 terms. Only six of these terms are 
nON-ZerO, SO Eijk€ijk = E123 + Eigse + C213 + a1 + €312 + €321 = 6. 


4.4 Relation between €;;, and 46;; 


An important relationship between ¢,;;, and 4,; is the following equation: 
€ijk€kim = Oi19j;m — Oimdjt- (4.12) 


This equation has four free suffices (2, 7, l and m) and therefore represents 81 
different equations! The left-hand side is summed over k, because the suffix k 
appears twice. 

The result (4.12) can be demonstrated by the following argument: since the 
three coordinate axes are equivalent, we need only consider the case 1 = 1. 
Consider now the possible values for j: 


1. If j = 1, €ijk = €11% = 0 and so the l.h.s. is zero; the r.h.s. is 61:d1m —Oi1mO11 
which is also zero since the two 6 terms cancel. 

2. If j = 2, €ijk = €124 = 0 unless k = 3, so only the k = 3 term contributes 
to the sum. When k = 3, the term €kım is zero unless l and m are 1 and 2. 
Therefore the |.h.s. takes the value +1 if l = 1 and m = 2, —1 if l = 2 and 
m = 1, and zero otherwise. Now the r.h.s. is 61:462m — 61m621. This is also 
equal to +1 when / = 1 and m = 2 (from the first term), —1 when / = 2 
and m = 1 (from the second term) and zero otherwise. 

3. If 7 = 3, an equivalent argument to the case j = 2 applies; the details are 
left to the reader. . 


Equation (4.12) is very useful for simplifying expressions involving two cross 
products. 


Example 4.10 


Derive the formula (1.9) for the expansion of the vector triple product ‘using 
suffix notation. 


(a x (b x c)); = €Eijkaj(b x c)k 


(writing the first cross product in suffix notation) 
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€:jk2;€kimOICm 

(writing the second cross product in suffix notation) 
(d:0jm — 9imdj1)a;bicm (using (4.12)) 

AmbiCm — ajbjci (using (4.3)) 

(a -c)b; — (a- b)c;, 


so we have shown that 
a x (bx c) = (a - c)b — (a: bic. (4.13) 


Notice one very important point in this analysis: in the second line, the k 
component of the vector b x c is required. When this is written down in suffix 
notation it is essential that ‘new’ suffices are used (l and m above) to avoid 
repeating the existing suffices 1 and j — recall the essential rule of suffix notation 
that, no suffix may appear more than twice. 


EXERCISES 


4.1 Write the vector equation a x b+ (a - d)c = e in suffix notation. 

4.2 Translate the suffix notation equation 6;;c; + €% 31046; = diemcibicm 
into ordinary vector notation. 

4.3 Use suffix notation to show that a x b = —b xa. 

4.4 Simplify the suffix notation expressions 
(a) ĝijEijk; 

(b) Eijk Eilm; 
(C) EijkEijm; 
(d) Eijk Eijk- 

4.5 Using suffix notation, find an alternative expression (involving no 
cross products) for a x b-c xd. 

4.6 If A and B are two N x N matrices, show that (AB)? = BT AT, 
where AT is the transpose of A defined by interchanging the rows 
and columns of A. 

4.7 Verify the formulae (4.9) and (4.10) for the determinant of a 3 x 3 
matrix. 

4.8 Use the formula (4.10) for the determinant of a 3 x 3 matrix M to 
show that 
(a) 6|M| = epqr€igk MpiMq; Mrk ; 

(b) |M7| = |M| ; 
(c) |MN|=|M||N|. 
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4.5 Grad, div and curl in suffix notation 


The differential operators grad, div and curl can be written using suffix nota- 
tion. To do this, the Cartesian coordinates (x,y, z) will be relabelled (z1, £2, 23). 
As in the previous section, the use of suffix notation results in a much more 
compact formulation and simplifies many of the computations. 

Consider first the gradient of a scalar field, Vf. This is defined by 


_ f{ Of of of 
ME (5: Or.’ E) 


The i component of V f is equal to the partial derivative of f with respect to 
z;, SO in suffix notation this can be written 


[V f]; = 5 (4.14) 


Thus the vector differential operator V defined in (3.12) can be written in 
suffix notation as 


[V]; = 5 (4.15) 
The divergence of a vector field u is 


Y. u = "u Ow ĝus _ uj 


Da ona On. Ba; (#140) 


where the summation convention implies the sum over j from 1 to 3. Note that 
the same expression results from taking the dot product of V defined by (4.15) 
with the vector u using the suffix notation formula (4.1) for the dot product 
of two vectors. 

The first component of V x u is, from (3.22), 


V xu], = 5 - 5 = 


where the repeated 7 and k imply a double sum, and so the suffix notation 


expression for V x u is 
Our 


This can also be obtained simply by taking the cross product of V with u using 
(4.8). 


[V x ul, = Cijk (4.17) 
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Example 4.11 
Let r be the position vector r = (£1, 22,23) and r = |r|. Use suffix notation to 
evaluate Oz;/Oz;. Hence find Vr, V -r and V xr. 

r = (2), 22,23), so in suffix notation, r; = z;. The three coordinate axes 
Z1,22,23 are independent. Thus the derivative of each of the z; with respect 
to one of the others is zero, while the derivative with respect to itself is 1. Thus 


Oz; a. 1 ift= Í, 
ðr; \0 ifižj 
= 0:3. (4.18) 


To find Vr, first write r = |r| = (r - r)!/? = (2;2;)!/, so 


ð ! 
[Vr]: Jz, (ita) 


1 R s, 
= z (2573) H? gg, (2525) 


1 Oz; 
= or 22i z, 


1 Ti 
= ~ Tjðij = e 
So 
Vr=r/r (4.19) 


as was shown in Exercise 3.3 without the use of suffix notation. Similarly, the 
divergence of r is 


Ox; 


and the curl of r is 


» 


Oz 
[V x r]; = cij e— = EijkÔjk = 0, 


Oz; 


using the result of Exercise 4.4(a). 


o -m a m“ me Coa 
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4.6 Combinations of grad, div and curl 


The operators grad, div and curl can be combined together in several different 
ways. Some of these combinations can be simplified or expanded into alternative 
expressions. These combinations are considered below, making use of suffix 
notation. All of these results can also be obtained without suffix notation, by 
writing out all the components, but in most cases the suffix notation method 
is much quicker. 


e Div grad: 
een O _ 0 (d0f\_ Ëf _ 
RV Oz; (v7) ~ Oz; (5) — ðzjðz; — wJ ey) 
This is the Laplacian of f introduced in Section 3.3.2, the sum of the second 
partial derivatives of f. 


e Curl grad: this combination was shown to be zero in Section 3.4.3. This result 
can be shown using suffix notation as follows: 


gw 
“ue Oz; Or; 
ð Of 
Eikj Oz, ðr; 
ein ge 5S (using Cikj = —€i;k) 
ð ðf 
Sik Dr; Orn 


= JQ, 


[V x (Vf); = 


(relabelling j + k) 


(as order of derivatives does not matter) 


since the expression has been manipulated to give minus itself. 
e Grad div: 


PRR (07), pa 
V wl], = 5 ( Ft) = Se (4.21) 


This quantity cannot be simplified further. 
e Div curl: 

ð Our 
az; I ax; 
a Ô Out 

a Oz; Oz; 

ð Ouk 


—e€E . ik ——a 
a Oz; Or; 


V-(V xu) 


(relabelling i + 7) 
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ak oTi Oz; 
= QO, (4.22) 


using exactly the same argument as for curl grad. 
e Curl curl: 


Um 
Eijk H — Ekim a 
t ar; | Oxy 


€: L€ Pum 
ijkEklm ðz;ðTı 

= (dindjm = dim jt) 
Au; Pi 07 u; 

Oz; O02; 02,02; 


= [V(V-u)-V7u],. (4.23) 


[V x (V x w)], 


um 
ðr j Oz; 


This result can be used to provide a physical definition for V? applied to 
a vector. The previous definition (3.19) was only defined in terms of the 
components of the vector in Cartesian coordinates. From the above result, 
Vu can be defined by 


V-u=V(V -u)-Vx(V xu). (4.24) 


These five combinations of grad, div and curl are the only ones that make 
sense. For example, the combination grad curl has no meaning since curl is a 
vector but grad can only act on a scalar. Combinations of three or more of the 
operators grad, div and curl can be evaluated using the above results, as in the 
following example. 


Example 4.12 
Show that 
V x (V°u) = V2(V x u). (4.25) 
Using the result (4.24), 
Vx(V7u) = Vx (V(V -u)- V x(V x u)) 


-V x (V x(V xu)) (since V x V =0) 
= -V(V. (V xu))+V’(V xu) (using (4.23)) 
= V?(V xu) (since V.V x u= 0). 


So the operators V? and Vx commute. 
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4.7 Grad, div and curl applied to products of 
functions 


Another useful application of suffix notation is in computing the action of grad, 
div and curl on products of vector and scalar fields. As in the previous section, 
these results can also be obtained by writing out all the components, but the 
suffix notation method is much more compact and elegant. Some of these results 
are straightforward applications of the usual rule for the differentiation of a 
product and can simply be written down without any calculation, but many of 
them are not so obvious. 

In the following, f and g are differentiable scalar fields and u and v are 
differentiable vector fields. 


(Vifg9)], = n (f9) = fe tose =(fVg+9Vf];, so 


Vifg)=fVo+g9VF. (4.26) 


V (fu) Aa jf ui 


= Vf-u+fV-u. (4.27) 


[V x (fu): = Eijk aT. Fay (Fu) 


ry, — Ou, 
Tj 


= Wh x. xut fV x uli. (4.28) 


ð 
V . (u x v) Jz, (Eik Yak) 


Ou; Ou, 
~ Stik gp MET ask BT 
3 


_ Ou; Ou, 
= Ekij a Oz; Uk = “ik Or, Uj 


= (V xu) v-(V xv)-u. (4.29) 
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[V x (u x v)]; 


ð 
Eijk oe (€ximUtUm) 


ð 
se (diudjm a a gy7 (UYm) 


O Oe Oe Ot 
7 ‘Oz; "Oz; 


= [u(V-v)+u-Vu-—u-Vuv—v(V-u)];, (4.30) 


where the operator u - V is defined by 


u: V (4.31) 


Oz; 
and can act on either a scalar or a vector. 
To find an expansion for the expression V (u.v), consider first the quantity 


Um 
[u x (V xv); = Eijk UjEkim a, 


ðv 
(b14jm — Ôimôji Juj Da, 
u i y OME 
70x, "On; 
Similarly, interchanging u and v, 


Ou; Ou; 


[v x (V xX u)); = "i Ba; -Yi pr; 


Adding these two equations gives 
u 2i — y OM gy DUE {Dti 

I Oa; ? Oz; 1 Ox; 3 Oz; 
[(V(u-v)-—u:-Vu-—v: Vul;. (4.32) 


[lu x(V xv)+v x(V xu); 


This can be rearranged to give 
V(u-v)=ux(V xv)+v x(V xu)+u. Vv+v: Vu. (4.33) 


The effect of applying grad, div or curl to products of more than two scalar 
or vector functions can be obtained either by the repeated use of the above 
results, or directly by suffix notation, as in the following example. 
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Example 4.13 
Find an expansion for V - (fgu). 


Vifg):u+(fg)V-u_ (using (4.27)) 
(fVg+gVf)-u+(fg)V-u 
f[Vg-ut+gVf-u+fgV-u. 


V -(fgu) 


Alternatively, using suffix notation, 


V -(fgu) ae 


om + af 


= Phy Fae Che utoVf- u. 


Example 4.14 


Show that u - Vu = V(|u|?/2) — u x (V x u). 
Apply (4.33) with v = u: 


V(u - u) = 2u x (V x u) + 2u . Vu. 
Rearranging this and dividing by 2 gives 
u: Vu = V(ju|?/2)— u x (V x u). (4.34) 


Example 4.15 


Use the results (4.30) and (4.33) to provide a definition of u - Vv that is not 
given in terms of Cartesian components. 
By subtracting (4.30) from (4.33), v - Vu is eliminated and we obtain 


V(u-v)—-Vx(uxv) =ux(Vxv)tux(V xu)+2u:-Vu—u(V-v)+0(V-u) 
which can be rearranged to give 
u: Vv = 5(V(u: v) - V x (u x v) - u x (Y x v) — v x (Y x u) 
+u(V -v)-v(V - u)). (4.35) 
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Summary of Chapter 4 


Suffix notation 


Suffix notation is a powerful tool for manipulating expressions involving vectors. 
The rules of suffix notation are as follows: 


Within any term in an equation, any suffix must appear either once or twice. 
No suffix may appear more that twice. | 

A suffix that appears once in any term is called a ‘free’ suffix. A free suffix 
takes the values 1, 2 and 3 and represents the components of a vector. For 
example a + b = c — d is written in suffix notation as a; + b; = c; — dj. 

In a vector equation, the free suffix must be the same in each term. The 
above equation may also be written a; + b; = cj — dj, or any other suffix 
may be used, provided the same suffix appears in each term. 

A suffix that appears twice in a term is called a ‘dummy’ suffix and is summed 
from 1 to 3. This is known as the summation convention. So ajb; means 
a,b; + a2b2 + a3b3 =a - b. 

A pair of dummy suffices can be changed. For example, a;b;, a,b; and ambm 
are all equal to a:b. 


e The order of terms in a suffix notation expression does not matter. 
e The Kronecker delta is defined by 6;; = 1 if i = j, 0 if i # j. Properties 


include Oi; = Oji, 0;;0; = a; and 0;;a;b; =a:b. 

The alternating tensor cijų is defined by €;;, = 0 if any of i, j, k are equal, 
€123 = €231 = €312 = 1, €132 = €213 = €321 = —1. Properties include Eijk = 
Cjki = Ckij, Cijk = —€jik- 


e The cross product of a and b can be written [a x b]; = €:jn0;b,- 
® Eijk and Oi; are related by CigkEklm = Oit0jm = dim Ojt- 
e Grad, div and curl can be written in suffix notation as follows: 


j (V x ul; = a a 


Of _ Ou; 
NU ae Ox; 


[V f]; = az,’ 


Combinations of operators and derivatives of products 


e V- (Vf) =V?f. 

eV x(Vf) =90. 

e V -(V xu) =0. 
eVx(Vxu)=V(V-u) — Vu. 


Co 
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Vifg) = fVg+9Vf. 


-(fu)=Vf-ut+fV-u. 


Vx(fu)=Vfxut+fV xu. 

V -(uxv)=(Vxu)-v-(Vxv)-u. 
Vx(uxv)=u(V-v)+u:-Vu-—u-Vv—-v(V u). 
V(u-v)=ux(Vxv)tux(Vxu)t+u- Vv+v. Vu. 


EXERCISES 


4.9 Write in suffix notation the vector equation a x b+ c = (a-b)b —d. 

4.10 Simplify the suffix notation expressions 
(a) dijójkôki; 
(b) EijkEklmEmni: 

4.11 Simplify the suffix notation expression ô;ja;bıCkķôn and write the re- 
sult in vector form. 

4.12 (a) Show that V x (f Vf) = 0. 
(b) Evaluate V - (fV f). 

4.13 Show that the vector u = V f x Vg is solenoidal. 

4.14 Verify the formula (4.34) for u : Vu by using (4.35). 

4.15 Show that V- V?u = V?V-u, 
(a) using suffix notation; 
(b) using (4.24). 

4.16 The vector fields u and w and the scalar field ¢ are related by the 
equation 

u+ V xw= Vọ+V’°u, 

and u is solenoidal. Show that @ obeys Laplace’s equation. 

4.17 Show that V f(r) = f'(r)r/r, where r is the position vector r = 
(z1, z2, £3) and r = |r|. 

4.18 The vector field u is defined by u = h(r)r, where h(r) is an arbitrary 
differentiable function. 
(a) Show that V x u = 0. 
(b) If V . u = 0, find the differential equation satisfied by h. 
(c) Solve this differential equation. 

4.19 A vector field u with the property that u = cV x u, where c is a 
constant, is called a Beltrami field. 
(a) Show that a Beltrami field is solenoidal. 
(b) Show that the curl of a Beltrami field is a Beltrami field. 
(c) A Beltrami field has the form u = (siny, f, g). Find the functions 
f and g and the possible values of c if it is given that g does not 
depend on z. 


3 


Integral Theorems 


This chapter describes two important theorems that link the material in Chap- 
ter 2 on line, surface and volume integrals with the definitions of the divergence 
and curl from Chapter 3. These theorems have great physical significance and 
are widely used in deriving mathematical equations representing physical laws. 


5.1 Divergence theorem 


Let u be a continuously differentiable vector field, defined in a volume V. Let S 
be the closed surface forming the boundary of V and let n be the unit outward 
normal to S. Then the divergence theorem states that 


J|, wav = H unas (5.1) 
Proof 


The volume V is divided into a large number of small subvolumes ôV; with 
surfaces ôS;, as shown in Figure 5.1. The proof of the divergence theorem then 
follows naturally from the physical definition of the divergence in terms of a 
surface integral (3.14). Within each of the subvolumes, V - u is defined by 
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Fig. 5.1. Division of a volume V into small subvolumes 6V for the proof of the 
divergence theorem. 


1 
Vite y: E (5.2) 


where the approximation becomes exact in the limit ôV; — 0. Now multiply 
both sides of (5.2) by ôV; and add the contributions from all the subvolumes: 


Evu yd Gg u-nds. (5.3) 
i i ôS; 


Now take the limit 6V; — 0. The l.h.s. becomes the volume integral of V -u over 
the volume V; this is just the definition of the volume integral. To simplify the 
r.h.s. consider two adjacent volume elements ôV and Vz (Figure 5.2). Since 


Fig. 5.2. Enlargement of two adjacent volume elements. 


the normal vector to each surface points outward, the normal vectors to the two 
surfaces along their common surface point in opposite directions: ny = —-nz2. 
Therefore the values of u-n cancel along the common surface: u-n; +u:ny = 0. 
This means that all the contributions to the sum on the r.h.s. of (5.3) from the 
interior of the region V cancel out, leaving only the surface integral over the 
exterior surface S. So in the limit dV; — 0, (5.3) becomes 
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[Jf S uav = ff unas, 
V S 


The divergence theorem is sometimes referred to as Gauss’s theorem. It has 
many important applications in physics, and it is important to develop a physi- 
cal intuition for the meaning of the theorem. Roughly speaking, the divergence 
theorem states that the total amount of expansion of u within the volume V 
is equal to the flux of u out of the surface S. This is essentially a conservation 
law, and the mathematical form of many physical conservation laws is derived 
from the divergence theorem. An example is given in the following section. 


0 


5.1.1 Conservation of mass for a fluid 


As an example of the application of the divergence theorem, this section 
presents the derivation of the law of conservation of mass for a fluid of variable 
density. 

Consider a fluid with density p(r,t) flowing with velocity u(r,t). Let V 
be an arbitrary volume fixed in space, with surface S and outward normal n 
(Figure 5.3). Then the total mass of the fluid contained in V is the volume 


integral of p: 
Mass of fluid in V = I p dv. (5.4) 
V 
Now the rate at which mass enters V is equal to the surface integral of the flux 
pu: 


™ ON 
cada i 


Fig. 5.3. Fluid flows with velocity u through a region V. 
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Rate of mass flow into V = - ff pu: nds, (5.5) 
S 


where the minus sign appears because n points outward, so mass enters V if 
u:n <0. 

We can now apply the physical law that mass is conserved: the rate of change 
of the mass in V must equal the rate at which mass enters V. Mathematically 


this is J 
s| ewv=-§ pu-nds. (5.6) 
dt JJ Jy s 


The surface integral on the r.h.s. can now be written as a volume integral using 
the divergence theorem. Also, the order of the derivative and the integral on 
the 1l.h.s. can be interchanged: 


Ils 9° wy =- - J| 5 ay, (5.7) 


where the time derivative has become a partial derivative since p is a function 
of space and time. These two integrals can now be combined into one: 


Whe P EN (pu) dV =0. (5.8) 


Now this result has been obtained without any restrictions on the volume V. 
Thus it is true for any arbitrary volume V. The only way that this can be true 
is if the integrand (the quantity inside the integral) is zero everywhere. If there 
were some point where the integrand were non-zero, a small volume could be 
drawn around that point, which would contradict (5.8). 

Therefore the law for conservation of mass of a fluid is 


ĉe + V- (pu) = (5.9) 


This conservation law takes the following form: the rate of change of the density 
plus the divergence of the flux is zero. Many other conservation laws can also 
be written in this form, for example conservation of energy or conservation of 
electric charge. 

By expanding the divergence of pu, (5.9) can be written in the form 


Phu: Vp+pV-u=0. (5.10) 


If the density of the fluid is constant and uniform, i.e. independent of time and 
space, then this equation simplifies to 


V-u=0. (5.11) 


A fluid obeying (5.11) is said to be incompressible. 
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5.1.2 Applications of the divergence theorem 


The divergence theorem has many important applications, in addition to the 
derivation of the mathematical form of conservation laws shown in the previous 
section. It can be used to simplify the evaluation of integrals, by converting a 
complicated volume integral into a simpler surface integral or vice versa. It can 
also be used to prove some important results, such as the uniqueness of the 
solution to Laplace’s equation 


V7o=0. (5.12) 
Some of these applications are illustrated in the following examples. 


Example 5.1 


Show that for any closed surface S, 


Gf (7 xu) nas =0. 
S 


Using the divergence theorem, the surface integral can be converted into a 


volume integral: 
p (V xu) -nas = ||) V-(V xu)dV. 
S V 


Since the combination div curl is always zero, this integral is zero. 


Example 5.2 


Find the relationship between the surface integral 


Pr nas 
S 


and the volume V contained within the closed surface S. 
Applying the divergence theorem, 


Pr ras= fff v-rav = fff sav =sv 


using the result (3.17) that V -r = 3. Thus the surface integral is three times 
the volume V. 


Example 5.3 


The scalar field @ obeys Laplace’s equation (5.12) in the region V and obeys 
@ = 0 on the surface S that encloses V. Show that the only possible solution 
for ¢ is ġ = 0 everywhere within V. 
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V7¢=0 => 6V?¢=0 > V-(6V4)- Vb: VG=0, 


using (4.27). Now integrate over the volume V and use the divergence theorem 
to convert the first term to a surface integral: 


ff ove-nas- fff wor av =o 


Since ¢ = 0 on S, the surface integral vanishes. The quantity |V¢@|? is always 
greater than or equal to zero, so its integral can only be zero if Vd = 0 
everywhere. This means that ¢ must be a constant, and since ¢ = 0 on S, this 
constant must be zero, so ¢ = 0 everywhere within V. 


Example 5.4 


The scalar field ¢ obeys Laplace’s equation (5.12) in the region V and the 
value of ¢ is given on the surface S that encloses V. Show that the solution to 
Laplace’s equation is unique. 

To prove uniqueness, suppose that there are two different solutions, ¢; and 
$2, obeying V?¢ı = 0 and V2¢_ = 0 in V. Since the value of ¢ is specified on 
S, d1 = ¢2 on S. Now consider the function Y = ¢, — ¢2. This function also 
obeys Laplace’s equation, since V?(¢, — ¢2) = V2d; — V2¢2 = 0. Moreover, 
p = 0 on S since ¢; = ¢2 on S. Now we can apply the result of Example 5.3 
to w: the only solution to V7 = 0 in V, » = 0 on S is 4% = 0 everywhere. 
Therefore ¢, = ¢2 everywhere, so the solution is unique. l 


5.1.3 Related theorems linking surface and volume 
integrals 


There are several other relationships between surface and volume integrals that 
can be derived from the divergence theorem by making different choices for the 
vector u. 


e Choose u = af, where a is a constant vector and f is a scalar field. Then 
V-u=fV-ata-Vf =a-Vf since a is constant. Applying the divergence 


theorem gives 
If a-Vf av = ff af- nds. 
V S 


Since a is constant, it can be taken out of the integrals: 


a-(fff vrav- ff mas) =o 


Now since a is an arbitrary constant vector, this holds for any a. This can 
only be true if the vector quantity within the large brackets is zero (for 
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example, choosing a = €1, €2, €3 in turn shows that each of the components 
of the vector in large brackets is zero). The resulting integral theorem is then 


uh Vidv= ff inas. (5.13) 


e Choose u = a x v, where a is a constant vector and v is a vector field. Then 
V-u=(V xa)-v—-—(V xv)-a=-(V x v)-a. The divergence theorem 


gives 
J -(V xv):a dV = ff axv-ndS= ff a-v x nas, 
V S S 


using the rule that the dot and cross may be interchanged in a scalar triple 
product. As in the previous example, the dot product with a can be taken 
out of the integral sign and then cancelled, giving 


J, -V xvdV = {fs x n dS. (5.14) 


e Choose u = fVg, where f and g are two scalar fields. Then V -u = Vf. 
Vg + fV7g, and the divergence theorem gives 


I Vf -V9a+fV79 dV = ff Iyo -nds. (5.15) 


This result is known as Green’s First Identity. 
e Choose u = fVg—gVf. By interchanging f and g in (5.15) and subtracting, 
we obtain 


| [V9 - Vf dV = fi] ({Vg-gVf)- ndS, (5.16) 
V S 


which is known as Green’s Second Identity. 


Historical note 


George Green (1793-1841) was a Nottingham miller who spent less than two 
years at school and learnt his mathematics by studying library books. In 1828 
he published privately his first and greatest work, ‘An essay on the application 
of mathematical analysis to the theories of electricity and magnetism’, which 
includes the two theorems above. As with many geniuses his work was not 
appreciated until several years after his death. Green’s mill in Nottingham 
has now been restored and is open to the public along with a Science Centre 
illustrating some of the applications of his work. 
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EXERCISES 


5.1 Use the divergence theorem to evaluate the surface integral 


Punas 
S 


where u = (zsiny,cos? x,y? — zsiny) and S is the surface of the 
sphere r? +y? + (z — 2)? =1. 

5.2 Verify the divergence theorem, by calculating both the volume in- 
tegral and the surface integral, for the vector field u = (y, z, z — x) 
and the volume V given by the unit cube 0 < z,y, z,< 1. 

5.3 An incompressible fluid is contained within a volume V with surface 
S and u-n = 0 on S. Using the divergence theorem, show that 


[[[ +: vew =0 


for any differentiable scalar field @. 
5.4 Two scalar fields f and g are related by Poisson’s equation, V? f = g. 


Show that 
[[[sav = 7t -nas 


5.0 Use the divergence theorem to evaluate the surface integral 


ET 


where v = (x + y,27,2”) and S is the surface of the hemisphere 
T? +y? +2? = 1 with z > 0 and n is the upward-pointing normal. 
Note that the surface S is not closed. 

5.6 Following the argument of Section 5.1.1, obtain the equation for 
conservation of electric charge relating the charge density q and the 
electric current density 3. 

5.7 Use (5.13) to obtain a definition for Vf as the limit of an integral, 
similar to the definitions of div and curl. 
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5.2 Stokes’s theorem 


Stokes’s theorem gives an alternative expression for the surface integral of the 
curl of a vector field. This is analogous to the divergence theorem, so Stokes’s 
theorem could be referred to as the ‘curl theorem’. The proof of the theorem is 
very similar to that for the divergence theorem, being based on the definition 
of cur! in terms of a line integral. 

Let C be a closed curve which forms the boundary of a surface S. Then for 
a continuously differentiable vector field u, Stokes’s theorem states that 


J| xu- nas= f$ u-dr, (5.17) 
S C 


where the direction of the line integral around C and the normal n are oriented 
in a right-handed sense (Figure 5.4). 


Fig. 5.4. Orientation of the curve C and the surface S for Stokes’s theorem. 


Proof 


To demonstrate the theorem we first divide the surface S into small pieces each 
with area ôS; and bounding curves ôC; (Figure 5.5). Within each piece of the 
surface, the definition (3.20) of V x u is 


where the approximation is exact in the limit 6S; — 0. Multiplying by 6S; and 
adding the contributions from all the surface elements, 
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Fig. 5.5. Division of the surface S into small elements 6S for the proof of Stokes’s 
theorem. 


Vx unde Df udr 


Now consider the limit 6S; — 0. The |.h.s. gives the surface integral of V x 
u: n over the surface S. On the r.h.s. the contributions to the line integrals 
from neighbouring elements cancel out, because the line elements dr point in 
opposite directions (Figure 5.6). Therefore only the curves that form part of C 


{e 5C, 
! 5S, ôS, 


Fig. 5.6. Line integrals along adjoining elements cancel out. 


contribute to the sum, so the sum simplifies to the line integral around C: 


J| 5 >u nas = f w-dr 
S Cc 
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5.2.1 Applications of Stokes’s theorem 


Stokes’s theorem can be useful for evaluating integrals, by converting line in- 
tegrals to surface integrals or vice versa. It can also be used to prove other 
theorems, as in Example 5.5 below, or to formulate physical laws (Example 
5.6). 


Example 5.5 


Show that any irrotational vector field is conservative. 
Suppose that u is irrotational, so V x u = 0. Then for any closed curve C, 


pu-dr= |] vxu-nds=o 
C S 


where S is any surface spanning C’. Thus u is a conservative vector field. Note 
that this result completes the demonstration of the statement in Section 3.4.3 
of the equivalence of the three properties (i) u = Vẹọ, (ii) V x u = 0, (iii) u is 
conservative. 


Example 5.6 


Ampère’s law states that the total flux of electric current flowing through a 
loop is proportional to the line integral of the magnetic field around the loop. 
Use Stokes’s theorem to obtain an alternative form of this law that does not 
involve any integrals. 

Let B be the magnetic field strength and j be the current density. The 
constant of proportionality is zo in SI units. Then Ampère’s law states that 


$ B-dr =p || j-nas 
C S 


for any surface S that spans the loop C’. Using Stokes’s theorem to transform 


the l.h.s. gives 
J| 5> B-ndS=m ff §-nas. 
S S 


Now if this is true for any loop C', and so any surface S, it follows that 
V x B= mĵ. 


Note the similarity between this argument and that used when applying the 
divergence theorem to the conservation of mass of a fluid in Section 5.1.1. 
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Example 5.7 


Use Stokes’s theorem to show that for any closed surface S, 


fi] (Y xu) -ndS =0. 
S 


Consider the case where a small hole is made in the closed surface. Then 
by Stokes’s theorem, the surface integral of (V x u)-n over the surface S is 
equal to the line integral of u - dr around the perimeter of the small hole. As 
the size of the hole shrinks to zero, so does the value of the line integral, giving 
the required result. Note that this result was obtained using the divergence 
theorem in Example 5.1. 


Example 5.8 


The surface S is defined by z? + 4y? = 1, —1 < z < 1. Use Stokes’s theorem to 
evaluate the surface integral 


I (xz?, —yz?,0) -n dS. 
S 


Note that this surface is not simply connected, but Stokes’s theorem can still 
be applied. By imagining a cut in the surface (Figure 5.7), the surface integral 
is equal to the sum of two line integrals around the two elliptical curves Ci 
and C» that form the ends of the cylindrical surface. In order to apply Stokes’s 
theorem, the vector field (xz?, —yz?,0) must be written as the curl of another 
vector field u. Seeking a solution of the form u = (0,0, A(z, y, z)), this can be 


Fig. 5.7. Stokes’s theorem can be used to transform the surface integral over the 
curved surface of the cylinder into two line integrals around the ends of the cylinder, 
by introducing a cut in the surface. The two line integrals along the cut cancel out. 
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achieved if 


= Oy’ -yz = Ax 


and these two equations are satisfied if h = zyz*. Now since the two curves Ci 
and C2 lie in planes z = constant, so dr = (dz, dy,0), u: dr =0 on C and C2 
so the value of the integral is zero. 


5.2.2 Related theorems linking line and surface integrals 


As in the case of the divergence theorem, Stokes’s theorem can be used to derive 
other theorems that relate line integrals to surface integrals by appropriate 
choices of the vector field u. 


e Choose u = af, where a is a constant vector and f is a scalar field. Then 
Vxu=Vfxa+fVxa=Vf xa, so Stokes’s theorem gives 


[[vixanas=$ af-dr 


Using the rules for manipulating the scalar triple product and taking out the 
constant vector a from the integrals gives 


a: (J[.-vtxnas) =a: ($ sar) 


As in Section 5.1.3, the constant a can be cancelled, giving 


J|.-5t*nas =$ sar, 


e Choose u = a x v, where a is a constant vector and v is a vector field. This 
case is more complicated, but provides a good example for the use of suffix 
notation. The line integral in Stokes’s theorem is 


f axv-dr=a-$ vxdr, 
C C 


interchanging the dot and the cross and taking the constant a outside the 
integral. Since V x (a x v) = a(V . v) — a : Vu, from formula (4.30), the 


surface integral is 
Ov; Ov; l 
J. (« Ore ag) ae 


Ov, 
[fo ie - Fim) so 


| (a(V -v) —a-Vv)-nds 
S 


96 


Vector Calculus 


where 7 and k have been interchanged in the second term. Applying Stokes’s 
theorem and cancelling the aj gives a relationship between a surface integral 
and a line integral in suffix notation: 


If. SE nj — Ja me dS = g ox dr (5.18) 


To obtain the form of this equation in vector notation, consider the quantity 
[(n x V) x v];. In suffix notation this is 
Ov; 
iklEkmnnm a — U = (dtm9in — O1ndjm)2m=— 
CjklEkmnlm Az, l ( lm9jn In jm) m Or 
Ov; Ov, 


= NZ nj. 
Oz; 1 Oz, 


This is now minus the quantity appearing in the surface integral (5.18), so 
the vector form of (5.18) is 


[[ -xv) xvas = $ v xar. (5.19) 


Note that suffix notation is the only secure method for obtaining such results, 
short of writing out all the components of the vector quantities longhand. 
Attempts to expand using the rules for a vector triple product generally give 
incorrect results. 
Choose the surface S to be a flat surface lying in the z, y plane, so n = 
(0,0,1) and choose u = (F(z, y), G(z, y),0). Then 

ðG OF ) 


V xu= (0,055 - Oy 


and u : dr = F dr + G dy. Then Stokes’s theorem gives 


[E-E deaf rasca 6 


a result known as Green’s theorem. 
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Summary of Chapter 5 


e The divergence theorem states that 


| V-ud = ff u-nds 
V S 


where S is the surface enclosing the volume V and n is the outward-pointing 
unit normal vector. 

e Geometrically, the divergence theorem follows naturally from the physical 
definition of the divergence. 

e The divergence theorem has many applications, including simplifying the 
evaluation of surface or volume integrals, deriving physical conservation laws 
and showing that Laplace’s equation has a unique solution. 

e Stokes’s theorem states that 


[fv xunds= $ war, 
S C 


where the curve C encloses the surface S and C and n are oriented in a 
right-handed sense. 

e Stokes’s theorem follows from the definition of the curl in terms of a line 
integral. 

e A number of other related theorems linking volume, surface and line integrals 
can be derived from the divergence theorem and Stokes’s theorem. 
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EXERCISES 


5.8 Show that 


$ r-dr=0 

C 
for any closed curve C. 

5.9 Verify Stokes’s theorem by evaluating both the line and surface in- 
tegrals for the vector field u = (22 — y, —y?, —y?z) and the surface 
S given by the disk z = 0, r? +y? < 1. 

5.10 Use Stokes’s theorem to show that 


$ f¥o-dr=—§ oVf-dr 


for any closed curve C and differentiable scalar fields f and g. 
5.11 If u is irrotational, express the surface integral 


|| ux Vf-ndS 
S 
as a line integral. 


5.12 The magnetic field B in an electrically conducting fluid moving with 
velocity u obeys the magnetic induction equation 


Show that the total flux of magnetic field through a surface enclosed 
by a streamline of the flow (a closed curve which is everywhere par- 
allel to u) is independent of time. 
5.13 Use (5.18) to show that the area A of a flat surface S enclosed by a 
curve C is 
A = 1/2 $ r xdr 
C 


6 


Curvilinear Coordinates 


6.1 Orthogonal curvilinear coordinates 


So far in this book we have used rectangular Cartesian coordinates. In many 
physical problems, however, these are not the most convenient coordinates to 
use. Consider, for example, the problem of finding the electric field produced by 
a charged sphere. In this chapter the general theory of non-Cartesian coordinate 
systems is introduced. Formulae for grad, div and curl in these coordinate 
systems are developed and the two most important examples, cylindrical and 
spherical polar coordinates are described. 

Suppose a transformation is carried out from a Cartesian coordinate system 
(x1,Z2,23) to another coordinate system (u1, u2,u3). This new system will 
be called a curvilinear coordinate system. It will be assumed that there is a 
one-to-one relationship between the z; and the u;, so that for example zı 
can be written as a function of the u;i, zt} = 2,(uj,U2,u3) and conversely 
u, = U1 (T1, T2, T3). 

The surfaces u; = constant are referred to as coordinate surfaces and the 
intersection of these surfaces defines the coordinate curves, so for example the 
u; coordinate curve is the intersection of the surfaces u> = constant and u3 = 
constant (Figure 6.1). 
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Fig. 6.1. The coordinate surfaces and coordinate curves of a curvilinear coordinate 
system. 


Now consider a small displacement dz = (dz,,dz2,dx3). Since the z; are 
functions of the u; this can be written as 


Ox Oz Ox 
dx = —d —d — duz, 
Bu, w+ Dus U2 + =— Dus U3 
or more compactly using suffix notation as 
Oz; 
dz: = du,, 


where the repeated suffix 7 on the r.h.s. implies summation from 1 to 3. Now 
the partial derivative 02 /Ou, means the rate of variation of z with u, while 
uz and ug are held fixed, so the vector 02 /Qu, lies in the uz and u3 coordinate 
surfaces and is therefore tangent to the u, coordinate curve. This enables a 
unit vector €ı to be defined in the direction of the u, curve, by 


€1 = Bur Êz Sh (6.1) 


where h, is a scale factor defined by 


(6.2) 


The unit vectors eg and eg are defined in a similar way, along with the scale 
factors hz and h3. 
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The displacement vector dz can then be written in terms of these unit 
vectors and scale factors as 


dx = h,e,du, + hzegduz + h3egdu3. 


Attention will be restricted to coordinate systems in which the unit vectors are 
orthogonal, so that 
€i: €j = 0i;. (6.3) 


Such coordinate systems are known as orthogonal curvilinear coordinates. This 
means that the coordinate curves are perpendicular to each other where they 
intersect. It will also be assumed that the coordinate system is right-handed, 
so that 

€1 X €2 = €63. (6.4) 


Locally, this coordinate system appears as a rectangular coordinate system 
with axes stretched by the factors h;, so that a change in u, of size du, leads to 
a change of distance hı du, in the e1 direction (Figure 6.2). Globally, however, 
the directions of the unit vectors e; vary in space, as do the scale factors h;. 


h, du, 


h,du, 
h, du, 


Fig. 6.2. An orthogonal curvilinear coordinate system appears locally as a rectangu- 
lar coordinate system. 


Having set up this general framework, formulae for various useful quantities 
can be derived in terms of the scale factors h; with reference to Figure 6.2. The 
length of a line element ds is found from 


ds? = dx - dz = hêdu? + hidu2 + Addu. (6.5) 


A surface element dS on the wu; coordinate surface generated by displacements 
duz, dug is rectangular and so has the area 


dS = hohg dus du3, (6.6) 
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and similarly for area elements on the uz and u3 coordinate surfaces. Finally, 
the volume element dV produced by displacements du, du2, du3 is again rect- 
angular and so its volume is 


dV = hı həh3 du, dus duz. (6.7) 


Example 6.1 


Show that the volume element in a right-handed orthogonal curvilinear coordi- 
nate system, h,h2h3 du, duz du3, can also be written in terms of the Jacobian 
of the transformation, J, defined to be the determinant of the matrix with 2,7 
element O2;/0u;. 

The determinant of a matrix can be interpreted as the scalar triple product 
of the vectors forming its rows or columns, so J can be written as the scalar 
triple product of three vectors: 


eA AO Lad 

Ou ðu» Oug 
= he, è hoe2 x h3€3 
= hıhzhz 


J 


since for a right-handed coordinate system, €2 xX €3 = e1 SO € : €2 X e3 = 1. 
Therefore the volume element can be written 


dV = J du duz duz. 


Example 6.2 


Parabolic coordinates (u,v, w) are related to Cartesian coordinates (z1, £2, £3) 
by the equations 


Ti = 2w, t2=u?-v*, 23=w. 
Sketch the u and v coordinate curves, find the scale factors hu, hy, hy and 
the unit vectors €u, €v, €w, and check that the (u,v, w) coordinate system is 
orthogonal. 

The u and v coordinate curves are the intersections of the v and u coordinate 
surfaces with the w coordinate surfaces. Since w = x3, the 21, x2 plane is a 
w coordinate surface. Consider the surface u = c, where c is a constant. In 
terms of the Cartesian coordinates this can be written 22 = c? — z? /4c? after 
eliminating v. Similarly the surface v = k can be written z, = —k? + r? J4k?. 
The u and v coordinate curves are therefore parabolas in the 2), z2 plane 
(Figure 6.3). 

The scale factors hu, hy, hy are just the magnitudes of the partial derivatives 
of the vector (1, 22,23) with respect to u, v and w: 
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X2 


x) 


Fig. 6.3. The coordinate curves of the parabolic coordinate system. 


hu = |(2v, 2u,0)| = 2Vu? + v?, 
hy = |(2u,—2v,0)| = 2V u? + v2, 
hw = |(0,0,1)|=1. 


The unit vectors are the vectors of the partial derivatives divided by the scale 
factors: 


eu = (v,u,0)/ vu? + v?, 
(u, —v,0)/ V u? + v2, 


(0,0,1). 


€v 


Cw 


Since the dot product of any two of these unit vectors is zero, the (u,v, w) 
system is orthogonal. 
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6.2 Grad, div and curl in orthogonal curvilinear 
coordinate systems 


In this section, formulae for the gradient of a scalar field and the divergence 
and curl of a vector field are derived for orthogonal curvilinear coordinate 
systems. In each case the physical definition is used, i.e. the definition which is 
independent of any coordinate system. 


6.2.1 Gradient 


The gradient Vf of a scalar field f is a vector perpendicular to the surfaces 
f = constant, defined by the equation 


df=Vf-dz 


where dz is an infinitesimal change in position and df is the corresponding 
change in f. Now from (6.1) we have 


dx = hie ,du; + hyeeduy + h3e3du3. 


If f is written as a function of the u; then 


_ Of Of Of 
Y = 3,0 + 5,7! © Bae duz 
10 ð ð 
= LBE pa + 2 BE miin + L B haie 


hy ðu," ha uz? hg Dug? 
Since this holds for any dz, the term in the large brackets is Vf: 


loaf. | 1 of 1 Of o, 
hi ðu, RA et i Ju” 


Z (Eile + 1 of 1 Of ) ae. 


Vf= (6.8) 
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6.2.2 Divergence 


The formula for the divergence of a vector field v = v,€1 + v2€2 + v3€3 in 
orthogonal curvilinear coordinates can be obtained using the definition (3.14) 


; 1 
V-v= im, yp? -nds. 


Since the coordinate system is orthogonal, the argument of Section 3.3 based 
on choosing ôV to be a small rectangular box can be repeated, with reference 
to Figure 6.2. The only difference in the argument is that the lengths of the 
sides of the box are scaled by the scale factors h;. On the u; surface on the 
right of the box, n = e1 So v : n = v, and 


v- -ndS = vi hz2h3 dugdu3 


where vı, ha and hg are evaluated at (u, + du;/2, ue, us). Similarly, on the 
opposite surface, 
v- ndS = —v hehe duzdu3 


evaluated at (u; — du; /2, u2, u3). Adding these two contributions and dividing 
by the volume h,h2h3du;du2du3 (exactly as in Section 3.3) gives the contribu- 
tion to V -v from these two surfaces as 


1 ð 
hı hə h3 Ou 


Note that h2 and hs cannot be cancelled out because in general they are func- 
tions of u,. The contributions to V -v from the other surfaces follow from cyclic 
permutation, so 


ae 
hihghs 


(vihzh3). 


ð ð ð 
V-v= (sn rhah) + Ju, 2h) + z (vshıha)) j (6.9) 


By combining the definitions of div (6.9) and grad (6.8), a formula for the 
Laplacian of a scalar field, V? f = V - (V f) can be obtained: 


TA A suL 
y f 7 hiheh3 E ( hı Ou 


EA ATEA (mer) 
Our he Our Ou3 hg u3 ' 


106 Vector Calculus 


6.2.3 Curl 


The curl of a vector field v in orthogonal curvilinear coordinates is found using 
the definition j 

n-Vxv = jim, gg f 9 4dr 
and following the argument of Section 3.4. To find the e3 component of V x v, 
consider a small rectangle in the u3 surface, with sides of length h,du, and 
hzduz (Figure 6.4). The line integral along the right-hand side of the rectangle 


€ 
(Ul), uu) À hidu, 


e; 
h, du, 


Fig. 6.4. Rectangle of four line segments for deriving the formula for V x v in 
orthogonal curvilinear coordinates. 


is approximately vzh2duz evaluated at (ui + duı/2,u2,u3) and the integral 
along the left side is approximately —vah2due evaluated at (u, — du; /2, ua, u3). 
Adding these two, taking the limit du, — 0, dug + 0 and dividing by the area 
gives a contribution 


l1 ð 
—— —— (ugh 
halo Ou ”? 2) 
to the eg component of V x v, and similarly the upper and lower sides of the 
rectangle generate a contribution 


1 ð 
Tate Bug 1) 


so the ea component of V x v is 


hy he Our 


The other components are found by permuting the indices. The determinant 
form of V x v is 


1 ð ð 
e3: V x v= —— (so orh) = Foz (hu) 


hi €i ha€2 h3€3 


a el! a 2 -2 
Vxv= hihzhz Ou uz ðus 


hiv, hove havz 
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EXERCISES 


6.1 Verify that for Cartesian coordinates the scale factors are all equal 
to 1. 
6.2 A coordinate system (u,v,w) is related to Cartesian coordinates 
(z1, 22,23) by 
Tı =uvw, T= uv(l -— w?) ?, 23 = (u? — v?)/2. 
(a) Find the scale factors hy, hy, hw. 
(b) Confirm that the (u, v, w) system is orthogonal. 
(c) Find the volume element in the (u,v, w) system. 


6.3 Find the scale factors and hence the volume element for the coordi- 
nate system (u,v,@) defined by 


Tı = ww cosh, zt2=uvsind, z3 = (u? — v?)/2, 


in which u and v are positive and 0 < 8 < 27. Hence find the volume 
of the region enclosed by the curved surfaces u = 1 and v = 1. 

6.4 Find the formula for Vf in a general orthogonal curvilinear coordi- 
nate system by writing V f in Cartesian coordinates and then finding 
the component of V f in the e, direction. 


6.3 Cylindrical polar coordinates 


Cylindrical polar coordinates (R,¢,z) are related to Cartesian coordinates 
(x1, 22,23) by 


zı = Reos¢, zo = Rsin@, r3 = Z. (6.10) 
The transformation in the reverse direction is 


R= 4z} +z},  ¢=tan™}(z2/%1) z= z3. (6.11) 


The coordinate system is shown in Figure 6.5. The scale factors are found 
using (6.2): 
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r 


Fig. 6.5. Cylindrical polar coordinates. 


hr = pr = |(cos ¢, sin ¢, 0)| = 1, 

h = |= p E Rsin ¢, Rcos ¢,0)| = R, 
Or 

h, —j= =], 
=| = 10.0.1) 


The unit vectors are found from (6.1): 
Ox 
eR = aR | be = (cos —, sin $, 0), 
Ox 
ep = ag/ he = (— sin ¢, cos ¢,0), 
Ox 
es = 5 [he = (0,0, 1). 


Note that the dot product of any two of the unit vectors is zero, so they obey 
the orthogonality condition. Also, since er x eg = ez, the coordinate system 
(R, ¢ġ, z) is right-handed. 

Using the scale factors, area elements on each of the coordinate surfaces can 
be found from (6.6). For example on the cylindrical R coordinate surface, 


dS = hgh, dġ dz = Rd¢dz. (6.12) 


Similarly, the volume element is, from (6.7), 
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dV = hrhgh, dRdġ dz = RdRdddz. (6.13) 


The formulae for the gradient, divergence, Laplacian and curl can be ob- 
tained by applying the results of Section 6.2: 


eo ees eh 
Vf= ARR T Rag tt pz E” (6.14) 
1 ug ae 
V-v= aap (Ron) + R Oo + (6.15) 
118 (pf), LEP, P 
Vl= RoR (Rap + Rg? D (6:16) 
_ 1 dv, Ovg Ovr v 
ane a E - 7) en & aa 
1 o Our 
R (Rvo oe Zer | €z. (6.17) 


Example 6.3 


Calculate the volume of a cone of radius a and height H. 

In cylindrical polar coordinates the equation of the cone is z = H R/a if the 
origin is chosen at the vertex of the cone. Then for a given value of z the range 
of Ris 0 < R < az/H. The volume is then 


H paz/H p2r 
J dV f J R dọ dR dz 
V 
az/ H 
N h 2r RdR dz 
ae 


Find the solution to Laplace’s equation in cylindrical polar coordinates that 
only depends on the distance R from the axis. 
If f obeys Laplace’s equation V? f = 0 and f depends only on R then from 


aie = na’ H/3. 


Example 6.4 


(6.16), 5 
ðf \ _ Of _ 
OR ( SE) =O ARS 
where c is a constant. Solving this differential equation gives 


f=clogR+d, 


where d is another arbitrary constant. 


110 l Vector Calculus 


6.4 Spherical polar coordinates 


Spherical polar coordinates (r,0,ġ) are related to Cartesian coordinates by 
Tı = rsin cos ¢, zə = rsin sin @, T3 = r cos ô, (6.18) 


and the inverse transformation is 


r=4/27+224+2?, 0= tan”! vi * 7 , ġ= tan"! a), (6.19) 
1 2 3 Ze i 


Note that ¢ here is equivalent to the angle ¢ for cylindrical polar coordinates. 
However, the reader should be aware that different authors use different nota- 
tion for the labelling of the angles in cylindrical and spherical polar coordinates. 
The variable r corresponds to its earlier use as the magnitude of the position 
vector T, representing the distance from the origin. 


x4 


A 


Fig. .6.6. Spherical polar coordinates. 


A sketch of the coordinate system is shown in Figure 6.6. Note that the 
ranges of the three coordinates are 


0O<r<o, 0<@O<7, 0< 6 < 2n. 
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A slightly modified form of spherical polar coordinates is used to measure 
position on the surface of the Earth. The longitude is the angle @ and the 
latitude is 7/2 — 0. 

The scale factors h,, hg and hg are 


h, = |(sin@cos ¢,sin@ sin ¢, cos@)| = 1, 
he = |(rcos@cos¢,rcos@sin ¢, —rsin6@)| =r, 
hg = |(-rsin@sin ¢,rsin@ cos ¢,0)| = r sinô. 


The unit vectors are 


e€, = (sin@cos@,sin@sin ¢, cos@), 
ee = (cos@cos¢,cos@sin ¢, —sin@), 
ep = (—sind¢,cos ¢,0) 


and it can be verified by taking the dot product of pairs of the unit vectors 
that the system is orthogonal. The area element on the spherical r coordinate 
surface is 


dS = hehg dô dọ = r° sin 0 dé do (6.20) 
and the volume element is 
dV = h,hehg dr dô dọ = r° sin 0 dr dô dọ. (6.21) 
The formulae for grad, div, Laplacian and curl are 
Of lof 1 OF o 
Vf = =e, ; ; 
f Or © T r 00 — r sinô O¢ re (6.22) 
V-v 3 Dy r)+ a 5 (sind vy) +— 60° (6.23) 


Of 1 ôf 1 æf 
2 277. ee ee 
Y= a ôr > (+ “4 HTT 0 aC mos) ý r? sin? 0 Og?’ on!) 


— _! (9 9 
Vxv = zaa (pnv) - aa 
TELET 
r \sin@ ð$ Or a i 
1/28 Ov, 
Les (Zw) = r) Cg. (6.25) 
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Example 6.5 


Suppose a sphere of radius a has a variable density p = po(1 — r/a) where po 
is a constant. What is the total mass of the sphere? 

The total mass M is just the volume integral of p. The limits on the integrals 
are0 <r<a,0<60< 7,0 < ¢ < 227 and the volume element is dV = 
r? sin @ dr dô dd. Choosing to do the integrals in the order r, 0, ¢, the mass is 


[ff oa —r/a)dV 


2n pr pa 
f f | po (1 — r/a)r? sin 6 dr dô do 
0 0 0 


2n pr 
f J po [r° /3 — rt /4a); sin 8 dô do 
o Jo 


M 


2n 
f poa? /12[ — cos 6] 5 do 
0 
T py a /3. 


Note that in triple integrals of this type, where the limits are constants, it is 
not strictly necessary to do the integrals one at a time. Thus it would have 
been possible to write after the second line, 


M = [r*/3 — rt /4a); [ — cos6]¢ (]o" = Tp a? /3. 


Example 6.6 


What proportion of the Earth’s surface lies further north than the 45° N lati- 
tude line? 

The required region of the surface is 0 < 0 < 1/4, 0 < ġ < 2m. Using the 
formula (6.20) for the area element, the area A is 


2r pr/ 
A= ff as = | i " »? sin 0d dé = 2ar? (1 — cos(1/4)) = mr? (2 — V2). 


As a proportion of the total area this is A/4ar? = (2 — //2)/4 = 0.15. 
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Summary of Chapter 6 


A curvilinear coordinate system (u4, u2, u3) is related to a Cartesian coordi- 
nate system (x1, 22,23) by z; = 2;(u1, u2, u3). 
e The unit vectors e; and the scale factors h; are defined for i = 1 by 


Ox 


ðu, 


hı = 


€1 = 1 
o] Ou ? 


and similarly for 7 = 2,3. 
e The system (u1, u2, u3) is orthogonal if e; : ej = 4;;. 
e The volume element in the (u1, u2, u3) system is dV = hı hoh du, duz duz. 
e Formulae for grad, div and curl in the (u1, u2,u3) system can be written 
down in terms of the scale factors and the unit vectors. 
e The two most important curvilinear coordinate systems are cylindrical polar 
coordinates, | 
zı = Rcosġ, z: = Rsin@, T3 =Z 


and spherical polar coordinates, 


zı = rsin ĝ cos, T2 = r sin ĝ sin @, £3 = r cos ô. 


EXERCISES 


6.5 A cylindrical apple corer of radius a cuts through a spherical apple 
of radius b. How much of the apple does it remove? 

6.6 Find the proportion of the Earth’s volume that is less than 30° away 
from the Equator. 

6.7 Find the divergence and curl of the unit vector eg in spherical polar 
coordinates. 

6.8 Find u: Vu for the vector u = eg in cylindrical polar coordinates. 

6.9 Find a formula for the R component of the Laplacian of a vector field, 
V?v, in cylindrical polar coordinates. Verify that the components of 
the Laplacian of v are not equal to the Laplacians of the components 
of v. 
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Cartesian Tensors 


7.1 Coordinate transformations 


At the very beginning of this book vectors and scalars were defined as ‘phys- 
ical quantities’. But what does this mean mathematically? In this chapter a 
precise mathematical statement is developed, using the idea that the physical 
quantity exists independently of any coordinate system that may be used. This 
new mathematical definition of vectors and scalars is generalised to define a 
wider class of objects known as tensors. Throughout this chapter attention is 
restricted to Cartesian coordinate systems. 

Consider a rotation of a two-dimensional Cartesian coordinate system z1, 
x2 through an angle 6 (Figure 7.1) to give a new coordinate system zr}, 25. 
Then by carrying out some simple geometrical constructions it can be seen 
that the coordinates of a point P in the 21, £2 system are related to those in 
the x}, £t} system by the equations 


A = 2),cosé4+272sin6, (7.1) 


Tə = 22c080- 72; sin, = (7.2) 
or in matrix form, 


zy \ cos sin6 Tı 
zh J) \ -sin@ cosé T2 J` 
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Fig. 7.1. Rotation of Cartesian coordinates through an angle 0. 


The 2 x 2 matrix relating (z},25) to (11,22) will be referred to as L: 


cos@ sinô 
L= , 7.3 
( —sin@ cosé ) (7:3) 
The matrix multiplication can be written in suffix notation, since 
zi = Lizi + Li2zr2 = L525, 
Ty = Laz + Lotz = L2;2;, 


where the repeated suffix j implies summation, so 
T; = Lij2;. (7.4) 


The rotation matrix L;; has one particularly important property. The inverse 
of the matrix is a rotation through —6, 


L-1- ( cos(—8) sin(—@) ) Z ( cos6é -sin ). 


—sin(—8@) cos(—6) sinô  cos@ 


which is the transpose of the matrix L. Thus LL? = I, or in suffix notation, 
L,;L4, = ôix. Since Lf, = Lgj, this can be written 


Lij Lez = Sik. (7.5) 


A matrix with this property, that its inverse is equal to its transpose, is said 
to be orthogonal. Using this property, the inverse of the transformation can be 
written down, simply by transposing the suffices: 


t= Libs: (7.6) 
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Another important property of the matrix L is that its determinant is 
|L| = cos? 6 + sin? 0 = 1. 


So far we have only considered a two-dimensional rotation of coordinates. 
Consider now a general three-dimensional rotation. For a position vector x = 
21€; + T2€2 + T3€3, the 2 component in the dashed frame is defined by 


T; = e; “c= €; °. €1T1 +e; * E212 +e; -€323 = e; -€jTj. 
This is of the form (7.4), where 
=al p. 
Lij = €e; ' ĉj, (7.7) 


so L;; is the cosine of the angle between e; and e;. By the same argument, the 
matrix which transforms from the dashed frame to the undashed frame has i,j 
element e; - ej = Lji, so again we see that the inverse of L is its transpose. 
Since LLT = I, the determinant of L obeys |L|? = 1, so |L| = +1. Orthogonal 


matrices with |L| = 1 represent rotations, while those with |L| = —1 are 
reflections. 
From (7.4) and (7.6), two further important properties of L follow: 
ðr; Oz; 
Oz; = Lij and ZA = Lji. (7.8) 


7.2 Vectors and scalars 


Now consider a vector v. Its components transform from one coordinate system 
to another in the same way as the coordinates of a point, so 


v; = Lijvj. (7.9) 


This equation gives the mathematical definition of a vector: v is a vector if 
its components transform according to the rule (7.9) under a rotation of the 
coordinate axes. 

Similarly, a scalar s is defined by the property that its value is unchanged 
by a rotation of coordinates, so 


s'=8. (7.10) 
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Using these new definitions of scalars and vectors, in terms of their trans- 
formation properties under a rotation of coordinate axes, a number of rigorous 
results can be proved, as illustrated in the following examples. Suffix notation 
and the summation convention are used throughout. 


Example 7.1 


Suppose that a and b are vectors. Show that their dot product a -b is a scalar. 
Since a and b are vectors, their components transform under rotation ac- 
cording to 
a; = Lijāj, b; = L;jbj. 
Now to show a-b is a scalar, we must show that its value in the dashed frame 
is the same as its value in the undashed frame. 


(a ` b)’ = a;b; = LijajLikbk = Lij Lina; by (7.11) 
= Ôjkajbk = akbk =a: b, (7.12) 


so a -b is a scalar. 


Example 7.2 


Suppose that f is a scalar field. Show that Vf is a vector. 
If f is a scalar then f = f'. To show that Vf is a vector we need to 
determine how it transforms under a rotation of coordinates. 


1 OF _ Of _ Of ds; 
(Vf); = ðr, Oz Ox; Oz! 
using the chain rule. Now making use of (7.8), 
of Of 
az’, = Lage Ox;’ 


so Vf obeys the transformation rule for a vector. 


Example 7.3 


A quantity is defined in a two-dimensional Cartesian coordinate system by 
u = (are, bz,)7. Show that this quantity can only be a vector if a = —b. 

If u is a vector, it must transform according to the rule u; = Liju; where 
Lij is the 2 x 2 rotation matrix (7.3). This gives 


i az, cos @ + br; sin 
— N —aze sin 8 + bz, cos 6 
but from the definition of u we also have 
»_ [f aza \ _ / -azı sinl + azz cose. 
br, / \ bricos +brzsin0 J` 
By comparing these two expressions we can see that they only agree if a = —b, 
so this is the condition for u to be a vector. 
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7.3 Tensors 


The definition of a vector as a quantity which transforms in a certain way 
under a rotation of coordinates can be extended to define a more general class 
of objects called tensors, which may have more than one free suffix. A quantity 
is a tensor if each of the free suffices transforms according to the rule (7.4). 
For example, consider a quantity 7;; that has two free suffices. This quantity 
is a tensor if its components in the dashed frame are related to those in the 
undashed frame by the equation 
Ti; = Lig LjmTem. (7.13) 
The rank or order of the tensor is the number of free suffices, so the quantity T;; 
obeying (7.13) is said to be a second-rank tensor. A tensor may have any number 
of free suffices. For example, a third-rank tensor P;jg transforms according to 
the rule 
Piję = LipLjqherPpgr- (7.14) 


The rule for a tensor of rank one is the same as the rule for a vector, so a vector 
can be regarded as tensor of rank one. Similarly, a scalar can be thought of as 
a tensor of rank zero. 

We have already met one second-rank tensor, 6;;, and a third-rank tensor, 
€ijz- Tensors can also be constructed from vectors, for example Ou,;/Oz; is a 
tensor. The demonstration that these quantities are indeed tensors is given in 
the following examples. 


Example 7.4 


Show that 6;; is a tensor. 

Consider the quantity LikLjmôkm. From the substitution property of 6,;, 
this is L;,L;,, which from the property (7.5) of L is ij. Now 6;, = 6;;, since 
6;; is defined the same way in any coordinate system. Thus 6;; obeys the tensor 
transformation law, 6: = Lit Ljmdkm- 


Example 7.5 


Show that €,;, is a tensor. 

Since ¢€,;, has three suffices, the appropriate transformation to consider is 
LipLjqLhkr€pqr- Using (4.10), this is €:;,|L| = €ijk, since |L| = 1 for a rotation. 
As for ij, €ijk is defined in the same way in all coordinate systems so E;j p= 
Eijk = LipLjqLkr€pgr- Therefore ¢,;, is a third-rank tensor. 
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Example 7.6 


If u is a vector, show that Ou;/0z,; is a second-rank tensor. 
Since u is a vector, u; = Li,ug. 


ðu; Our uk O21 Our 
ao = Ligo = Lie a = Lie biz: 
ðr; ðr; Oz Oz; Oz 


which is the transformation rule for a second-rank tensor. 


7.3.1 The quotient rule 


Tensors often appear as quantities relating two vectors, for example 
Qi = T,;;. (7.15) 


The quotient rule states that if (7.15) holds in all coordinate systems and for 
any vector b the resulting quantity a is a vector, then T;; is a tensor. 


Proof 


The quotient rule is proved as follows: Since a is a vector, 
a; = pupae = Lit T 55 b;. 


Since b is a vector, it obeys b; = Lmjbm (note that this is the inverse trans- 
formation, from the dashed to the undashed frame, so the suffices of L are 
transposed). Substituting for b; gives 


a; = LikTkjLmjbm. 
But since (7.15) holds in all coordinate systems, 
a = Timb: 


Subtracting these two results, 
(Tim g LikTkjLmj)bm = 0. 


If this result holds for any vector b, then the quantity in brackets must be zero, 


SO 
Tim = Lik LmjTkj. 


Therefore, T;; is a second-rank tensor. D 
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A more general form of the quotient rule also holds: if an mth rank tensor 
a is linearly related to an nth rank tensor b through a quantity T with m+n 
suffices, then T is a tensor of rank m+n. 


EXERCISES 


7.1 Show that the definition Li; = e; - e; is consistent with the matrix 
given in (7.3). l 

7.2 If u is a vector field, show that V - u is a scalar field. 

7.3 Given that a and b are vectors, show that the quantity a,b; is a 
second-rank tensor. 

7.4 Show that in a two-dimensional Cartesian coordinate system (z1, £2) 


the quantity 
2 
T7122 -T 


T —TıT? 
is a tensor. 
7.5. If @ is a scalar field, show that the quantity 
8? o 
Tix E Or; Oz; 


is a second-rank tensor. 

7.6 If Tj; is a tensor, show that T;; is a scalar. 

7.7 Write the divergence theorem in the form of suffix notation and 
hence obtain the analogue of the divergence theorem for a second- 


rank tensor T;j: 
If ea dV = fi] Tinj dS. (7.16) 
y Oz; s 


7.8 Write down the transformation rule for a tensor of rank four. 
7.9 If Qijkı is a tensor of rank four, show that Qijjı is a tensor of rank 
two. 
7.10 A quantity u; has the property that for any vector a, u;a; is a scalar. 
Show that the u; are the components of a vector. (This is a form of 
the quotient rule.) 
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7.3.2 Symmetric and anti-symmetric tensors 


A second-rank tensor 7;; is said to be symmetric if Ti; = Tj; and anti-symmetric 
if T;; = —T;;. A tensor of rank greater than two can be symmetric or anti- 
symmetric with respect to any pair of indices. For example ð;; is a symmetric 
tensor, while cij is anti-symmetric with respect to any two of its indices. 

It is important to verify that symmetry is a physical property of tensors, 
i.e. that if a tensor is symmetric in a Cartesian coordinate system it is also 
symmetric in other Cartesian coordinate systems. This can be confirmed as 
follows: suppose that A;; is a symmetric tensor, so A;; = Aji. Then in a 
rotated frame, 

Ay; = LigLjmAkm = LimLit Amk = Aji: 


so Aj, is also symmetric. 


Example 7.7 


Show that any second-rank tensor T;; can be written as the sum of a symmetric 
tensor and an anti-symmetric tensor. 

For any tensor T;;, the tensor S;; = Ti; + Tj; is symmetric. Similarly, 
Ai; = Ti; _ Li is anti-symmetric. Since Sij + Ai; = 2T;j, Tij can be written 
as Tij = ij/2 + A;; /2. 


Example 7.8 


The second-rank tensor T;j obeys cijkTjk = 0. Show that T;; is a symmetric 
tensor. 

By expanding out the implied double sum, for i = 1 we have €123T23 + 
€132132 = 0, which gives T23 = T32. Similarly the other required results follow 
from taking 2 = 2 and i = 3. 

The same result may be obtained more elegantly by multiplying the given 
equation ijk Tjk = 0 by Emni: 


0 Emni€ijkl jk 
(OmjOnk ~ OmkOnj) 15k 


Tmn TZ LANM) 


SO Tmn = tnm- 
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7.3.3 Isotropic tensors 


The two tensors 6;; and €:;, have a special property. Their components are 
the same in all coordinate systems. A tensor with this property is said to be 
isotropic. Isotropic tensors are of great importance physically, and it turns out 
that there are very few examples of isotropic tensors. This is illustrated by the 
following results. 


Theorem 7.1 


There are no non-trivial isotropic first-rank tensors. 


Proof 


Suppose that there exists an isotropic first-rank tensor (i.e. an isotropic vector), 
u = (u1, U2, u3). Now consider a rotation through 7/2 about the r3-axis, which 
is given by the matrix 


0 1 0 
L=| -1 0 0}. (7.17) 

0 0 1 
If u is a first-rank tensor then u; = Liju; = (u2, —u1, u3). Now if u is isotropic, 
u; = Ui, SO U1 = U2 and uz = —uı: Therefore u; = uz = 0. By considering a 
rotation about the z-axis in a similar way, it can be shown also that u3 = 0, 
so the only solution is u = (0,0,0). 0 


Theorem 7.2 


The most general isotropic second-rank tensor is a multiple of 6;;. 


Proof 


Suppose that a;i; is an isotropic second-rank tensor. Consider the rotation 
through 7/2 about the z3-axis given by (7.17). aij must obey a;; = LimLjn@mn, 
which in terms of matrix multiplication is a’ = LaLT. Carrying out these ma- 
trix multiplications gives the result 


Q22 —Q21 023 
LaL? = —ai2 Qil —Q@13 . (7.18) 
a32 —Q31 a33 


This must be equal to a;; if the tensor a;; is isotropic. The terms on the 
diagonal give a;; = a22. The other terms give a ;3 = a23 and a23 = —a3, from 
which a3 = a23 = 0. By considering the analogous rotations about the other 
coordinate axes it follows that a1; = a22 = a33 and that all the off-diagonal 
terms are zero, so a;; = A0;;, where À is an arbitrary constant. o 
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Theorem 7.3 


The most general isotropic third-rank tensor is a multiple of €;;,. 


Proof 


If aijk is an isotropic third-rank tensor, then 


Consider the same rotation (7.17), for which the only non-zero elements of L 


are Li = 1, Lo, = —1 and L33 = 1. Therefore for any choice of i, j and k 
in (7.19), only one term on the r.h.s. is non-zero. Choosing (i,j,k) = (1,1, 1) 
gives @i111 = a222 and the choice (i,j,k) = (2,2,2) gives a222 = —a111, SO 


a111 = @222 = 0. A different choice of rotation matrix would yield a333 = 0. 
By making further choices of (7,7,k) the following equations can be ob- 
tained: Q112 = —4@221, 0221 = 4112, 4122 = Q211, 0211 = —@122, Q121 = —4212; 
a212 = Q121. From these and the analogous equations involving the suffices 2 
and 3 it follows that all 18 elements with two suffices equal are zero. 
Finally, by considering the cases when i, j and k are all different, (7.19) 


gives a123 = —@213, 4231 = —Q132, 4312 = —Q32,. The analogous equations for 
rotations about the other axes can be used to show that a123 = a231 = a312 = 
—a321 = —Q132 = —4213, SO that ijk = A€ijx for some constant A. D 
Theorem 7.4 


The most general isotropic fourth-rank tensor is 
Qijki = ÀAijóki + HÔikôji + VÕuÔjk, (7.20) 
where A, p and v are constants. 


Proof 


An isotropic fourth-rank tensor must obey 
Gijkl = LipLjqherLisapgrs- (7.21) 


Using the rotation (7.17), only one of the 81 terms in the implied sum on the 
r.h.s. is non-zero. Since Liz = 1, Lz; = —1 and L33 = 1, a suffix 1 on the l.h.s. 
becomes a suffix 2 on the r.h.s., a suffix 2 on the |.h.s. becomes a suffix 1 on the 
r.h.s. and changes the sign, while a suffix 3 remains unchanged. By applying 


these rules, a3113 = @2223 = —@1113, SO @1113 = @2223 = 0. Similarly, any other 
term with three suffices equal and the fourth one different must be zero. Also 
@2113 = —Q1223 = —@2113 SO A2113 = 1223 = 0 and all similar terms with only 


one pair of equal suffices are zero. 
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The only remaining terms are those with two pairs of equal suffices and 
those with all four suffices equal. Applying the rotation (7.17) to terms in 
which the first two suffices are equal and the last two suffices are equal gives 
@1122 = 42211, 41133 = @2233 and a3322 = @331:1. Using the rotations about the 
other coordinate axes it follows that these six terms are all equal. Similarly, 
Q1212 = 42121 = 41313 = Q2323 = @3131 = 3232 and Q1221 = @2112 = @1331 = 
42332 = @3113 = 43223. The terms with all four suffices equal must obey a111 = 
Q2222 = 23333. Thus there can be at most four independent components of the 
tensor, @1122, @1212, @1221 and Q1111- 

To proceed it is necessary to consider a different rotation, for example the 
rotation through an arbitrary angle 0 about the z3-axis given by 


cos sin 0 
L= | —sinð cos 0 |. (7.22) 
0 0 1 


Using this rotation, a@1111 is related to all the terms with suffices equal to 1 or 
2. Applying (7.21) gives 


Q1111 = cos* 0 a111 + sint 0 a2222 
- 2 2 
+ sin” 8 cos” 6(@1122 + @2211 + @1212 + @2121 + @1221 + @2112). 


Simplifying this equation and using the relations above, the trigonometric fac- 
tors cancel out leaving 


@1111 = @1122 + @1212 + 41221, (7.23) 


so in fact there are only three independent components, which can be labelled 
Q1122 = A, 1212 = H, Q1221 = V. The tensor aiji can therefore be written in 
terms of À, p and v in the form (7.20). Note that this ensures that (7.23) is 
satisfied. 

o 


126 Vector Calculus 


7.4 Physical examples of tensors 


Tensors appear in many contexts, including fluid mechanics, solid mechanics 
and general relativity. Some of these applications will be described in Chapter 
8. The following two sections briefly consider two other examples of tensors. 


7.4.1 Ohm’s law 


Ohm’s law states that there is a linear relationship between the electric current 
j flowing through a material and the electric field E applied to the material. 
This can be written 

jJ =0E, (7.24) 


where the constant of proportionality ø is known as the conductivity (an inverse 
measure of electrical resistance). Note that (7.24) forces the vectors j and E to 
be parallel. For some materials, this may be true, but consider a substance with 
a layered structure made of different materials (Figure 7.2). For this material, 


Fig. 7.2. For a material made up of layers, the electric field E and the electric current 
j may not be parallel. 


current may flow more easily along the layers than across them. For example, 
if the substance is made of alternate layers of a conductor and an insulator, 
then current can only flow along the layers, regardless of the direction of the 
electric field. 

It is useful therefore to have an alternative to (7.24) in which j and E do 
not have to be parallel. This can be achieved by introducing the conductivity 
tensor, Cik, which relates 7 and E through the equation 
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Ji = On Ex. (7.25) 


Since 7 and E are vectors, it follows from the quotient rule that o,, is a tensor. 

The values of oip depend on the properties of the material. For example, 
suppose that there are alternating layers of a conductor with conductivity o 
and an insulator. If axes are chosen such that the z3 direction is perpendicular 
to the layers, then in this coordinate system 


do 0 0 
Ok = 0 oo 0 
0 00 


Now suppose that the material has no such layered structure, so that there 
is no preferred direction and is made of a uniform material with conductivity 
go. Such a material is said to be isotropic, meaning ‘the same in all directions’. 
In this case dik = do biz, SO 


Ji = Cn ER = O0 dik Ek = 00 Ei 


and so the simple rule 
3J = 09 E 


holds. This is why 4;, is said to be an isotropic tensor: it represents the relation- 
ship between two vectors that are always parallel, regardless of their direction. 


7.4.2 The inertia tensor 


Consider a body rotating with angular velocity §2. Then, as shown in Section 
1.3.1, the velocity vector at the position vector r is 


v=Q xr. 


The angular momentum of a particle of mass m is h = mr xv. The total angular 
momentum of a rotating body can then be determined as a volume integral, by 
considering dividing the body into small volume elements dV each with mass 
pdV, where p is the density of the body. The total angular momentum H is 


therefore given by 
J J il p(r xv), dV 
V 


IRG x (R x r)); dV 
J| oe. — (r-§2)r;) dV 


H; 
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I, p (776i; N; = Tj Nri) dV 
J p (r76;; — rirj) Rj dV. 
V 


Since 92; is a constant it can be taken out of the integral, leaving the equation 


A; = 1;;22;, (7.26) 


where I;; is called the inertia tensor of the body and is defined by 


nE Il p (r25,j — rar) aV. (7.27) 


Note that as in the previous example, the tensor appears as a quantity relating 
two vectors, and the quotient rule confirms that J;; is a tensor. The inertia 
tensor is an example of a symmetric tensor, since it is clear that I;; = Iji. 


Example 7.9 


Find the inertia tensor for a cube with sides of length 2a and constant density 
p, for rotations about its centre. 
To find J;; we need to compute two volume integrals. First, 


If pr? dV = f J J p (£? +y? + 2°) dz dy dz 
V —a J—a J —a 
3p | / J z? dz dy dz 


3p (2a) (2a) r? dr 


—a 


8pa° = Ma’, 


where M = 8pa? is the mass of the cube. The second volume integral is 


I, prir; dV. 


For i # j this is zero, since for example the integral of zy is zero since this is 
an odd function of z and y. For i = j, for example i = j = 1, we have 


J| o” dV = Ma’/3 


from the working of the first integral. Putting the two parts together, 
EEE 
Iij = Ma’ ôi; = Ma? ôi; /3 = z Ma ij. 


Note that the inertia tensor is isotropic. This means that for a cube rotating 
about its centre, the rotation vector and angular momentum vector are always 
parallel. 
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Summary of Chapter 7 


e Under a rotation of coordinate axes from a frame with unit vectors e; to a 
frame with unit vectors e;, the coordinates of a point are related by 


1 ee. ». o . 
Ti = Lijz; 


where Lij = eí - ej. 
e The inverse of the transformation is 


r ° — ee ! 
Li — LjiTj, 


so L7! = LT. Such a matrix is said to be orthogonal. In suffix notation, this 
result is written 
Lig Les = ôik. 


e A scalar s has the same value in each frame, s' = s. 

e A vector v transforms according to the rule v; = Lijvj. 

e If a quantity T;; transforms according to the rule Ti; = LikLimTkm then Tj; 
is a tensor of second rank. The rank of a tensor is the number of free suffices. 
Thus vectors are tensors of rank one and scalars are tensors of rank zero. 

e The quotient rule says that if a; = T,;5; and a is a vector for any choice of 
the vector b, then T;; is a tensor. 

e A tensor Tij is symmetric if T;; = 7}; and anti-symmetric if Tj; = —T}i. 

o 6;; and €;;, are tensors of a special type known as isotropic tensors. This 
means that their components do not change when the coordinate axes are 
rotated. A second-rank isotropic tensor must be a multiple of 6;; and a third- 
rank isotropic tensor must be a multiple of €;;,. 

e In physical systems, tensors frequently arise as quantities relating two vec- 
tors. This allows two vectors to be linearly related to each other without being 
parallel. Examples include the conductivity tensor and the inertia tensor. 
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EXERCISES 
7.11 B,, is an anti-symmetric tensor, so B,, = —B;,. Show that the 
anti-symmetry persists in a rotated frame, i.e. BY, = — B$,- 


7.12 If B,, is an anti-symmetric tensor, show that B,, = 0. 

7.13 The third-rank tensor A;;, is symmetric with respect to its first 
two suffices but anti-symmetric with respect to the second and third 
suffices. Show that all elements of A;;, must be zero. 

7.14 A quantity Aj; is related to a vector B by Aj; = €:j4.Bx. 

(a) Show that A;j is a tensor and describe its symmetry property. 
(b) Find an equation for B in terms of A;j. 

7.15 Find an isotropic fourth-rank tensor that can be written in terms of 
Eijk. 

7.16 Write down an isotropic fifth-rank tensor. Show that the most gen- 
eral isotropic fifth-rank tensor must have at least ten independent 
components. 

7.17 Show that the kinetic energy E of a body rotating with angular 
velocity §2 is related to its inertia tensor I;, by E = Ijk Rj R, /2. 
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Applications of Vector Calculus 


This chapter provides a brief introduction to some of the many applications 
of vector calculus to physics. Each of these is a vast topic in itself and is the 
subject of numerous books and a great deal of current research, so it is not 
possible to go into any detail in this book. However. a number of important 
governing equations and results can be obtained using the methods described in 
the previous chapters. In particular, it will be seen that the equations describing 
the behaviour of physical quantities such as electric fields and the velocity of a 
fluid are written in terms of the gradient, divergence and curl operators. 

The following sections discuss the flow of heat within a body, the behaviour 
of electric and magnetic fields. the mechanics of solids and the mecharics of 
fluids. There are however several other subjects which use the language of vector 
calculus, including the theories of quantum mechanics and general relativity. 
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8.1 Heat transfer 


In this section the equation describing the flow of heat within a solid body 
is derived. The argument is based on the law of conservation of energy, so is 
similar to the argument for the conservation of mass for a fluid given in Section 
5.1.1. 

Consider a solid with a temperature T which depends on space and time 
and a thermal conductivity K. Then the heat flows from hot to cold at a 
rate proportional to the temperature gradient, so the heat flux q is given by 
q = —KVWT. The minus sign appears here because the vector VT points in 
the direction of increasing temperature but the heat flows in the direction of 
decreasing temperature. 

Now consider an arbitrary region within the solid, denoted by a volume V 
with surface S and outward normal n. The thermal energy or heat content of 
a volume element dV is T cpdV where p is the density of the material and c is 
its specific heat. So the total heat content H of the volume V is 


H= [ff Teo: 


The rate of change of this heat content must equal the rate at which heat flows 
into the volume V, assuming that there are no sources of heat within V. This 
rate of inflow of heat is the integral of the heat flux —q-n over the surface 
S, where the minus sign appears since for heat to flow in, g must point in the 
opposite direction to n. Equating the rate of change of heat content with the 
rate at which heat flows into V gives 


Ma E cpd = f -q-nds= ff KYT- nas. 
S S 


The surface integral on the r.h.s. can be converted to a volume integral using 
the divergence theorem, giving 


M S LE ae If Vv -(KVT) dV (8.1) 


Finally, since the volume V is arbitrary, the volume integrals can be cancelled, 
giving 
OT 
cpa = V -(KVT), (8.2) 
since if (8.2) were not true at any point in space, then introducing a small 
volume V around this point would contradict (8.1). 
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If the thermal conductivity K is independent of position, then (8.2) can be 
simplified to give 


which is known as the heat seein or diffusion equation, and is of great im- 
portance since it occurs in many other contexts besides heat conduction. Other 
applications include processes of molecular diffusion such as the transport of 
chemicals within a living cell. The parameter k = K/cp is known as the ther- 
mal diffusivity of the material. Note that the units of the diffusivity k are 
length? /time. 

If the system is steady, so that there is no dependence on time, then T 
obeys Laplace’s equation, V?T = 0. 

The effect of the heat equation is to smooth out the distribution of temper- 
ature within a body. This is illustrated by the following examples. 


Example 8.1 
The surface S of a body is maintained at a constant temperature To. Show that 
the temperature T within the body approaches To as t —> oo. 

Define a new temperature 6 by 6 = T — To. Then since To is fixed, 6 obeys 


the heat equation 5 
a 2 
a kV“0 


within the volume V of the body and 0 = 0 on the boundary S. Multiplying 
through by @ gives 


1 062 
2 Ot 


Now integrating over the volume V and applying the divergence theorem gives 


2 
J, ae dV = Fig 9V6- ndS - k ff V0]? dV. 


The surface integral is zero since 6 = 0 on S, so 


2a 2 
zs M, 9? dV k fff [V0]? dV < 0, 


so the volume integral of 8? decreases unless @ is a constant. Since 0 = 0 on S 
the only possible value for this constant is zero, so the temperature decreases 
until 6 = 0, i.e. T = To everywhere. 


= k0 V?0 = kV - (0V0) — k|Y0}?. 
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Example 8.2 


Find a solution to the heat equation which is proportional to sin(az), where a 
is a constant. 

Seek a solution to (8.3) of the form T(z,t) = f(t) sin(az). Substituting into 
(8.3) gives 


df . _ EOR df _ 2 
Er sin(az) = —ka*fsin(at) => u ka* f. 
The solution to this differential equation for f is f = fo exp(—ka?t), where fo 
is a constant, so 
T (x,t) = fo exp(—ka7t) sin(az). 


Note that the amplitude of the solution decreases exponentially with time, and 
that shorter waves (larger a) decay more rapidly than longer waves (smaller a). 


8.2 Electromagnetism 


The fundamental equations describing the behaviour of an electric field E and 
a magnetic field B are written in terms of the divergence and curl of these 
vector fields: 


vE = £, (8.4) 
€0 
V.B = 0, (8.5) 
OB 
VxE = OR’ (8.6) 
; OE 
VxB = oj + Moto (8.7) 


where p is the density of electric charge, €o and uo are positive constants and 
J is the electric current density. These are known as Maxwell’s equations and 
from these many important properties of electric and magnetic fields can be 
derived. 

Each of the equations (8.4)—(8.7) can alternatively be written in an integral 
form. By integrating (8.4) over a volume V with surface S and applying the 
divergence theorem, we obtain 


ff e-nas = |)f Pav 
S v € 
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which is known as Gauss’s law. In words, this states that the total flux of 
electric field out of the surface S is proportional to the amount of electric 
charge contained inside the surface. Similarly, from (8.5) it follows that 


ff B-nds =o, 
S 


so that the net flux of magnetic field through any closed surface is zero. By 
integrating (8.6) over a surface S and using Stokes’s theorem, Faraday’s law of 
electromagnetic induction is obtained: 


ð 
f Ear =-2 ff B-nas 


Finally, from (8.7) we obtain 


$ B-dr =u || (5 +05) -ndS 
C s Ot 


which reduces to Ampeére’s law when there is no time dependence. 


8.2.1 Electrostatics 


Electrostatics is the study of steady electric fields. In a steady state, it follows 
from (8.6) that E is irrotational: 


VxE=-—0. 


Recall from Example 5.5 and Theorem 3.1 that this means that E is conserva- 
tive and that £ can be written as the gradient of a potential, E = -V &. The 
choice of the minus sign here is merely a convention. 

Now applying (8.4), the potential 6 must obey 


Vé=-V-E=-£. (8.8) 
€0 


This equation is known as Poisson’s equation. If the charge density p is zero, 
then # obeys Laplace’s equation V? = 0. The problem of determining the 
electric field and potential due to a stationary distribution of electric charges 
reduces to the problem of solving Poisson’s equation or Laplace’s equation in 
the appropriate geometry. These problems are generally referred to as ‘potential 
problems’ and the associated theory is known as ‘potential theory’. Some simple 
examples are given below. 
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Example 8.3 


A sphere of radius a contains a uniform distribution of electric charge p. De- 
termine the electric field E and the potential ® both inside and outside the 
sphere. 

Since the charge distribution is uniform within a sphere, it is most appropri- 
ate to use spherical polar coordinates and there is no preferred direction so both 
E and ¢@ only depend on the radial coordinate r. Referring back to (6.24) for 
the formula for the Laplacian in spherical polar coordinates, Poisson’s equation 
for the potential inside the sphere is 


12 „232 = —p/e 
r2 Or ar) 
Multiplying by r? and integrating gives 
> OF 
298 
" Or 


where C is a constant of integration, so 


OS 2 
Or pr [3€9 + C/r = E,, 


where Æ, is the component of E in the r direction. Now since the electric field 
is a physical, measurable quantity, the constant C must be zero to avoid a 
singularity at the centre of the sphere r = 0. Within the sphere therefore the 
electric field and potential are: 


= —pr? /3e9 + C, 


E, = pr /[3€0, $ = —pr* /6e9 + D, 


where D is an arbitrary constant which always appears in potentials. 
Outside the sphere, Laplace’s equation holds, 


19 (e2) = 0, 


f 


r? Or Or 
SO 
>» Of ð ' 
2 — = — = 2 = = 
"3, 7 Cir E, 


In this case the constant C cannot be set to zero since r = 0 is outside the 
region of consideration. Instead, C can be found by imposing that the electric 
field is continuous across the surface r = a. Equating the values of E, at r = a 
from the formulae for E, inside and outside the sphere gives 


pa/3eg = —C/a® => C = -—pa?/360. 


Hence, outside the sphere 
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= 3cor? j 3€9r’ 
where the constant of integration associated with # has been set to zero. Note 
that the electric field obeys the inverse square law: E, is proportional to r~?. 


Example 8.4 
Find the electric field Æ and the potential @ due to a point charge Q. 

This can be done using the result of the previous example, by considering 
a point charge as a small sphere of uniform charge density. The total electric 
charge Q contained within a sphere of radius a with uniform charge density p 
is just the density multiplied by the volume, Q = 47a°p/3. The electric field 
and potential can then be written in terms of Q as 


Q Q 


Anegr2’ Aner 


rT 


8.2.2 Electromagnetic waves in a vacuum 


In a vacuum, where there is no electric charge and no electric current, Maxwell’s 
equations take a simple and symmetric form: 


V-E = 0, (8.9) 
V-B = 0, (8.10) 
0B 
VxE = TE’ (8.11) 
OE 
VxB = HotoAE (8.12) 
Taking the curl of (8.11) and using (4.23) gives 
Vx(Vx BE)=V(V-E)-V E=-Vx -r 
Now using (8.9) and (8.12) this can be written 
; ðE 
-V? FE = -— Z 
V ð (uoco ot ) ’ 
so E obeys : 
ð E 
aa eV E (8.13) 


where ¢? = 1/po€p. This equation is called the wave equation since its solutions 
are waves travelling at speed c, as shown in the examples below. The constant 
c is the speed of light and the waves are known as electromagnetic waves. 
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Radio waves, light and X-rays are all examples of electromagnetic waves. These 
different types of waves have different frequencies but they all travel at the same 
speed c. 


Example 8.5 


Suppose that E = (f(z,t),0,0) in a Cartesian coordinate system so that the 
wave equation (8.13) becomes 


Of _ .ef 
ð? Oa?" 


Show that f(z,t) = sin k(x — ct) is a solution for any value of the constant k 
and interpret this solution physically. 
For the function f(z,¢) = sin k(x — ct), 
O° f 2: O° f 
2 = —(kc)* sin k(x — ct) and an? = 


so f(z,t) = sin k(x — ct) obeys (8.14). Physically, this solution corresponds to 


(8.14) 


—k* sin k(x — ct), 


f(x,0) fixt) 


Fig. 8.1. The solution f(z,t) = sin k(x — ct) of the wave equation, at t = 0 (solid 
line) and at some later time (dashed line). 


a sine wave travelling to the right at speed c. At t = 0, f(z,0) = sin kz which 
is a sine wave which has f = 0 at z = 0. At a later time t, the point at which 
f = 0 has moved to the position z = ct, so the wave has moved to the right 
a distance ct (Figure 8.1). The speed of the wave is therefore c. Similarly, it 
can be shown that the function g(z,t) = sin k(x + ct) is a solution of (8.14), 
representing a sine wave travelling to the left at speed c. 
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Example 8.6 


Show that E = Eof(k-x — wt) is a solution to the wave equation (8.13), where 
Eo and k are constant vectors, f is any function and w = +c|k|. 

Let u = k -x — wt so that the solution under consideration is E = Eof (wu). 
Then 


ðE ,, df Ou _ af 
Ot Eo du Ot wes du 
and similarly > P P 
OE f Ou f 
Of omoa ae 


To find VE, note that since Eo is constant, V°? (Eof(u)) = EoV?f(u) = 
EoV -(Vf(u)). The gradient of f(u) is 

_ df _ df 

Vf(u) = q Ye z — wt) = que 


and taking the divergence of this gives 


vsa = 0 (F). pe yuks TS py 


d du? 
Therefore, the wave equation is on provided that 
w? Eo ai = Bot Liki > w= c*|k|’. 


So the function f is arbitrary and the only condition is that the frequency of 
the wave w must be related to the speed c and the constant vector k (which is 
known as the wave vector) by w = +clk|. 


EXERCISES 


8.1 Use dimensional analysis to determine how the time taken for heat 
to diffuse through a body depends on the size L of the body and its 
thermal diffusivity k. Hence answer the following questions. 

(a) If it takes six hours to defrost a frozen chicken, how long. would 
it take to defrost. a woolly mammoth? 

(b) Cookery books state that the time taken to cook meat is, for 
example, twenty minutes per pound plus twenty minutes. Is this 
correct? 

8.2 Show directly from Maxwell’s equations that the charge density p 
and the electric current density 7 obey the conservation law 
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8.3 Derive the formula for the electric field Æ due to a point charge Q 
using Gauss’s law. 

8.4 Show from Maxwell’s equations in a vacuum that the magnetic field 
B obeys the wave equation. 

8.5 For the electromagnetic wave in which the electric field is given by 
E = Eo f(k - x — wt), calculate the corresponding magnetic field B. 
What can be deduced about the directions of the vectors E and B? 

8.6 Using Maxwell’s equations in a vacuum, obtain an equation in the 
form of a conservation law for the rate of change of the energy w = 
|B|? /2 + |E|? /2c? of an electromagnetic wave. 


8.3 Continuum mechanics and the stress tensor 


Continuum mechanics is the study of continuous media including solids, liquids 
and gases. Solids have the property that when acted on by a force they deform 
but then reach an equilibrium in which the internal forces within the material 
balance the imposed force; this behaviour is described in Section 8.4. Fluids, 
which include liquids and gases, move continuously when subjected to a force 
and so require a different treatment (Section 8.5). Common to both solids and 
fluids is the concept of stress and the stress tensor described below. 

Consider a small section of a surface, with area 6S and unit normal n, 
within a material (Figure 8.2). The two sides of the material are labelled side 1 
and side 2, with n pointing in the direction of side 2. The material on side 2 
exerts a force 6F on the material on side 1 through the surface 6S, due to 
interactions between the molecules of the material. Since the force is exerted 
through ôS, the magnitude of the force is proportional to the area 6S. The 
force is a vector quantity, and the only available vector is n, so dF must be 
some quantity multiplied by n. However, 6F need not be parallel to n. The 
angle between ôF and n can be allowed to vary by using a tensor to relate 6F 
and n, 

OF; = Pynj 6S, (8.15) 


where P;; is known as the stress tensor of the material. Since 6F and n are 
vectors, it follows from the quotient rule that P;; is a tensor. 

Consider now the force exerted by side 1 on side 2. This is obtained by 
reversing the direction of the normal and using (8.15), so the force is —P,;n; 6S. 
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OF 


Fig. 8.2. A force óF is exerted on the material on side 1 through the surface 6S. 


This is just —dF;, so this is consistent with Newton’s third law as the forces 
are equal and opposite. 

Now consider a region V within the material, with surface S and outward 
normal n (Figure 8.3). The total force exerted on this region by the surrounding 


Fig. 8.3. The total force on a volume V is found by integrating over the surface S. 


material is found by integrating the force due to all the surface elements over 


the surface S: ap 
F= ff Pynjas = |j] 3i dV 
ge v 02%; 


using the divergence theorem for tensors (7.16). In general there may be other 
forces acting on the volume V, for example the force due to gravity; such long- 
range forces are known as body forces. Suppose, however, that body forces are 
negligible so that the only forces acting on the volume V are those exerted 
through the surface S. If the material is in equilibrium, then the total force 
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acting on the region must be zero. Since this must be true for any arbitrary 
volume V, the stress tensor P;; for a system in equilibrium must obey 

OP. 

t="), 

Oz; 

An additional constraint on P;; comes from the consideration of the moment 

of the forces acting on a volume V. The moment of a force F about a point O 
is r x F where r is the position vector of the point of application of the force 
relative to O. The total moment of the forces acting through the surface S is 


therefore 
fi] EijkTj dF, = fi] EijkLjPkmnNm dS. 
S S 


Using the tensor form of the divergence theorem again, this can be written as 


= 9Pem 
I Eijk Jg, (24P em) dV = If Eijk (5m Pim + 2; a) dV. 


Now for a body in equilibrium we have already shown that the second term in 
the bracket is zero. The total moment of the forces must also be zero (otherwise 
the body would start to rotate), and since the volume V is arbitrary it follows 
that 


Eijk Pej = 0. 


This means that the tensor Pk; is symmetric, as shown in Example 7.8. There- 
fore the stress tensor for a material in equilibrium must obey the two constraints 


i = and Py = Pj. (8.16) 


Example 8.7 
Suppose that a material is subjected to a body force b per unit volume. Find 
the total force f acting per unit volume. 

For a volume V the total force is the volume integral of f over V. This is 


made up of the volume integral of the body force b plus the surface integral of 
the forces due to the stress tensor: 


[ffs = [Jf av ff Panas = fff n+ Bear 


Since this is true for any volume V, 


OP; 


fi =b) + On; 


(8.17) 
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8.4 Solid mechanics 


In the theory of the mechanics of a solid material it is generally assumed that 
there is a relationship between the stress tensor of the material and the strain 
of the material. The strain is the amount of stretching or deformation of the 
material. Suppose that, as a result of imposed forces, a material is deformed 
so that the material originally at position vector r is moved a distance u(r). 
Then the strain tensor Ei; is defined by 


_il ðv; OU; 
Ei; — 9 (Z + Sa) - (8.18) 


The simplest assumption is that there is a linear relation between the stress 
tensor and the strain tensor. This is known as Hooke’s law. Since the stress 
tensor and the strain tensor are both second-rank tensors, the quantity relating 
them will in general be a fourth-rank tensor. If a material obeys the linear 
relationship 

Pij = cij Eri (8.19) 


then it is said to be an ideal elastic solid. 
If the further assumption is made that the material is isotropic, then the 
tensor C;;4: is isotropic, so from Section 7.3.3, 


Pi; = (bij der + bin jt + Vôuðjr) Eri 
= NOiz Ena + pEiy + vEji 
Abi; Ekk +(u+ v)Eij (8.20) 


since from its definition E;; is symmetric. There are therefore only two inde- 
pendent constants in the relationship between stress and strain for an isotropic 
elastic solid. Since v is arbitrary, we may take v = p to obtain 


Piz = Oij Ekk + 2pEij. (8.21) 


The two constants À and p are known as Lamé’s constants. The inverse rela- 
tionship, giving the strain in terms of the stress, can be found by setting i = j, 
giving 

Py, = 3A Lg, + 2pk = (3A + 2p) Exr 
and then substituting Ex, = Px, /(3A + 2u) in (8.21), giving 


1 Nbi5 Pkk 


) EEE Pee. 
7 2p ” (2p) (3A + 2p) 
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Example 8.8 


Find the strain tensor and the stress tensor for an isotropic elastic material 
when it is subjected to 

(a) a stretching deformation v = (0,0, az3); 

(b) a shearing deformation v = (bz3,0,0). 

The two deformations are illustrated in Figure 8.4. 


(a) (b) 


X3 


-æ = æ æ a q 


Xi 


Fig. 8.4. Deformation of a cube of material (dashed line) by (a) v = (0,0,az3), (b) 
v = (bz3, 0,0). 


(a) Here the only non-zero element of E;; is E33 = a. Hence Ex, = a and 
so the non-zero elements of P,; are Py; = Pz2 = Aa and P33 = (A + 2y)a. 

(b) In this case E3; = Eis = 6/2 and other elements of E;; are zero. Since 
Ekk = 0 the only non-zero elements of Pj; are P3, = Py3 = pb. 

These two simple examples show that the diagonal components of P;; rep- 
resent forces of stretching or compression of the material, while the off-diagonal 
components represent shearing forces. 


Example 8.9 


Show that the anti-symmetric tensor 


Git (2 2 se | 
we aD Oz; Oz; 


represents a rotation of the material. 
Since S;; is anti-symmetric, $1; = S22 = S33 = 0 and S;; has only three 
independent components, S12, S23 and S31. The first of these is 


1/o ð 1 1 
Sa=3 (52-52 = -3 [V x v]; = -5612 [V x v], - 


Similar equations hold for the other components, so S;; is related to V x v by 
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1 
Sij = Tgĉijk [V x v]; x (8.22) 


Thus the components of S;; depend only on V x v, which from Section 3.4.2 
represents a rotation of the body. 


Example 8.10 


Show that in an isotropic material P;; cannot depend on the anti-symmetric 
tensor Sij. 

Following the arguments above, if Pj; = CijkiSkı and the material is 
isotropic, then the analogue of (8.20) is 


Pij = Abig Skk + (u = v)Sij. 


Since S;; is anti-symmetric, Sk = 0. This means that the symmetric tensor 
P,; is equal to a constant multiplied by the anti-symmetric tensor Sij, which 
is clearly a contradiction. 


8.5 Fluid mechanics 


An important property which distinguishes a fluid (such as a liquid or a gas) 
from a solid is that a fluid is unable to support a shear stress. This means 
that if a fluid is subject to a shear stress, such as in Example 8.8(b), the fluid 
will move continuously. For a fluid at rest, therefore, the stress tensor can only 
include diagonal elements. Furthermore, if the fluid is subjected to a stretching 
force, as in Example 8.8(a), the fluid will stretch and continue to move for as 
long as the force acts. This motion can only be avoided if the forces acting on 
the fluid are the same in all directions. Hence the stress tensor in a fluid at rest 
must be isotropic, 

Pij = —pdi;, (8.23) 


where p is the hydrostatic pressure of the fluid, which in general is a function 
of position. The pressure can then be defined in terms of the stress tensor by 


p= —P,;/3. (8.24) 


Now suppose that the fluid is in motion, with a velocity given by the vector 
field u(r). In this case, there is an additional contribution to the stress tensor 
due to the motion, so 

P= —pi; + dij (8.25) 
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where d;; = 0 when the fluid is at rest. Since (8.25) does not uniquely define p 
and dij, it can be imposed that (8.24) holds for a fluid in motion, which means 
that d;; = 0), 

The tensor d;; represents the forces due to the fluid motion. These forces 
arise due to friction between adjacent molecules of the fluid that are moving 
at different velocities. It is generally assumed that there is a linear relationship 
between d;; and the tensor representing the difference in velocities, Ou;/Oz;. A 
fluid behaving in this way is said to be Newtonian, and experiments show that 
this is a good approximation for most fluids. The analysis of the relationship 
between d;; and Ou;/0z,; now follows exactly as for the case of a solid (Section 
8.4). If it is assumed that the fluid is isotropic, then d;; can only depend on 
the symmetric part of Ou,/Oz;, 


peso Conga 
7" 2\@2; Ox; )” 


which is known as the rate-of-strain tensor. The argument of Example 8.10 
shows that d;; cannot depend on the anti-symmetric part of Ou,/Oz;. The 
form of the dependence of d;; on e;; is 


dij = NOij kk + 2Mei;, (8.26) 
where, as in (8.21), À and yp are constants. The constraint dj; = 0 means that 


3A + 2p = 0. Eliminating A, the form of the stress tensor for a Newtonian fluid 
is then 


2 

Pij = —pois + 2peij — z HôijEkk. (8.27) 

The constant p is the viscosity of the fluid, which represents the fluid’s ‘sticki- 
ness’ or friction. : 


8.5.1 Equation of motion for a fluid 


Having obtained the formula for the stress tensor in a fluid, it is now possible 
to write down the equation of motion. The equation of motion is derived from 
Newton’s second law, force = mass x acceleration. The force can be obtained 
from the stress tensor using (8.17). The force f per unit volume due to body 
forces b and the stress tensor P;; is 


= OP ij 
fi = bi+ Oz; 
Op Ou; 3u; 2 Ou, 
beste Saba ae pee) j 
Ox; PH (a 4 ob) 3” Ox,0z; 
Op Ou; 1 Ou; 
bp = 2 i SRI Ii a | pine eam 
Oz; TH 5s On, j 3 ðr;ðz; 
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In vector notation this can be written 


1 
f=b-Vp+uV’u+ guVV u. 

To obtain the acceleration of a fluid particle, consider the change in veloc- 
ity of a particle in a short time interval ôt. A particle at position r moving 
with velocity u(r,t) moves a distance u(r, t)ôt in this time. Therefore its new 
velocity is u(r + u(r,t)ót,t + ôt). The velocity is a function of four variables, 
u = u(z,y,2,t), so expanding this new velocity using Taylor’s theorem gives 

Ou Ou Ou ðu 


u(r, t) + Uz dt a + uy dt + usdt— + it 


a] 
= u(r,t)+dtu-Vut it. 
The acceleration a of the particle is the limit of the difference between the new 
and old velocities divided by ôt, so 


P haere Prt 


Ot 


Since the mass per unit volume of the fluid is just the density p, the equation 
of motion is 


u(r + u(r,t)dt,t + dt) 


p (F +u-Vu) =b- Vp+uV’ u+ SUWY u. (8.28) 
This important equation is known as the Navier-Stokes equation. 

From now on, assume that the fluid is of constant density, so that p does 
not depend on time or space. Then, from Section 5.1.1, where the equation for 
conservation of mass for a fluid was derived, the velocity field u obeys V -u = 0. 
There are two more important equations of fluid mechanics that can be derived 
from the Navier-Stokes equation; these are described in the following sections. 


8.5.2 The vorticity equation 


Assume that the body force b is conservative, so that it may be written as the 
gradient of a potential, b = -V $. By using the vector identity 


u:Vu=V(\|ul?/2)-uxVxu 


which was derived in Example 4.14, the Navier-Stokes equation can be written 


p (F + V(|ul?/2) — u x w) = -Vp — VĒ + pV*u, (8.29) 


148 Vector Calculus 


where w = V x u is the vorticity of the fluid, which as shown in Section 3.4.2 
is proportional to the local rate of rotation of the fluid. Now by taking the curl 
of (8.29), the equation can be simplified considerably because the three terms 
involving grad will then disappear. The pressure and the body forces are then 
eliminated and an equation for the rate of change of the vorticity is obtained: 


= -Vx(uxw)= Vw, (8.30) 
where we have made use of the fact that the order of the curl and Laplacian op- 
erators may be interchanged. By expanding the curl of the cross product using 
(4.30) and noting that both u and w are solenoidal, (8.30) can be alternatively 
written 3 

w 


E a Vosa Vus V 
ry +u-Vw-w- Vu ra w. (8.31) 


This is known as the vorticity equation. 


Example 8.11 


Obtain the simplified form of the vorticity equation for two-dimensional flow, 
u(r) = (u(x, y), v(z, y),0), of an incompressible fluid. 
If u = (u(x, y), v(z, y),0), then the vorticity w is 


ðv Ou 
w=Vxu= (0,0, = am By? J (0, 0, w(z, y)). 


Thus the vorticity vector has only one component, which is perpendicular to the 
plane of motion of the fluid. This means that the term w - Vu in the vorticity 
equation is zero, since w : Vu = wOu/0z = 0. The vorticity equation simplifies 
to the scalar equation 


Ow _ Hp? 
ry +u-Vw= ~ W. (8.32) 


Example 8.12 


An incompressible fluid flows along a straight two-dimensional channel and 
the velocity u is parallel to the walls of the channel. Show that the vorticity 
equation reduces to the diffusion equation. 

Choosing the z-axis to be parallel to the channel walls, the velocity u has 
the form u = (u, 0,0). As the fluid is incompressible, V -u = 0 so 0u/Oz = 0. 
Since u is independent of z, so is w, so the second term in (8.32) is u- Vw = 
uðw/ðx = 0. The vorticity equation is then 


Ow B2 


which is the same as the diffusion equation (8.3). 
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Ow 
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ee Voca Nua y 
a tu Vw-w-Vu a w. (8.31) 
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ðv ðu 
w=Vxu= (0,0, 53 = By? a (0, 0, w(z, y)). 


-Vx(uxw)= TVu, (8.30) 
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equation is zero, since w- Vu = wðu/ðz = 0. The vorticity equation simplifies 
to the scalar equation 


Ow _H~2 
Y +u-Vw= ~ w. (8.32) 


Example 8.12 


An incompressible fluid flows along a straight two-dimensional channel and 
the velocity u is parallel to the walls of the channel. Show that the vorticity 
equation reduces to the diffusion equation. 

Choosing the z-axis to be parallel to the channel walls, the velocity u has 
the form u = (u, 0,0). As the fluid is incompressible, V -u = 0 so 0u/Oz = 0. 
Since u is independent of z, so is w, so the second term in (8.32) is u- Vw = 
uðw/ðx = 0. The vorticity equation is then 

Ow _ u 


_ Hoe 
ares 


which is the same as the diffusion equation (8.3). 
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8.5.3 Bernoulli’s equation 


A second important equation can be obtained by taking the dot product of u 
with (8.29). If the flow is steady, so the time derivative is zero, and the viscosity 
of the fluid is negligible, then 


u-V(plul?/2+p+9) =0. (8.33) 


This means that the quantity p|u|*/2 +p + does not change in the direction 
of u, so that it is constant along the path of a fluid particle. If in addition, the 
vorticity is zero, then (8.29) becomes 


V(plul?/2+ p+) =0 (8.34) 


so p|u|*/2+p+@ is constant everywhere. These two results are different forms 
of Bernoulli’s equation. In words, Bernoulli’s equation means that where the 
velocity of a fluid is high, its pressure is low. This important result can be used 
to explain a wide range of physical observations, including the swerving of a 
spinning ball and the lift generated by an aeroplane’s wing. 


Example 8.13 


Why does a spinning ball swerve? 

Consider a ball travelling to the left with velocity u. Equivalently, the ball 
may be considered to be stationary, with air flowing past it with velocity v = 
—u. Now suppose that the ball is spinning anticlockwise, as shown in Figure 
8.5. As the ball spins, it drags air with it, rotating around the ball. The two 
components of the fluid motion due to the translation and the rotation are 
oppositely directed above the ball but point in the same direction below the 
ball. Therefore the total speed of the fluid is lower above the ball than below 
the ball. According to Bernoulli’s equation, the pressure is greater above the 
ball than below. This pressure difference exerts a force on the ball which is 
directed downwards in Figure 8.5, causing the ball to swerve. 


Example 8.14 


To model the effect described qualitatively in the previous example, consider a 
cylinder of radius b and assume that in the frame of the cylinder the velocity 
at the surface of the cylinder is (v sin @ — a)ég in polar coordinates. Assuming 
that the viscosity is negligible and that Bernoulli’s equation is valid, find the 
pressure at the surface and hence calculate the total force acting on the cylinder. 

If Bernoulli’s equation holds and there are no body forces, then p = c — 
p|u|?/2 where c is constant, so 


p = c — plv? sin? 0 — 2avsin 8 + a”)/2. 
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Fig. 8.5. A ball moving to the left and spinning anticlockwise experiences a net 
downward force. 


Now if there is no viscosity the stress tensor is simply P;; = —p ði; so the total 


force is 
r= f PynjdS = f} —pn, a8. 
S S 


In a Cartesian coordinate system (z1, 22), n = (cos@,sin@) and the area ele- 
ment dS is bd@ per unit length of the cylinder. The force per unit length in the 
z> direction is 


2r 
Fy = f —(c — plv? sin? 6 — 2av sin 0 + a*)/2) sin 0 b d8 
0 


27 
f — pav sin? 6 b dO 
0 
= -—m7paub. 


The force in the zı direction is zero since the functions sin? 0 cos @, sin @cos0 
and cos@ all integrate to zero. This means that there is no drag force on the 
cylinder, which shows that this simple model is inadequate in some way. In fact 
it is the neglect of viscosity which is not valid. Even though the viscosity of a 
fluid such as air is small, it can have a large effect on the solution. 
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Summary of Chapter 8 


Many important equations of mathematical physics are written in terms of 
grad, div and curl. 
The heat equation, or diffusion equation, is 


OT ; 
= AVL. 


This describes heat transfer within a body and other physical processes in- 
volving molecular diffusion. 
Electric and magnetic fields obey Maxwell’s equations (8.4)—(8.7). In a steady 
state the electric field is irrotational and its potential & obeys Poisson’s 
equation V? = —p/éo where p is the charge density and €o is a constant. 
In a vacuum, Maxwell’s equations lead to the wave equation which describes 
the motion of electromagnetic waves. 
In a continuous material such as a solid, a liquid or a gas, molecular forces 
are transferred through a surface ôS with normal n according to the formula 
ôF; = Pijn; 6S, where P;; is the stress tensor of the material. If the material 
is in equilibrium, P;; must obey 
ðP.. 
Ox; =0 and Pj = Pii 
In an isotropic, ideal elastic solid with a deformation v(r) the strain tensor 
E;;j and stress tensor P;; are given by 
1 f Ou; | dv; 
Beet (= etd) 
2 \ðzr; Oz; 
where À and u are constants. 
In a Newtonian fluid moving with velocity u(r), the rate-of-strain tensor e,; 
and stress tensor P;j are 
pe 1 Ou; Ou; 
Y 2\ 02; Oa; 
where p is the pressure in the fluid and p is its viscosity. 
If the fluid is incompressible then V -u = 0 and its equation of motion is the 
Navier-Stokes equation 


) j Pi; = Aĝ;ij Ekk + 24Eij, 


2 


) ’ P;; = —p oi; + 2 ues; — 3 


Uòi; Ekk, 


p(k +u- Vu) =b- Vp+ uV’u, 


where p is the density and b is the body force acting. 
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8.7 A material has a stress tensor P,; that takes the form 
0 a T2 
P;; = 0 zy  bz2T3 


T CT2T3 T3 


in a Cartesian coordinate system (21, £2, £3). 

(a) If there are no body forces acting and the material is in equilib- 
rium, find the values of the constants a, b, c, m and n. 

(b) For these values of the constants, find the magnitude and direc- 
tion of the force exerted on the surface zı = 1, 0 < z2,23 < 1, by 
the material in the region zı > 1. 


8.8 Verify (8.22) by expanding the r.h.s. using suffix notation. 
8.9 Show that for an isotropic elastic solid in equilibrium, the deforma- 


tion v must obey (A+ u)V V -v + pV720 = 0. 


8.10 An isotropic elastic solid with Lamé constants À and yp is subjected 


8.11 


to a deformation v, = ax, 22, v2 = b(z? — z2), v3 = 0. 
(a) Find the strain tensor E;;. 
(b) Find the stress tensor P;;. 
(c) Determine whether it is possible for the material to be in equi- 
librium. 
A compressible fluid with negligible viscosity is initially at rest with 
uniform density po and pressure pọ, with no body forces. A small 
perturbation is then introduced: so that there is a velocity u(r, t) 
and the density becomes po + pı (r,t). 
(a) Assuming that products of the small quantities u and pı can 
be neglected, show that the equation for conservation of mass (5.9) 
becomes 

Opi 


Ot + poV-u=0. (8.35) 


(b) Assuming that the perturbation pı to the pressure is related to 
pi by pı = ap, where a is constant, show that the Navier-Stokes 
equation reduces to 


Ou 
pose = —aVp. (8.36) 


(c) Hence show that the density perturbation pı obeys the wave 
equation and interpret this result physically. 


Solutions 


Solutions to Exercises for Chapter 1 


1.1 Density and power are scalars; all the others are vectors. 

1.2 |a| = 44049 = V13, |b| = VI +0 +1 = V2, a +b = (3,0,2),a-b= 
(1,0,4) and a -b = 2+0 — 3 = -1. To find the angle between a and b, use 
the formula a - b = |a||b|cos@. Substituting the values already obtained, 
this becomes —1 = V13/2cos@, so cos@ ~ —0.196 and 8 = 101°. 

1.3 The component of u in the direction of v is u-v/|v|. For u = (1, 2,2) 
and v = (—6,2,3), u-v = —6 +4 + 6 = 4 and |v| = /36+4+9 = 7, 
so the answer is 4/7. Similarly the component of v in the direction of u is 
u-v/|u| = 4/3. 

1.4 The plane perpendicular to a = (1,1, —1) is r-a = constant, sox+y—z = 
constant. If the plane passes through z = 1, y = 2, z = 1 then the value of 
the constant is 1 + 2—1 = 2. 

1.5 Let two adjacent sides of the rhombus be the vectors a and b (as in Figure 
1.10). Then since the sides are of equal length, |a| = |b|. The diagonals 
of the rhombus are a + b and a — b. Taking the dot product of the two 
diagonals, (a + b) - (a — 6) =a-a+b-a—a-b—b-b= |a|? — |b|? =0, 
so the diagonals are perpendicular. 

1.6 Consider the unit cube, which has sides of length 1 parallel to the coordinate 
axes. Two opposite vertices are the points (0,0, 0) and (1, 1, 1), so one diag- 
onal is the vector (1, 1,1). Another diagonal runs from (0,0,1) to (1, 1,9), 
so another diagonal is the vector (1,1, —1). Both diagonals have magnitude 
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V3 and their dot product is 1. Using the formula a - b = |a||b| cos@ gives 
1 = 3 cosĝ, so 0 = cos™! (1/3) = 70.5°. 

1.7 Choose the origin at one vertex of the triangle and let. two of the sides be 
represented by the vectors a and b. Then the third side is represented by the 
vector b — a. The midpoint of the side opposite the origin has the position 


vector a + (b—a)/2 = (a +b)/2, so the line from the origin to the midpoint 
of the opposite side is given parametrically by r = A(a +b) /2. Similarly the 
line from A to the midpoint of the opposite side is given parametrically by 
r=a+yp(b/2—a). These two lines meet when A(a+b)/2 = a+p(b/2-a). 
In this equation we can equate the terms in a and b. This gives the solution 
u = À = 2/3 so the lines meet at the point with position vector (a + b)/3. 
Similarly, it can be found that the line from B to the midpoint of the 
opposite side meets the line from the origin to the midpoint of the opposite 
side at (a + b)/3, so all three lines pass through this point. 

1.8 The position vector (1,1, 1) lies on the line and a vector along the line is 
(2, 3,5)—(1, 1,1) = (1, 2, 4). The parametric form is therefore r = (1,1,1)+ 
A(1, 2,4). The cross product form is found by taking the cross product of 
this equation with (1, 2,4): r x (1, 2,4) = (1,1,1) x (1, 2,4) = (2, -3, 1). 

1.9 To prove the sine rule, use the fact that the area of the triangle is |a x b|/2. 
Since this area is also equal to |b x c|/2 and |c x a}/2, 


absinC = bcsin A = casin B, 
and dividing by abc gives the sine rule. 
The cosine rule follows from the magnitude of the vector c: 
e = |e}? =c-c = (a—b)-(a— b) = a? +b? — 2a - b =a? +b? — 2abcosC. 
1.10 (a) (i) For any two vectors a and b, a + b is a vector; (ii) a+ (b+ c) = 
(a+ b) +c; (iii) a+ 0 = 0 +a = a; (iv) a+ (—a) = (-a) +a = 0. Hence 
vectors and addition form a group. 
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1.15 


(b) Vectors and the dot product do not form a group since a -b is not a 
vector. 

(c) Vectors and the cross product do not form a group because a x (b xc) # 
(a x b) x c. Also, there is no identity element. 

(a) |a x b|? + (a - b}? = (|a||b| sin 8)? + (\a||b| cos@)* = |a|?|b|?. 

(b) a x (bx (a x b)) = a x ((b- b)a — (b- a)b) = -a -ba x b. 

(c) (a — b)- (b — c) x (c-— a) = (a — 5b) - (bx c-—bxa+ecxa) = 
a-bxc-—b-cxa=0. 

(d) First consider (b x c) x (c x a). This can be written d x (e x a), where 
d = bxc, which using (1.9) is (d-a)ce—(d-c)a = (bx c-a)c since d-c = 0. 
Now take the dot product with (a x b), which gives the result ((a x b)-c)?. 
Take the cross product of the second equation with a. Using (1.9) this 
can be written |a|/?z —a-za =a x b. Now using the first equation this 
simplifies to |a|?a2 — a = a x b, which can be rearranged to solve for z: 
gz = (a xb+a)/|a|?. Geometrically, the two equations given in the question 
represent the equations of a plane and a line. The solution therefore is the 
point at which the line and plane meet. 

Multiply the first equation by b and the second by a, giving r-ab = b and 
r-ba =a. Subtracting these two equations gives r-ab—r-ba=b-a, 
which can be written r x (b x a) = b — a. This is now in the form of the 
equation for a straight line. 

(a) Set a x b = aa + 8b + ye. The coefficient a is found by taking the dot 
product with b x c, and similarly for 8 and y, giving 


a = (a-bb-c—a-clbl?)/a-bxe, 
B (a-cb-a—b-cla|*)/a-bxe, 
y = (lal*|bl’ - (a-b)’)/a-bxe. 


(b) Dot the equation obtained above with c and multiply through by ax b-e: 
(axb-c)*? = (a-bb-c-—a-clb\’)a-c 

+(a-cb-a—b-cla|*)b-c + (|a|?|bl? — (a - b)?)|el? 

|a|?|b|?|e|* — (a - c)°|b]? — (ab)? |e)? — (b - c)? Ja]? 

+2(a-b)(b-c)(a-c). 


(c) Since the faces of the tetrahedron are equilateral triangles, the dot 
product of any two vectors forming the edges is 1/2. The volume V is 
|(a x b-c)|/6, so 36V? = 1-1/4—1/4-1/44 2/8 = 1/2 and V = 1/62 = 
V2/12. 


The rate of change of h is 
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The first term is zero since Or /Ot = v and v xv = 0. The second term is also 
zero since O0v/Ot is the acceleration of the particle, which from Newton’s 
second law is proportional to the force F which is proportional to r. 

1.16 The contour lines are z? — y = constant, i.e. y = z? + constant. These are 
parabolas, shifted by different amounts in the y direction, as shown below. 


EA 


1.17 At the point (0,1), u(x, y) = (xz + y,—z) = (1,0). At any other point on 
the y-axis, z = 0 so the y component of u is zero and the vector field points 
horizontally, to the right if y > 0 and to the left if y < 0. Similarly, on the 
line z + y = 0 the vector field points vertically. By considering a few other 
points the vector field can be sketched, as shown in the following figure. 
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Solutions to Exercises for Chapter 2 


2.1 The parametric form of the curve is z = t, y = t, z = t?, so dr = 
(dx, dy, dz) = (1, 2t, 2t) dt and in terms of the parameter t the vector field 
F is F = (5t*, 2t,t + 2t?). So F - dr = 5t* + 6t? + 4t? and the value of the 
integral is 


1 
[Fare | 5t* + 6t? + 4t° dt = 4. 
C 0 


2. 


o 


First the parametric form of the line joining the points (0,0,0) and (1,1,1) 
must be found. The general form of the equation of a line is r = a+ tu (1.6). 
Choosing t = 0 to correspond to (0, 0,0) gives a = 0 and choosing t = 1 to 
correspond to (1, 1, 1) gives u = (1,1,1), so the equation of the line is r = 
t(1,1,1) and hence dr = dt(1, 1, 1). The vector field F = (5z?, 22, x + 2y) 
can be written in terms of t as (5t?, 2t, 3t), so F - dr = 5t? + 5t dt and the 
integral is 


1 
[Fear | 5t? + 5tdt = 25/6. 
C 0 


Notice that the answer is different from the previous answer, even though 
the same vector field was used and the start and end points of the curve 
C are the same in each case. In other words, for this vector field the value 
of the line integral depends on the path chosen between the two points. 
Hence F is not a conservative vector field. 

2.3 Using z as the parameter to evaluate the integral, dr = (1,e”,e*) dz and 
u = (y’, z, z) = (e*”,z,e7), so u - dr = (2e?? + ze") dz. The integral is 


1 

o: 

[u-ar= | 2e?7 + ze” dr = [e° + ze” —e*], = e? 
C 0 ` 


using integration by parts. 

2.4 The parametric form of the ellipse is z = acosĝ, y = bsinĝ, z =0,0<08< 
2x, so dr = (—a sin ĝ,bcos0,0)dô and r x dr = (0,0,ab) d0. This means 
that the integral f, r x dr only has a component in the z direction and 
the magnitude of the integral is 2rab. Note that this is twice the area of 
the ellipse. 

2.5 Since the surface is z = 0, the normal is n = (0,0,1). Hence u-n =z+y 


and the integral is 
1 p2 1 oe: 
f f| s+udar = f [zy +y*/2], dz 
o Jo 0 
1 


J| nas 
S 
f 2r + 2dr = (x? + 22]; = 3. 


O 
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Note that the y integral was carried out first. However, it is equally possible 
to carry out the z integral first, and this gives the same result. 

This surface integral has six parts, corresponding to the six faces of the 
cube shown below, which must be considered separately. On the surface 


zr = 0, n = (-1,0,0) and u = r = (2,y,z) = (0,y,z) so u:n = 0 and 
there is no contribution to the integral from this surface. By symmetry 
the same result holds for the surfaces y = 0 and z = 0. On the surface 
zr =1,n = (1,0,0) and u -n = z = 1. Since this is a constant, the 
integral over this surface is just the value of the constant multiplied by the 
area of the surface, which gives 1. The same result holds on the surfaces 
y = 1 and z = 1, so the total value of the surface integral, adding the six 
contributions, is0+0+0+1+1+41=3. 

The two curves y = z? and z = y? meet at (0,0) and (1,1). Doing the z 
integral first, the limits are y? < z < /y and 0 < y <1, so 


1 pyy 1 
J unas = f l a? dedy = f (y3/? — y®)/3 dy = 3/35. 
S 0 Jy? 0 


The surface is written parametrically as (x, y, x+y”), so two vectors parallel 
to the surface are (1,0, 1) and (0, 1, 2y). By taking the cross product of these 
two vectors, n dS = (—1, —2y, 1) dz dy. Since this has positive z component 
the direction of n must be changed, so ndS = (1,2y,—1) dz dy and u - 
ndS = z dz dy. In terms of z and y the region of integration is z +y? < 0, 
xz > —1, so doing the z integral first, -1 < z < —y? and —1 < y < 1, and 
the integral is 


1 p—-y? 1 
[[u-nas= [ J zdzdy = f yt/2 — 1/2dy = 1/5 - 1 = —4/5. 
S —-1J-1 -i 
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2.9 The required volume integral is 


1 p2 73 
I dV f f| | ++ dzdyas 
V o Jı Jo 
1 p2 
f f 327 + 3y? + 9 dy dx 
o Ji 


1 D 
f [32y + y? + 9y]? dz 
0 


1 
l 3r? + 16dz = 17. 
0 


2.10 Doing the integrals in the order z, y, x, the range of z between the two 
planes is -z — 1 < z < z + 1 and y ranges from —vV1 — z? to v1 — r? at 
a fixed value of z. Finally, the limits on the outer z integral are —1 and 1. 


The volume is 
p [= om dz dy dx J [Z 2(x + 1) dy dz 
—z-1 
[ 4V1—2? (4+ 1)dzr = 2r. 
—] 


2.11 The edge of the pond is where z = 0, so z? + y? = 1. The limits of the 
volume integral are obtained as follows. Doing the z integral first, at a fixed 
value of z and y, z ranges from 0 to 1 — z? — y?. The limits on z and y are 
the same as in the previous exercise, so the volume V is 


(er =[ TLT” ewt 
Í ‘an 1 — £? — y? dy dz 


[e y°/3] ae dr 
$| 0-2) as = 2/2, 


where the results of Section 2.1.2 have been used to find the last integral. 
The volume of the pond is therefore approximately 1.57 mê. 

In the case of a hemisphere, the volume is 27/3, so the volume of the 
hemisphere is greater by a factor 4/3. 


V 


160 Vector Calculus 


Solutions to Exercises for Chapter 3 


3.1 f = ryz, so Vf = (yz,2z,ry). At the point (1,2,3) this has the value 
(6, 3, 2). The directional derivative in the direction of the vector (1, 1,0) is 
found by taking the dot product of V f and the unit vector in this direction. 
This is (6, 3,2) - (1,1,0)/V2 = 9/vV2. 

3.2 First write the equation of the surface in the form f = constant, so f = 
y — x — z? = 0. A vector normal to the surface is Vf = (—1, 1, —32°). At 
the point (1, 2,1) this is Vf = (—1, 1, —3). To find the unit normal, divide 
by the magnitude, so n = (—1,1,—3)/V/11. The normal pointing in the 
opposite direction, (1, —1,3)//11, is an equally valid answer. 

3.3 ¢ġ =r = |r| = (z? +y? +27)”, so 


an 


1. ; 
ag Te ty tt) aE = 


(x? + y? + z2)!/2 ° 
Similarly, 


Op _ y Da 
ðy (zw? +y? 4 22)1/2" Əz (z? +y? + 22)1/2 


Hence Vé = (x,y, z)(z? + y? + z7)~1/2 which can also be written as r/r 
or as 7, the unit vector in the direction of r. 
Geometrically, the level surfaces ¢ = constant are concentric spheres cen- 


tred at the origin. The vector V¢ points in a direction perpendicular to 
these surfaces, i.e. radially away from the origin. 


f= constant 
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3.4 The sphere and the cylinder intersect when z? +y? +z? = 2 and z? +y? = 1. 


3.0 


Subtracting these two equations gives z? = 1, so the points of intersection 
are the two circles z? + y? = 1, z = +1. 

The angle between the two surfaces is the angle between the normals to the 
surfaces. By taking the gradients of the two surfaces, these normal vectors 
are nı = (22, 2y,2z) and nz = (2x, 2y,0). The angle 0 between the normals 
is found using the dot product, nı -n2 = |n,||n2|cos@. This gives 


Aa? + 4y? = 2y T? + y? + z? 2V1? + y? cosd. 


At a point of intersection this simplifies to 4 = 2\/22cos6, so cos = 
1/\/2 and hence 6 = 45°. Note that this is the same for all the points of 
intersection. 

Vf = (2ry,xz? + 3y? — 1). This is zero when zy = 0, so either z or y 
must be zero. Thus f has maxima, minima or saddle points at the points 
(+1,0), where f = 0; (0,1/V3), where f = —2/3V3 and (0, —1/V3), where 
f = 2/3V3. Since f = y(x? + y? — 1), the contour f = 0 includes the line 
y = 0 and the circle z? + y? = 1. This means that the points at (+1, 0) 
must be saddle points. Putting together all this information, the sketch of 
f and its gradient is as shown below. Lines are contours f = constant and 
arrows are Vf. 


3.6 f=a-:r = az + ayy + a3z, so Vf = (a1,a2,a3) = a. Geometrically, 


f =constant is the equation of a plane (1.3) which is perpendicular to the 
vector a. 


3.7 Proceeding as in Example 3.6, the vector field F = (sin y, z,0) is conser- 


vative if and only if there is a function ¢ satisfying F' = V4, so 
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a a, Ob _ 


—=siny, By L, az 


Oz 
The first equation gives ¢ = zsiny + h(y, z), where h is an arbitrary func- 
tion of y and z. The second equation then becomes 


0. 


zr cosy + AL T 
y ðy TS 
This equation cannot be satisfied for all values of z, since the terms in- 
volving z do not match and h does not depend on z. Hence F is not a 
conservative vector field. 
3.8 If F can be written as V¢ then 
O¢ 2,2 O¢ 2,,2 
— = — Z= — ’ — = 0. 
uleh Ee a 
Integrating the first of these equations (using the substitution z = y tan 8) 
gives @ = tan™!(z/y) + h(y,z) for any function h. From the second and 
third equations it follows that h may be taken to be zero, so F = V¢ where 
$ = tan—*(z/y). 
Now consider the line integral of F around the unit circle z? + y? = 1, 
z = 0, given parametrically by z = cos 0, y = sin ĝ, z = 0,0 < 0 < 2r. 


$ F-dr 
C 


27 
J (sin 0, — cos 8, 0) - (— sin 0, cos @, 0) d8 
0 


2r 
l —1d@ = —2r. 
0 


At first sight this appears to contradict Theorem 3.1, since we have a non- 
conservative vector field which can be written as the gradient of a potential. 
The resolution of the conflict is that both F and ¢ are undefined at the 
origin, so Theorem 3.1 does not apply. 

3.9 The components of V¢ are just the partial derivatives of ¢ with respect 
to x, y and z, so for @ = z? + zy + yz*, Vo = (22 + y, 2 + 27, 2yz). The 
Laplacian can be found either by taking the divergence of this vector, or the 
sum of the second partial derivatives of ¢, giving the result V7 = 2 + 2y. 

3.10 The gradient V¢ is 


Vo = (k cos(kz) sin(ly) exp( vV k? + l?z), 


lsin(kz)cos(ly)exp(v k? + ł?z), 
V k2 +l? sin(kz) sin(ly) exp(V k? + Pz) . 


Now take the divergence of this vector: 
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V*¢ = —k*sin(kz) sin(ly) exp(/k2 + 122) 
—I? sin(kz) sin(ly) exp( v k? + 122z) 
+(k? + l°) sin(kz) sin(ly) exp(V k? + 12z) 
== O 


This means that the function ø is a solution to the equation V7¢ = 0, 
which is Laplace’s equation. 

3.11 The unit normal to the surface @ = constant is V¢/|V@|. Here, ¢ = 
ry? + 2yz so Vd = (y*, 2ry + 2z,2y). At the point (—2,2,3), Vo = 
(4, —2,4) which has magnitude V16 + 4+ 16 = 6, so the unit normal is 
n = (2/3, —1/3, 2/3). 

3.12 For ¢(2z,y,z) = T? +y? +2? + zy —32, Vo = (2x +y — 3, 2y +2, 2z). $ has 
a minimum or maximum where V¢ = 0, which gives 2x + y = 3, 2y = —z, 
z=0,soz=2, y = —1, z = 0. At this point the value of @ is —3. Since 
$ becomes large and positive when z, y or z become large, this must be a 
minimum value. 

3.13 First find a normal to the surface. If f = r? + y? — 2z? then Vf = 
(22, 2y, -6z7) so a normal to the surface at the point (1,1, 1) is (2, 2, —6). 
The equation of the plane is therefore 2z + 2y — 6z = constant, using (1.3). 
Imposing that the plane must pass through (1, 1,1) gives the value of the 
constant to be —2. 

3.14 (a) For u = (y,2,2), 


_ (Ox Oz Oy Ox Oz dOy\ | 
& az’ Oz Be ba SERE 
(b) For v = (xyz, 27, g y), 


_ 3(zyz) | Ə(2*) | Əz -y) 


V-v Dx "On 0 Oe. =yz+0+0= yz. 
vxy = [Atay _ 2) O(eyz) _ Ae—y) 3?) _ (zyz) 
Oy Oz” @z Or” ôr Oy 


= (-1-22z,ry —1,-zz). 


3.15 These results follow directly from the linearity of the differential operator: 


V . (cu + dv) 


0 ð ð 
zz +dv) + gy + dvg) + z, (cus + dv3) 
= goui ðv „ĝu , J22 ou | g? 


Ox Ox Oy Oy Oz Oz 
= cV-u+dV-v, 
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and similarly for curl. 

3.16 u is irrotational if V x u = 0. V x u = (1-1, az — 2z,bcosz — cos z), so 
the solution is a = 2, b= 1. 

3.17 (a) V x u = (—2zsiny + 2yr + 2zsiny — 2ry, y? — y?, 2yz — 2yz) = 0 so 
u is irrotational. 
(b) The potential function ¢ for which u = V@ can be found using the 
step-by-step method of Example 3.6 or as follows: 0¢/0x = y?z suggests 
@ = xzy*z. Taking V¢ gives all the terms in u except the trigonometric 
terms. The z component suggests ¢ = z? cosy + ry?z. Taking the gradient 
of ġ gives all the terms in u, so this (plus an arbitrary constant) is the 
potential function. 
(c) The line integral of u along the curve is just the difference between the 
values of ¢ at the endpoints. The endpoints are (0,0,0) and (1,0, 1), so the 
line integral is #(1,0, 1) — ¢(0,0,0) = 1. 


Solutions to Exercises for Chapter 4 


4.1 cijkajbk + ajdjcj = ei. 

4.2 First, tidy up and rearrange the equation using (4.3) and the symmetry 
property €kji = €ixnj: Ci + Eikjakbj = dpbj}emcmc;. Now this can be written 
as the vector equation c + a x b = (d - b)(e- c)e. 

4.3 [a x b]; = €ijkâjbk = —€;4j0;b, = —cikjbkaj = —[b x a]. 

4.4 (a) ĝij€ijk: this expression is always zero, since if i = j then €;;j = 0, while 
if ¿ Æ j then ĝ;j = 0. 
(b) Eijk Eilm = éjkiĉilm using (4.6). Using (4.12) this is jkm — Ôjm Ôk- 
(c) EijkEijm = ÔjjÓkm = Ôjmôkj» using (b) with | = J. Using 05; = 3 and the 
substitution property of 6;; this can be simplified to Siin = Ôkm = 2ðkm. 
(d) Using (c) with m = k, €i;x€ijn = 2644 = 6, as obtained in Example 4.9. 

4.5 Using suffix notation, 


axb-cxd 


Eijk QjbkEilmCidm 

(Oj10km — Ôjmôkı)ajbkcrdm 
QI bmn Cl dm = ambicidm 
(a-c)(b- d) — (a-d)(b-c). 


4.6 Let C = AB, so in suffix notation C;; = Aj,B,;. Then 
(AB)}, = CE = Cji = Ajr Bri. 
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4.7 


4.8 


Now consider BT AT: 
(BTA’);; = Bi At; = BriAjr- 


These two expressions are the same (recall that ordering of terms does not 
matter in suffix notation) so (AB)! = BT A’. 
The determinant of a 3 x 3 matrix M is 


Mu Mia Miz 
[M| =| Ma Meo. Mzz 
Mz, M32 M33 
Expanding this determinant, 
IMI. = Mii (Me22M33 — Mo3M32) + Mı2(M23M31ı — M21 M33) 
+M13(M2; M32 — M22 M31). (8.37) 


Now expand the suffix notation expression €;i4.M@1;M2;M3,. Since i, j and 
k are repeated, there is a sum over all three indices, so 


3 3 3 


Eijk MiiM2;M3k = 3 > `S Eijk Mii M25; Mx. 


i=1 j=1 k=1 


Since only six of the 27 elements of €;;, are non-zero, there are six terms 
in this sum, and writing them out gives (8.37), so we have shown that 
|M]| = cij Mı: Mo; Me. 

Now turn to the formula (4.10), €pgr|M| = cije MpiMq; Mrk. First note that 
the formula is true for p = 1, q = 2 and r = 3, since in this case it reduces 
to the result shown above. Now consider the effect of interchanging p and 
q. The l.h.s. changes sign, since €pgr = —€gpr. The r.h.s. becomes 


€ijkMgiMpj Mrk = EjikMajMpiMrk (relabelling ¿i + j) 
F. —€ijk Mpi Mg; Mrk, 


so the r.h.s. also changes sign when p and q are interchanged. Similarly, 
both sides change sign when any two of p, q and r are interchanged. This 
suffices to prove the result, since both sides are zero when any two of p, q 
and r are equal and all permutations of 1, 2,3 can be achieved by a suitable 
sequence of interchanges. 

Make use of Epar| M| = Eijk MpiMqj Mrk. 

(a) Multiplying both sides by e,,, and using the result of Example 4.9, 
6|M| = €pgr€ijz MpiMq;M,x. Note that because of the six repeated suffices 
there are 3° = 729 terms in this sum! 

(b) Using the above result, 
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6|MT] = epareij M Mi, May, = €pqr€ige MipMjy Mrr. 

Now since 7, j, k, p, q and r are all dummy suffices, we can relabel z + p, 
j 4 q, k 471, so that the formula for |M7| is identical to that for |M]. 
(c) The formulae for |M| and |N| are €pgr|M| = €:j4 MpiMgjMrk, €pqr|N| = 
Elmn Np Ngm Nrn. By multiplying these together, we obtain 

6| M||N| = €ijk€lmnMpiNptMgj Nam Mrk Non. 
Now MyiNpi = MZ Np = (MTN), so 

6|M||N| = €:jn€imn(M?N)ia(M™N)jm(M™N) en = 6|MT N], 

using the result of part (a). Thus we have shown that |M||N| =|M7N|, or 
equivalently |M7||N| = |MN|. Applying the result of part (b), it follows 
that |M||N| = |M N]. 

4.9 Since each term in the equation is a vector, we first introduce a free suffix 
i for each term: (a x b); + c; = (a- b)b; — d;. Now we introduce dummy 
suffices for the dot and cross product, making sure that i is not reused: 
€4jk 0; Dy + C; = ajb;bi — d;. 

4.10 (a) Using (4.3) twice and (4.4), 


ÔijÓjkÂki = ÔikÔki = Ôi = 3. 
(b) Using (4.12), 
EijkEklimEmni = (Ôuðjm — Ôimôjt)Emni = Ejnl — Einiĝji = Ejnt- 


4.11 Using (4.3), ô:jajbiCkðni = aibick. Here the 7 is a dummy suffix and the k 
is a free suffix, so the result is the k component of the vector (a - b)c. 

4.12 (a) Vx(fVf)= Vix VF+fV x (Vf), using (4.28). Each of these terms 
is zero since any vector crossed with itself gives zero and the combination 
curl grad is always zero. 
(b) V-(fVf) = VF-VE+fV-(VF), from (4.27). This can be simplified 
to V-(fVf) =|VFIP? + fV7F. 

4.13 u is solenoidal if its divergence is zero. 


V-u V -(Vf x Vg) 
(Vx Vf)-Vg-(V x Vg): Vf using (4.29) 
0 using (3.23). 


4.14 Applying (4.35) with u = v, the second term is zero and the fifth and sixth 
terms cancel, leaving 


u-Vu=(V(u-u) —2ux(V x u))/2 
which is (4.34). 
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4.15 (a) In suffix notation, 


ð _. ð Ou; 
° 2 = _ — 2 = passin N 
Mae At pe ee, 
Similarly, ; i 
VV.u= OV:u_ ð Oui 


02,02; ~ O2;02; Az; 
so the two expressions are equal since the order of partial derivatives can 
be interchanged. 
(b) Using (4.24), 


V-V7u 


a (V x u)) 
V-(V(V-u)) (since div curl is zero) 
V7(V-u) (since V? = V-V). 


Note that the first V? acts on a vector but the second acts on a scalar, so 
they must be interpreted differently. 
4.16 Take the divergence of the equation: 


V-utV- -Vxw=V-VGtV-Vu. 


The first term is zero as u is solenoidal. The second term is zero because 
the combination div curl is always zero. The last term is also zero since 
from the previous exercise, V-V2u = V? V -u = 0. So the equation reduces 
to 0 = V . Vo = V?¢ which is Laplace’s equation. 


Of(r) _ df(r) Or _ apti 
[V f(r)]: = Dr; = “ae oe f (r)—, 
using the usual rule for differentiating a function of a function together 
with the result ge = 2;/r from (4.19). Thus V f(r) = f'(r)r/r. 
4.18 u = h(r)r. 
(a) V xu = V x (A(r)r) = VAX r+hV xr = h'(r)r x r/r = 0, using 
the results of the previous exercise for Vh. 
(b) V-u=V-(A(r)r) = VA rt+thV -r = h'(r)jr-r/r+3h=rh'(r)+3h. 
So if V -u = 0, h(r) obeys the differential equation 
dh 


ay +3h= 0. 


(c) Using the method of separation of variables, 


J$- dh _ dr 


which gives log h = —3logr+c= ME +c for some constant c. Taking 
the exponential of both sides, h = A/r3 where the constant A = expe. 


4.17 
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(a) For a Beltrami field, V-u = V-(cV x u) =cV-(V x u) = 0. 

(b) Let v = V x u, so u = cv. Taking the curl of this equation, V x u = 
V x (cv) sov=cV xv. 

(c) If u = (siny, f,g), then u =cV x u gives the three equations 


ape Zeo p20! -e(?i_ 
sny =e (2 3) fel a), g=o( 3 cosy). 


Given that g is independent of z, it follows that f = 0, g = —ccosy and 
sin y = c(csiny). Hence either c= 1, g = — cosy or c= —1, g = cosy. 


Solutions to Exercises for Chapter 5 


o.1 


5.2 


9.3 


The surface integral is equal to the volume integral of V - u, but V-u = 
sin y + 0 — sin y = 0, so the value of the integral is zero. 
u = (y,z,z2 — 1), so V . u = 1. The volume integral is therefore 


1 pl fl 
J| v-uw = | | f l dz dy dz = 1. 
V o Jo Jo 


The surface integral has six parts from the six faces of the cube. On the 
face where z = 0, n = (—1,0,0) and so u : n = —y. Similarly, on the face 
where z = 1, n = (1,0,0) and u-n = y, so the surface integrals from 
these two faces cancel. The same argument holds for the faces y = 0 and 
y = 1. On z = 0, n = (0,0, —1) and u-n = —z + 2 = z, while on z = 1, 
n = (0,0,1) and u-n = z-z = 1- z. The integrals over these two surfaces 


then give Eoi 
f f| z+1-zdzqy=1. 
o Jo 


Therefore, both the surface integral and the volume integral give the answer 
1 so the divergence theorem is verified. 

In order to use the divergence theorem, the volume integral must first be 
written in terms of a divergence. This can be done using (4.27): 


If i Ify ($u) - V -u dV. 


Now since the fluid is incompressible, V - u = 0. Applying the divergence 
theorem then gives 


J| = vow = H gunas. 


Since it is given that u: n = 0 on S, the value of the integral is zero. 
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5.4 This result follows directly from applying the divergence theorem to the 
vector field Vf, since V- Vf = V? f =g. 

5.5 To apply the divergence theorem, a closed surface must be used. Let S’ 
be the surface z = 0, z? + y? < 1, forming the base of the hemisphere. 
The divergence theorem can now be applied over the entire closed surface 


S + S', giving 


JJ, 5 va = [f v- nas+ ff v-nas. 


The surface integral over S can therefore be found by subtracting the sur- 
face integral over S’ from the volume integral. 

For v = (2+ y,2?, 27), V -v = 1, so the volume integral is just the volume 
of the hemisphere, 27/3. The surface integral over S’ is 


J| v-nas= ff -zas 


since n = (0,0,—1) on S’. This integral can be evaluated using polar 
coordinates (r,0), where x = rcosð, 0 <r <1,0 <0 < 2r and dS = 


r dð dr (see Section 2.3.2). 
1 2n 
f f —r? cos? 6 r dô dr 
o Jo 


I -r?° dS 
1 
[ Pra 
0 


=r j4. 


The required integral is then 


I v : ndS = 2r/3 + 7/4 = 117/12. 
S 
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5.6 The argument follows that of Section 5.1.1 with p replaced by q and pu 
replaced by j, so the law of conservation of electric charge is 


5.7 Consider (5.13) in the case where the volume ôV is small. The volume 
integral of Vf is then approximately equal to Vf times the volume: 


VfIV x p fnds. 
6S 


Dividing by the volume and taking the limit dV — 0 gives the definition 


1 
Vi 7 ami wp, tn dS, 
which is analogous in form to the original definitions of div and curl. 
5.8 Using Stokes’s theorem, 


$ r-dr= || vxr-nas=c 
C S 
since V xr = 0. 


5.9 First consider the line integral around the circle z? + y? = 1, evalu- 
ated parametrically using z = cosĝô, y = sin, 0 < @ < 2r, with 
dr = (— sin, cos, 0). The value of the line integral is 


2n 
[ (2cos0 —sin@, — sin? 6,0) - (— sin 0, cos 0,0) dé 
0 
2r , 
= f —2cos6sin@ + sin? 0 — sin? 0 cos 0 d0 
0 


2n 
= i — sin 20 + (1 — cos 20) /2 — sin? 0 cos 0 dô 
0 


= T, 


since all terms except the 1/2 give zero when integrated between 0 and 2r. 
Now to compute the surface integral we need (V x u)- n. The line integral 
was taken in an anticlockwise sense in the z,y plane, so the right-hand 
rule means that n points in the positive z direction, so (V x u)-n = 
(V x u): e3 = 1. The value of the surface integral is then just the area of 
the surface, which is 7 since the surface is a disk of radius 1. Thus Stokes’s 
theorem is verified since the surface integral and the line integral are both 
equal to r. 
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5.10 Applying Stokes’s theorem, 


$, [Vo-ar [fv «#¥9)-n as 


J (Vf x Vot+fV x(Vg)) ndS 
S 


J (Vf x Vg) -n ds. 
S 


Similarly, 


pavf-ar= || (vax VH ndS=- || Fx Vo) -n as. 


Therefore 


$, [Vo-dr=— $ oVf-dr. 


Alternatively, this result can be obtained by applying Stokes’s theorem to 
the line integral of V (fg). 

5.11 If u is irrotational, V x u = 0 and so u x Vf = -V x (fu). Applying 
Stokes’s theorem, the surface integral is equal to ¢, — fu - dr. 

5.12 The rate of change of the total magnetic flux through a surface S is 


$ [| B-nas= |] Y x (ux B)-nds =f ux B-dr. 
dt S S C 


The curve C is a streamline so u is parallel to dr and hence the scalar 
triple product that appears in the line integral is zero. Therefore the flux 
of B through S does not change with time. 

5.13 Applying (5.18) with v = r gives 


I ÔT nj- Gaon dS = [fv x ar] ; 
Ir, d C j 


Using O2;,/02r; = 6;, and 6;; = 3, this simplifies to 


[|mas = [f oxar] . 


For a flat surface nj is constant so the 1.h.s. has the value 2An,, giving the 
required result. 
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Solutions to Exercises for Chapter 6 


6.1 For Cartesian coordinates, u; = z;. The scale factor hı is, from (6.2), 


Or Oz) Ox. 0x3 
Or} 


hı = = |(1,0,0)| = 1. 


62,’ Or,’ O21 


By the same reasoning hz and hz are also 1. 
6.2 (a) Using the definition of the scale factors, 


he = (vw, v(1 —_ w?)1/2, u)| = (u? +v7)/2 
hy = |(uw,u(1 — w?)/?, -v)| = (u? + v2), 
hy = (uv, —uvw(1 — w?)-/2,0)| = uv/(1 — w?)!/2, 


(b) To show that the (u,v, w) system is orthogonal we first need to compute 
the unit vectors. From (6.1) and the above results for the scale factors, these 


are 
eu = (vw, v(1 — w?)/?,u)/(u? + v7), 
€x = (uw,u(l— w?)!/2, —v)/(u? + v?)!/?, 
Cw = ((1 z w*)/?, —W, 0). 


Now take the dot product of these unit vectors to check that they are 
orthogonal: 


€u ` Cy = (uvw? + uv(1 — w?) — uv) /(u? +v?) =0, 
and similarly ey -€» = 0 and ey - €w = 0. 
(c) The volume element in the (u, v, w) system is 
dV = hyhohg du dv dw = (u? + v”)uv/(1 — w?)!/? dududw. 


6.3 Following the approach of the previous exercise, the scale factors are found 
to be hy = (u? + v?)!/?, hy = (u? + v?)!/? and he = uv, so the volume 
element is dV = (u? + v?)uv. The volume V between the surfaces u = 1 
and v = 1 is therefore 


1 sl p2r 1 
V = J fl J (u? + v?)uv dO du dv = an | [u*v/4 + uv /2], dv 
o Jo Jo 0 


which yields the result V = 7/2. 
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6.4 In Cartesian coordinates, V f and e are 


vi= (5 of of) as dz hs = + (# Ax. azn) 


Ax,’ ðT2' T3 ~ Our Ou,’ Ou,’ Ou, 
so the component of V f in the e, direction is 


_ 1 (Of Ox, | Of Oz2 | Of Or3\ _ 1 OF 
eiry a = (5258 Ou, Oro Ou, + nae l= Roe 


Similar expressions for the e3 and eg components follow, giving the formula 
(6.8) for Vf in an orthogonal curvilinear coordinate system. 

6.5 Assume that the cylinder is aligned with the centre of the apple and that 
a < b. The volume integral can be carried out using either cylindrical 
or spherical coordinates. Using cylindrical coordinates, the limits on the 
coordinates R and @ are 0 < R <a,0 < ¢ < 2r. The limits on z are 
determined by the radius of the sphere, z? + R? = b?, so the limits are 


—V/b? — R2 < z < Vb? — R?. The required volume is 
I D " RdġdzdR = 2r f oR — R? dR 
Vo? —R? 0 


= 4n[—(b? - R?)3/?/3)° 
= 4r(b? — (b? — a?)3/?) /3. 


The proportion of the apple removed is this volume divided by 47b?/3, 
which is 1 — (1 — a?/b)3/2. 

6.6 The limits on the coordinates are 1/3 < 0 < 27/3,0< @ < 27,0 <r <a, 
where a is the radius of the Earth. The volume V is 


27/3 p2n pa 
= J f f r° sin 0 dr dọ d8 = 2r [r3/3], [- cos6],74. = 2ra? /3, 
o Jo 


so exactly half of the Earth’s volume is less than 30° away from the Equator. 
6.7 In spherical polar coordinates the divergence of eg is zero, using (6.23). 
The curl, from (6.25), is 


1 Osind, _ lôr _ cotd 1 
rsinĝð 06 ” E oF co 


6.8 The formula u- Vu = V (|u|?/2) — u x (V x u) is used to find u- Vu. 
Since u is a unit vector, its magnitude is constant so V(|w|?/2) = 0. Using 
(6.17), V x u = ez/R, so u x (V x u) = (0,1,0) x (0,0, 1/R) = (1/R, 0,0). 
Therefore u : Vu = —er/R. 


V x e = 
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6.9 Recall that the definition of the Laplacian of a vector field is V?v = VV - 
v- V x V xv. Using the formulae (6.14) and (6.15), the R component of 
VV vis 

a (un 4 Oem , 120% , 202) 
OR\R OR RO 
which can be expanded to give 
lðvr _ vr? , Our _ 1 dug | 1 Pug 4 ves 8v: 
RƏR R OR? R? ð¢ tR ðpðR ` OROz 
Computing the R component of V x V x v using (6.17) gives 
l1 ð Rove _ 1 Gur 7 O*uR } 8v, 
R? ð¢ ðR R? ôg? ðz? @ROz 
Subtracting these two quantities gives the R component of V°v: 
lðvr vR A ur 2 dug m 1 vR , ZR 

RƏR R` ƏR R ðf R? ðf Tz l 

Note that this cannot be the Laplacian of the R component of v, since it 

involves vy. In fact (V?v)pr and V?(vp) are related by 


(V°v)r = 


(V?v)R = V? (vr) e e 


Solutions to Exercises for Chapter 7 


7.1 The definition L;; = e;-e,; states that L;; is the cosine of the angle between 
the e; and ej unit vectors. Referring to Figure 7.1, the cosine of the angle 
between ei and e, is cos ĝ, and the cosine of the angle between e, and ez is 
the same. The cosine of the angle between e} and ez is cos(7/2—9) = sin@, 
and the cosine of the angle between e} and e, is cos(7/2 + 6) = — sinô. 
These results agree with the matrix in (7.3). 

7.2 u is a vector, so u; = L;;u;. To show V -u is a scalar we need to compute 
its value in the dashed frame: 

Ou! ð Ou; Or} 


(Vu) = ax! T Bat Lists) = Di ae Ox!’ 


using the chain rule. Now using (7.8), 


Ou; ðuj _ Ou; 


A A ee 
(Vou)! = Liga Lie TOA Gry Be; 


=V-u. 


7.3 


7.4 


7.9 


7.6 


7.0 
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Thus V -u has the same value in the dashed and undashed frame, so V u 
is a scalar. 

Since a and b are vectors, a; = Li,a, and b; = Ljmbm. Note that the 
suffices are carefully chosen to avoid repetition. Then 


(a;b;)’ = a;b; = LikLjmakbm. 


This agrees with the transformation rule (7.13), so a;b; is a second-rank 
tensor. 

Note that 7;; can be written as T;; = a;bj where a is the position vector 
(11,22) and b = (r2,—2,) which was shown to be a vector in Example 
7.3. Therefore T;; is a tensor by the result of the previous question. Al- 
ternatively, the fact that T;; is a tensor can be confirmed by following the 
method of Example 7.3. 

Since ¢ is a scalar, ¢' = ġ. The transformation rule for T;, is 


Bg! Ip OtmOtn G6 


jk = Oz'Or, TmT Ox; OT, — OLmOEn 


so Tjk is a second-rank tensor. 
If T; is a tensor, Ti; obeys the rule 


T;; = Lig LjmT em. 
Setting j = ti in this formula, 
T;; = LikLimTkm = Okml km = Tkk = Ty. 


So the value of T;; is the same in the dashed and undashed frames, i.e. Ti: 
is a scalar. 
In suffix notation, the divergence theorem (5.1) becomes 


Ou; = Ae 
I, 9x, dV = fp win dS. 


This result holds if u; is replaced by Tij, T2; or T3;, giving the required 
result (7.16). 


7.8 If Qijxe is a tensor of rank four, 


7.9 


Qi jx = LipLjqhirLisQpars- 
Qijxt obeys the above rule, so by setting k = j, 
Qijj = LipLjqLjrLisQpqars = LipdgrLisQ pars = LipLisQprrs; 


so Qij;1 obeys the transformation law for a tensor of rank two. 
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7.10 Given that u;a; is a scalar, 
od d 
u,a; = uja; = ujLija; 


so (u; — u;L,;)a, = 0. If this holds for any a;, then u; = Liju;, so u; is a 
vector. 
7.11 B,s is an anti-symmetric tensor, so 


B; = LruLsıBuv = —LsvLruBvu = -Bir 


Hence B is also anti-symmetric in the dashed frame. 
7.12 If B,, is anti-symmetric, B,, = —B,, so B,, = —B,, and hence B,, = 0. 
7.13 Aijk has the properties Åijk = Ajik and Aijk = — Åikj. Repeatedly apply- 
ing these rules alternately gives 


ijk = Åjik = —Ájki = —AÅkji = Arij = Åikj = —Aije- 


So any element of A;;, is equal to minus itself, hence all elements are zero. 
7.14 (a) Aij = cijk Bk, SO 


f z> l 
Aij = Eijk B, 
= Lin LjnLkpEmnplkr Br 


= LimLjndpr€mnpBr 
= LimLjn€mnpBp = LimLjnAmn, 


so A;; obeys the transformation rule for a second-rank tensor. Since €;;, is 
anti-symmetric with respect to any two indices, A;; is antisymmetric. 
(b) Multiply through by €;jm: 


EijmÅij = €ijmeigk Be = 20emBe = 2Bm, 


making use of Exercise 4.4(c). Hence Bm = €:jm Ai; /2. 

7.15 The most general isotropic fourth-rank tensor is a;ijpı = A6i; 9x1 + HÔikôjı + 
y6i10j;~, from Theorem 7.4. From (4.12), the difference between products of 
6;; terms can be written in terms of a products of €;;, terms. For example, if 
A=1, u = —1, v = Q, we have the isotropic tensor aijki = 6:;0K1 — ÔikÔji = 
€CilmEmyk- 

7.16 Since 6;; and €;;, are isotropic tensors, the combination Aijkim = 9ij€kim 
is an isotropic fifth-rank tensor. Since any two of the five suffices can be 
chosen for the 6, for example Qijkim = Oik Ejim OF Aijkim = Oil€jkm> the total 
number of different tensors of this type is the number of ways of choosing 
two objects from five, which is 5!/2!(5 — 2)! = 10. So there are at least ten 
different components in the most general isotropic fifth-rank tensor. 
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7.17 Following the approach of Section 7.4.2, the kinetic energy of a volume 
element dV is p|v|?dV/2 and v = 2 x r, so 


1/2 fff eer dV 


= 1/2 If P Eijk 257 kEilm SIT m AV 
V 

= 12 [ff P (Qj NjTkrE = (257; Tk) dV 
V 


= 1/2 Jff p (jrr? — Tirk) 23Q dV 
y 


= 154252 /2. 


E 


Solutions to Exercises for Chapter 8 


8.1 The size of the body L is a length and the diffusivity k has units of 
length? /time. The only combination of these which has the units of time is 
L* /k. Therefore the time for heat to diffuse through a body is proportional 
to the square of the size of the body. 
(a) Assume that the mammoth is the same shape as the chicken, made of 
the same material (so the wooliness is ignored) and 20 times the length of 
the chicken. The defrosting time is therefore 20? x 6 hours which is 100 
days. 
(b) Assume that to cook properly, a certain temperature must be reached 
in the interior. The time required for this is proportional to L?. The mass 
M is proportional to L3, so L œ M'/3. Therefore the cooking time should 
be proportional to the two-thirds power of the mass: t x M?/3. 

8.2 Using (8.4), 


Op _ OE 

a OY OE 
= TAA $ (V xB- Hoj) /couo from (8.7) 
= -V-j 


since the combination div curl is always zero. 

8.3 Gauss’s law says that the total flux of electric field through the surface is 
the total charge within the surface divided by €9. Taking the surface to 
be the surface of a sphere of radius r, with area 4rr?, Gauss’s law gives 
Arr? E, = Q/co, which agrees with the result of Example 8.4. 
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8.4 


8.5 


8.6 


8.7 


Vector Calculus 


Take the curl of (8.12): 


VxE 
V x(V x B) = poco 
Expanding the r.h.s. and using (8.10) and (8.11) gives 
ð? B 
-V’?B = -Hoco => oe 


so B obeys exactly the same wave equation as E. 
Given the electric field E = Egf(k-x—wt) = Eof(u) where u = k-z — wt, 
the magnetic field B can be found using (8.11): 


VxE = Vx (Eof(u)) 
= Vf x Eo 


= Fis x Eo. 
Integrating with respect to t and changing the sign to find B gives 
1 1 
B = =k x Eof(k -x - wt) = =k x E. 
wW W 


Hence the magnetic field B is perpendicular to the electric field E. 
The energy can be written w = B-B/2+ E-E/2c’, so the rate of change 
of energy is 


ðw 8B 1_ ðE 

a | P ata” a 
= -B-VxE+E-VxB using (8.11) and (8.12) 
= —V-(E~xB). 


This gives the conservation law 


Ow 

Ot +V-P=0 
where P = Ex B is known as the Poynting vector, representing the energy 
flux of the electromagnetic wave. 
(a) If the material is in equilibrium then the stress tensor must be symmet- 
ric (8.16), so a = 0 and c= b. Also, OP;,;/0z; = 0. ae the first row (i = 1) 
this is satisfied. The second row (i = 2) gives mzz’~' + bz = 0, som = 2 
and b = —2. From the third row, cz3 + nz3™! = Zo. so n = 2 and c = —2. 
(b) For the surface z; = 1 the normal is (1,0,0) (the normal points in 
the direction of the side which is exerting the force). Thus the force is 
F; = Pijn; dS = Pii dS. Since Pii = Pazı = 0 and P3ı = ZT, the force is 
only in the z3 direction. Its magnitude is 


1 pl 
r= | f T2 dT dzz = 1/2. 
0 vO 
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8.8 The r.h.s. of (8.22) is 


1 1 Ov 
~ 5 Stak [V xv] = — gEikEkm a, 


1 OUm 
— 5 (fiim z imit) BE 


— 1/9: _9\_¢ 
J 2 Oz; Oz; ct 


8.9 Applying the condition ôP;;/ðx; = 0 and using (8.21), 


= i OE xr OE;; 
= Abi; Oz; $ au Ox; 
T ð Ov; Ov; 070; 
7 Aan, Oz; i (aes zs seats | 


ð 2 
AEH a V TREY vi. 


8.10 (a) From the definition of the strain tensor, Ey; = az2, Ei? = En = 


8.11 


az, /2 + br,, E22 = —2bxz and the other components of E;; are zero. 

(b) Using (8.21), Pii = A(a = 2b) x2 + 2pare, P,2 = Po = 2u(az, /2+ bz), 
Poo = (a — 2b)z2 — 4ubr2, P33 = A(a — 2b)r2 and other components of P;; 
are zero. 

(c) Applying the equilibrium condition OP;;/Oz; = 0, this is identically 
satisfied for i = 1 and ż = 3, but for i = 2, pa + 2ub + A(a — 2b) — 4ub = Q. 
This is satisfied if a = 2b. 

(a) The term V- (pu) in (5.9) can be written u: V (po +p1)+ (po +p1)V -u. 
Now since po is a constant and terms involving products of u and pı can be 
neglected, this simplifies to po V - u, so (5.9) becomes 0p; /Ot + po V -u = 0. 
(b) In the Navier-Stokes equation (8.28) the product term u- Vu can be 
ignored and there is no body force b or viscosity u, so the remaining terms 
are (po + pi)(Ou/Ot) = —V(po + pi). The term involving pı and u can be 
ignored, po is constant and p, = ap;, so ppdu/dt = -aV pı. 

(c) The velocity u can be eliminated from (8.35) and (8.36) by taking the 
time derivative of (8.35): 


Op, _ Ou 
7 = pov: Ot 


Hence the density perturbation obeys the wave equation (8.13). Physically, 
the waves are sound waves travelling through the fluid. 


=aV - Vp; = aV’ pi. 
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